function Pcall=Pcall(S,K,T,r,q,cc)

d1=(log(S/K)+(r-q+cc^2/2)*T)/(cc*sqrt(T));
d2=d1-cc*sqrt(T);

Pcall=S*exp(-q*T)*N(d1)-K*exp(-r*T)*N(d2);

end