close all
clear all

cv = [.5:.5:2.0];

SG = 1;

A = pi*(2.5^2)/4;  %in^2

V = 0:.1:8;  %in/sec

Q = A*V;  %in^3/sec

q = 60*Q/231;  %gal/min

figure(1)
hold all
figure(2)
hold all
figure(3)
hold all

for i = 1:length(cv)

  dP = ((q/cv(i)).^2)/SG;
  F = A*dP;
  
  figure(1)
  plot(V,F)

  figure(2)
  plot(V,(F./V)*(4.45*39.4))

  figure(3)
  plot(V,Q.*dP*4.45/39.4);
  
end