blob: 9190a6ca6e8ac1a4042ba1a1c164b4c8b02388d7 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
clc
//initialisation of variables
clear
A= 1/16 //mile^2
d= 2 //ft
h= 18 //ft
h1= 5 //ft
f= 0.006
l= 200 //ft
h2= 10 //ft
g= 32.2 //ft/sec^2
//CALCULATIONS
X= sqrt(1/((1.5+(4*f*l/d))/(2*g)))
function [y]=fun(H)
y=A*5280^2*H^-0.5/(%pi*d^2*X/4)
endfunction
vec2=intg(h-h1,h,fun)
T= vec2
//RESULTS
printf ('time for the channel to fall = %.f sec',T)
|