summaryrefslogtreecommitdiff
path: root/1268/CH3/EX3.2/3_2.sce
blob: 3212ffdbd17f72a6c99dc7100e87e2c904cd902e (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
clc;
disp("Example 3.2")
// the formula used is u=2U(1-(r/R)^2)
// In the first part u=U/2 and in second part u=U
// first step
//(r/R)^2=3/4
R=5; // in cm
r1=5*((0.75)^0.5);
// second step
// (r/R)^2=1/2
r2=5*((0.5)^0.5);
disp(" At r= ")
disp(r1)
disp(" cm we have half the avergae velocity and at r= ")
disp(r2)
disp(" we have axial velocity equal to avergae velocity.")