function v=hh(x) % option price as a function of a standard normal random variable x T=.25;K=10;S0=10;r=.05;sig=.3; ST=S0*exp(r*T-.5*sig^2*T+sig*sqrt(T)*x); v=exp(-r*T)*max(0,ST-K)