close all
clear all

cv = 1;  %[1.25:.25:1.75];

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);
  
end