An Accurate Numerical Integrator for the Solution of Black Scholes Financial Model Equation

In this paper the Black Scholes differential equation is transformed into a parabolic heat equation by appropriate change in variables. The transformed equation is semi-discretized by the Method of Lines (MOL). The evolving system of ordinary differential equations (ODEs) is integrated numerically by an L-stable trapezoidal-like integrator. Results show accuracy of relative maximum error of order 10–10.