syms D Do g H(t) t
Do=10
D=0.1
g=9.81;
eq=diff(H,t)==pi()*D^2/4*(2*g*H)^0.5/(pi()*Do^2/4);
cond=H(0)==2;
h(t)=dsolve(eq,cond)
When i try to solve this equation, my output is below
h(t) =
(981*(t/10000 + (20*109^(1/2))/327)^2)/200
(981*(t/10000 - (20*109^(1/2))/327)^2)/200
İ only want to get (981*(t/10000 + (20*109^(1/2))/327)^2)/200 value, i dont want second value( (981*(t/10000 - (20*109^(1/2))/327)^2)/200) and i want to plot in 0-100. How can i ?
I tried t=linspace(0,10,100)
plot(t,h)
it gives error,