How can i solve and plot this problem

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,

Answer 1

[Mario Malic]

Use function fplot

fplot(f,[xmin xmax]); % plots f over the interval [xmin xmax].

[esat gulhan]

.