| (%i23) | kill(all); D1:40*10^-3,numer; l1:3; Lamda1:0.055; Qtr1:0; Qk1:0.0003; n:1.0; Hp:3.5; |
| (%i8) | _jg_SQH_TV_z1_1:filename; |
| (%i9) |
Qp1:Qk1-Qtr1; Qr1:Qtr1+(Qp1/2); m1:Qtr1/Qr1; A1:m1^2+((4*m1*(1-m1))/(n+1))+((4*((1-m1)^2))/(2*n+1)); |
| (%i13) | w1:(3.14*(D1)^2)/4; |
| (%i14) | DeltaH:(8*(1-m1)+Lamda1*(l1/D1)*A1)*((Qr1^2)/(2*9.81*(w1^2))); |
| (%i15) | fH(l):=(8*(1-m1)+Lamda1*(l/D1)*A1)*((Qr1^2)/(2*9.81*(w1^2))); fH(3); |
| (%i17) | wxplot2d([Hp-fH(l)], [l,0,l1], [y,3.48,3.5])$ |
