function y = call(s,sigma,r,K,t)
            % Black-Scholes formel for call
            y = s.*normcdf( d1(s,sigma,r,K,t) ) - exp(-r*t)*K*normcdf( d2(s,sigma,r,K,t) );