function plotcalloptlogist(r,sigma,T)
% use for plotting either calloptlogist or calloptloggam. For the latter, remove line 1.
%u=rand(1,1000000);   % rremove for loggam
y=gamrnd(2,1,1,200000);
p=[];
for ex=.75:.05:1.4
   [v,x,CALL]=calloptloggam(y,1,ex,.05,.3,.25);

   %[v,x,CALL]=calloptlogist(u,1,ex,.05,.3,.25);
p=[p;[v CALL]];
end
%plot((.75:.05:1.4),p(:,1),'r',(.75:.05:1.4),p(:,2),'g')
plot((.75:.05:1.4),p(:,2)-p(:,1),'r')
xlabel('strike price/initial price')
%ylabel('Black Scholes -Loglogistic Price')
ylabel('Black Scholes -Loggamma Price')