A robust and accurate finite difference method for a generalized Black-Scholes equation. (English) Zbl 1214.91130

Summary: We present a numerical method for a generalized Black-Scholes equation, which is used for option pricing. The method is based on a central difference spatial discretization on a piecewise uniform mesh and an implicit time stepping technique. Our scheme is stable for arbitrary volatility and arbitrary interest rate, and is second-order convergent with respect to the spatial variable. Furthermore, the present paper efficiently treats the singularities of the non-smooth payoff function. Numerical results support the theoretical results.


91G60 Numerical methods (including Monte Carlo methods)
65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs
65M12 Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs
65M50 Mesh generation, refinement, and adaptive methods for the numerical solution of initial value and initial-boundary value problems involving PDEs
Full Text: DOI


[1] Black, F.; Scholes, M.S., The pricing of options and corporate liabilities, J. political econ., 81, 637-654, (1973) · Zbl 1092.91524
[2] Cox, J.C.; Ross, S.; Rubinstein, M., Option pricing: a simplified approach, J. fin. econ., 7, 229-264, (1979) · Zbl 1131.91333
[3] Hull, J.C.; White, A., The use of control variate technique in option pricing, J. fin. econ. quant. anal., 23, 237-251, (1988)
[4] Rogers, L.C.G.; Talay, D., Numercial methods in finance, (1997), Cambridge University Press Cambridge · Zbl 0867.00036
[5] Schwartz, E., The valuation of warrants: implementing a new approach, J. fin. econ., 13, 79-93, (1997)
[6] Courtadon, G., A more accurate finite difference approximation for the valuation of options, J. fin. econ. quant. anal., 17, 697-703, (1982)
[7] Wilmott, P.; Dewynne, J.; Howison, S., Option pricing: mathematical models and computation, (1993), Oxford Financial Press Oxford · Zbl 0797.60051
[8] Vazquez, C., An upwind numerical approach for an American and European option pricing model, Appl. math. comput., 97, 273-286, (1998) · Zbl 0937.91053
[9] Ikonen, S.; Toivanen, J., Operator splitting methods for American option pricing, Appl. math. lett., 17, 809-814, (2004) · Zbl 1063.65081
[10] R. Zvan, P.A. Forsyth, K.R. Vetzal, A general finite element approach for PDE option pricing models, University of Waterloo, Canada, 1998. · Zbl 0945.65005
[11] Forsyth, P.A.; Vetzal, K.R., Quadratic convergence for valuing American options using a penalty method, SIAM J. sci. comput., 23, 6, 2095-2122, (2002) · Zbl 1020.91017
[12] Zvan, R.; Forsyth, P.A.; Vetzal, K.R., Penalty methods for American options with stochastic volatility, J. comput. appl. math., 91, 199-218, (1998) · Zbl 0945.65005
[13] Angermann, L.; Wang, S., Convergence of a fitted finite volume method for the penalized black – scholes equation governing European and American option pricing, Numer. math., 106, 1-40, (2007) · Zbl 1131.65301
[14] Wang, S., A novel fitted finite volume method for the black – scholes equation governing option pricing, IAM J. numer. anal., 24, 699-720, (2004) · Zbl 1147.91332
[15] Cho, C.; Kim, T.; Kwon, Y., Estimation of local volatilities in a generalized black – scholes model, Appl. math. comput., 162, 1135-1149, (2005) · Zbl 1079.91022
[16] O.A. Ladyzenskaja, V.A. Solonnikov, N.N. Ural’ceva, Linear and quasilinear equations of parabolic type, in: Amer. Math. Soc. Transl., vol. 23, Providence, RI, 1968.
[17] Kangro, R.; Nicolaides, R., Far field boundary conditions for black – scholes equations, SIAM J. numer. anal., 38, 1357-1368, (2000) · Zbl 0990.35013
[18] Kellogg, R.B.; Tsan, A., Analysis of some difference approximations for a singular perturbation problem without turning points, Math. comp., 32, 1025-1039, (1978) · Zbl 0418.65040
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.