function y = put(s,sigma,r,K,t)
            % Black-Scholes formel for put
            y = K*exp(-r*t) + call(s,sigma,r,K,t) - s;