function dh_dt=modelo(t,h)

F_in=400; %m3/h
area=10; %m2
area_sal=0.1; %m2
dh_dt=(F_in/3600-sqrt(20*h)*area_sal)/area; %m/s

    