summaryrefslogtreecommitdiff
path: root/1187/CH10/EX10.1/1.sce
blob: 66375d01c7f2e0a2344e996f2d6db1ccfcadfc6c (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
clc

Q=400; // m^3/s
b2=20; // m
g=9.81; // m/s^2
b1=25; // m

h2=(Q/b2/sqrt(g))^(2/3);
// Since energy is conserved
// h1 + u1^2/2g = h2 +u2^2/2g = h2 + h2/2 = 3h2/2

// h1 + 1/2*g*(Q/(b1h1))^2 = 3*h2/2;

// h1^3-5.16*h1^2+13.05 = 0;

// By solving this cubic equation

h1=4.52; // m

disp("(a) The depth of the water under the brigde =")
disp(h2)
disp("m")

disp("(b) the depth of water upstream =")
disp(h1)
disp("m")