P=read('P.dat',-1,1); //- Prix Black et Scholes;
Pc=read('Pc.dat',-1,1); //- Prix avec Cout de Transaction;

xbasc();
printf('Payoff\n');
plot2d(s,P0(s),1,rect=RECT);

scanf('%c');

xbasc();
printf('comparaison des valeurs\n');
plot2d(s,P,2,rect=RECT);
plot2d(s,Pc,5,rect=RECT);

scanf('%c');

xbasc();
printf('comparaison des delta\n');
J=diag(sparse(ones(I,1)),1);
D=(J-speye(J))/h;
//plot2d(s,D*Pbs,2,rect=RECT);
plot(s,P0(s),'black--');
plot(s,D*P,'b--');
plot(s,D*Pc,'r');
//xtitle('Deltas:');
legend('Prix','Delta sans cout','Delta avec cout');