clc;
disp("Example A.15")
dratio=4/5
pratio=dratio^-5
disp(pratio,"Pressure drop is ")