close all
clear all

Ap = 1.37;

Pc = 2000;  %psi
m = 50/400;  %GPMperpsi.


dP = 0:3000; %(psi)

S = 10;  %Ram speed (in/s)

for i = 1:length(dP)
   Qv(i) =  ReliefValveFlowRate(dP(i),Pc,m);  %GPM
end
Qv = Qv*231/60;  %in^3/s;

Qt = (Ap*S);


Qm = Qt-Qv;


SpeedPercentage = Qm./Qt;
idx = find(SpeedPercentage < 0);
SpeedPercentage(idx) = 0;

plot(dP,SpeedPercentage,'x');