close all
clear all

A = pi*2.5^2/4-pi*1.75^2/4;
F = 125:3000;

gamma = 1.0;


D = [58 231/2 231];
P = F/A;
I = [P(1) 200 500];

gamma = 1.4;
Vol_tot_adi = zeros(1,length(P));
for i = 1:length(D)
Vol_tot_adi = Vol_tot_adi+Vol(P,I(i),D(i),gamma);
end

gamma = 1.0;
Vol_tot_iso = zeros(1,length(P));
for i = 1:length(D)
Vol_tot_iso = Vol_tot_iso+Vol(P,I(i),D(i),gamma);
end


figure
hold all
plot(Vol_tot_iso/A,F,Vol_tot_iso/A,F(1)+Vol_tot_iso*(F(end)-F(1))/(Vol_tot_iso(end)-Vol_tot_iso(1)))
plot(Vol_tot_adi/A,F,Vol_tot_adi/A,F(1)+Vol_tot_adi*(F(end)-F(1))/(Vol_tot_adi(end)-Vol_tot_adi(1)))
xlabel('Stroke(in)')
ylabel('Force (lbs)')
P(1)
aa = axis;
aa(2) = 72;
axis(aa);