blob: 20eaab5abec2bf244895a302169bfe5a74ca3340 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
// determine discharge and head with two identail
clc
f=0.025
l=70
D=0.3
k=2.5
g=9.8
m=((f*l/D)+k)/(2*g*(((%pi*D*D)/4)^2))
disp(m)
mprintf('\n H1 =15 +%d Q^2',m)
b=5.35
a=112.8
c=7.9
Q=(1/(2*a))*(b+sqrt((b^2)+(4*a*c)))
mprintf('\n Q= %f m3/s',Q)
H1=15+85*Q^2
mprintf('\n H1 = %f m',H1)
|
