blob: a8c19f359bba51b51b83334d0bf2dd678098f0a8 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
//Example 2.6 On Velocity and Stream Function
// Initialisation of variables
function[z] = shi(x,y)
z = x^2 - y^2;
endfunction
// Calculations
h = 0.00001;
u = (shi(3,2+h)-shi(3,2))/h; // Partial derivative wrt y
v = -(shi(3+h,2)-shi(3,2))/h;
Velo = sqrt(u^2+v^2);
theta = atand(v/u);
//Results
disp(theta,"Anticlockwise angle at P(3,2) (Degree)", Velo, "Magnitude of velocity at P(3,2) (ft/sec)");
|
