AJCM  Vol.5 No.3 , September 2015
An Accurate Numerical Integrator for the Solution of Black Scholes Financial Model Equation
Abstract: 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.
Cite this paper: 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.

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