Untitled

R = 1;       
r = 0.05;   
h0 = 2;      
C = 0.65;    
g = 9.81;    


A_cilindro = pi * R^2;
A_foro = pi * r^2;

% Soluzione analitica
t_svuotamento_analytical = (2 * A_cilindro / (A_foro * C * sqrt(2 * g))) * sqrt(h0)

% Soluzione numerica con il metodo di Eulero
dt = 0.01;  
h = h0;     
t = 0;      
while h > 0
    v_out = C * sqrt(2 * g * h);  
    dh_dt = - (A_foro / A_cilindro) * v_out;
    h = h + dh_dt * dt;
    t = t + dt;
end

t_svuotamento_num = t  % Tempo numerico

% Differenza tra le due soluzioni
t_diff = abs(t_svuotamento_num - t_svuotamento_analytical)