function hp = nadrzdif (t, h) %volani: t-doba simulace, h-pocatecni podminka
q1 = 0.1;    %zadani konstant
A = 0.1;     %zadani konstant
fi = 0.9;    %zadani konstant
g = 9.81;    %zadani konstant
hp = (q1 - A*fi*sqrt(2*g*h))/2; %dif. rovnice
