//Example 20.1
clc
Cv=4;
G=1.26;
P=100;//psi
q=Cv*sqrt(P/G);
disp('gpm',q,'q=')