clc //initialisation of variables p=0.1//m p1=50//m a=1.000//m^1/3/s R=501.45*10^-4//m^2 s=8.73*10^-2//m n=0.030 m=1.672*(0.087)^2/3//m //CALCULATIONS S=p/p1 Q=((1)^1/3)*m*(S)^1/2//m^3 //RESULTS printf('The flow of rate trough the pipe=% f m^3',Q)