Tihomir Nikolov14th May 2016
BS_call <- function(St,K,r,sigma,t){ d1=(1/(sigma*sqrt(t))*(log(St/K) + (r + (sigma^2)/2)*(t))) d2=(1/(sigma*sqrt(t))*(log(St/K) + (r - (sigma^2)/2)*(t))) St*pnorm(d1) - K*exp(-r*(t))*pnorm(d2) }