##C(S,t) = C(S(t), E, T-t, r) = SN(d_{1}) - Ee^{-r(T-t)}N(d_{2}),##
where,
##N(d) = \frac{1}{\sqrt{2\pi}}\int^{d}_{-\infty}e^{{-\frac{1}{2}}s^{2}}ds,##
is the cumulative distribution function for the standard normal distribution,
##d_{1} = \frac{\ln\frac{S}{E} + (r +...