function p=optionprice(K,mu,sig)
% price of a European option when the stock has N(mu,sig) distribution
% K may be a vector
p=(mu-K).*(1-normcdf((K-mu)/sig))+sig*normpdf((K-mu)/sig);