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 1010.

Share and Cite:

Akpan, I. and Fatokun, J. (2015) An Accurate Numerical Integrator for the Solution of Black Scholes Financial Model Equation. American Journal of Computational Mathematics, 5, 283-290. doi: 10.4236/ajcm.2015.53026.

Conflicts of Interest

The authors declare no conflicts of interest.


[1] Dewynne, J., Howison, S. and Wilmott, P. (1995) Option Pricing: Mathematical Models and Computation. Financial Press, Oxford.
[2] Black, F. and Scholes, M. (1973) The Pricing of Options and Corporate Liabilities. Journal of Political Economy, 81, 637-659.
[3] Kangro, R. (2011) Computational Finance.
[4] Bensaid, B., Lesne, J., Pages, H. and Scheinkman, J. (1992) Derivative Asset Pricing with Transaction Costs. Math. Finance, 2, 63-86.
[5] Merton, R.C. (1973) Theory of Rational Option Pricing. Bell Journal of Economics and Management Science, 4, 141-183.
[6] Boyle, P. and Vorst, T. (1973) Option Replication in Discrete Time with Transaction Costs. Journal of Finance, 47, 271-293.
[7] Davis, M., Panis, V. and Zariphopoulou, T. (1993) European Option Pricing with Transaction Fees. SIAM Journal on Control and Optimization, 31, 470-493.
[8] Frey, R. (1998) Pefect Option Hedging for a Large Trader. Finance and Stochastics, 2, 115-141.
[9] Frey, R. (2000) Market Illiquidity as a Source of Model Risk in Dynamic Hedging. In: Gibson R., Ed., Model Risk, RISK Publications, London.
[10] Genotte, G. and Leland, H. (1990) Market Liquidity, Hedging and Crashes. American Economic Review, 80. 999-1020.
[11] Jarrow, R. (1992) Market Manipulation, Bubbles, Corners and Short Squeezes. Journal of Financial and Quantitative Analysis, 27, 311-336.
[12] Platen, E. and Schweizer, M. (1998) On Feedback Effects from Hedging Derivatives. Mathematical Finance, 8, 67-84.
[13] Schonbucher, P. and Wilmot, P. (2000) The Feedback Effect of Hedging in Illiquid Markets. The SIAM Journal on Applied Mathematics, 61, 232-272.
[14] Hodges, S.D. and Neuberger, A. (1989) Optimal Replication of Contingent Claims under Transaction Costs. Review of Futures Markets, 8, 222-239.
[15] Whalley, A.E. and Wilmot, P. (1997) An Asymptotic Analysis of an Optimal Hedging Model for Option Pricing with Transaction Costs. Mathematical Finance, 7, 307-324.
[16] Bohner, M. and Zheng, Y. (2009) On Analytical Solutions of the Black-Scholes Equation. Applied Mathematics Letters, 22, 309-313.
[17] Emamirad, H., Goldstein, G.R. and Goldstein, J.A. (2012) Chaotic Solution for the Black-Scholes Equation. Proceedings of American Mathematical Society, 140, 2043-2052.
[18] Lamberton, D. (1998) Error Estimates for the Binomial Approximation of American Put Options. The Annals of Applied Probability, 8, 206-233.
[19] Lai, Y. and Spainer, J. (1998) Applications of Monte Carlo/Quasi-Monte Carlo Methods in Finance: Option Pricing. Springer, Berlin.
[20] Pironneau, O. and Hecht, F. (2000) Mesh Adaption for Black and Scholes Equations. East-West Journal of Numerical Mathematics, 8, 25-35.
[21] Paras, A. and Avellanneda, M. (1994) Dynamic Hedging Portfolios for Derivative Securities in the Presence of Large Transaction Costs. Applied Mathematical Finance, 1, 165-193.
[22] Goldstein, J.A., Mininni, R.M. and Romanelli, S. (2008) A New Explicit Formula for the Solution of the Balck-Merton-Scholes Equation. In: Sengupta, A. and Sundar, P., Eds., Infinite Dimensional Stochastic Analysis, World Scientific, Singapore, 226-235.
[23] Strikwerda, J.C. (1989) Finite Difference Schemes and Partial Differential Equations (Wadsworth & Brooks/Cole Mathematics Series). Wadsworth & Brooks/Cole Advanced Books & Software, Pacific Grove.
[24] Schiesser, W.E. (1991) The Numerical Method of Lines: Integration of Partial Differential Equations. Academic Press, San Diego.
[25] Haq, S., Hussain, A. and Uddin, M. (2011) RBFs Meshless Method of Lines for the Numerical Solution of Time-Dependent Nonlinear Coupled Partial Differential Equations. Applied Mathematics, 2, 414-423.
[26] Lambert, J.D. (1993) Numerical Methods for Ordinary Differential Systems. John Wiley, Chichester.
[27] Coppex, F. (2009) Solving Black Scholes Equation: A Demystification.
[28] Ankudinova, J. and Ehrhardt, M. (2007) On the Numerical Solution of Nonlinear Black Scholes Equations.
[29] Fatokun, J.O. and Akpan, I.P. (2013) L-Stable Implicit Trapezoidal-Like Integrators for the Solution of Parabolic Partial Differential Equations on Manifolds. African Journal of Mathematics and Computer Science Research, 6, 183-190.
[30] Fatokun, J.O. and Akpan, I.P. (2014) A Trapezoidal-Like Integrator for the Numerical Solution of One-Dimensional Time Dependent Schrodinger Equation. American Journal of Computational Mathematics, 4, 271-279.

Copyright © 2022 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.