//example 8.1
//page 288
clc; funcprot(0);
//initialisation of variable
L=10;//length
R=0.002;//radius
P1=200;
P2=260;
pi=3.14;
rho=0.81*1000;
mu=19.1/10000;
Sp=-(P1+0.81*6*9.81-P2)/L;//slope of pressue
taumax=R/2*Sp*1000;
disp(taumax,"maximum shear force(N/m^2)");
vmax=R^2/4/mu*Sp*1000;
disp(vmax,"maximum velocity(m/s)");
Q=pi*R^4/8/mu*Sp*1000;
disp(Q,"flow rate (m^3/s)");
R=vmax/2/mu*2*R*rho;
disp(R,"reynolds no.");
clear