A FAST, STABLE AND ACCURATE NUMERICAL METHOD FOR THE BLACK–SCHOLES EQUATION OF AMERICAN OPTIONS

In this work we improve the algorithm of Han and Wu [SIAM J. Numer. Anal. 41 (2003), 2081–2095] for American Options with respect to stability, accuracy and order of computational effort. We derive an exact discrete artificial boundary condition (ABC) for the Crank–Nicolson scheme for solving the Black–Scholes equation for the valuation of American options. To ensure stability and to avoid any numerical reflections we derive the ABC on a purely discrete level.Since the exact discrete ABC includes a convolution with respect to time with a weakly decaying kernel, its numerical evaluation becomes very costly for large-time simulations. As a remedy we construct approximate ABCs with a kernel having the form of a finite sum-of-exponentials, which can be evaluated in a very efficient recursion. We prove a simple stability criteria for the approximated artificial boundary conditions.Finally, we illustrate the efficiency and accuracy of the proposed method on several benchmark examples and compare it to previously obtained discretized ABCs of Mayfield and Han and Wu.

[1]  Peter A. Forsyth,et al.  Analysis of the stability of the linear boundary condition for the Black–Scholes equation , 2004 .

[2]  R. C. Merton,et al.  Theory of Rational Option Pricing , 2015, World Scientific Reference on Contingent Claims Analysis in Corporate Finance.

[3]  Matthias Ehrhardt,et al.  Discrete transparent boundary conditions for the Schrödinger equation: fast calculation, approximation, and stability , 2003 .

[4]  Ghada Alobaidi.  American options and their strategies , 2000 .

[5]  Antonella Basso,et al.  A two-step simulation procedure to analyze the exercise features of American options , 2004 .

[6]  Matthias Ehrhardt,et al.  Discrete transparent boundary conditions for parabolic systems , 2006, Math. Comput. Model..

[7]  Matthias Ehrhardt,et al.  Fast calculation of energy and mass preserving solutions of Schrödinger-Poisson systems on unbounded domains , 2006 .

[8]  Peter A. Forsyth,et al.  Quadratic Convergence for Valuing American Options Using a Penalty Method , 2001, SIAM J. Sci. Comput..

[9]  Robert M Corless,et al.  Compact finite difference method for American option pricing , 2007 .

[10]  Christophe Besse,et al.  A Review of Transparent and Artificial Boundary Conditions Techniques for Linear and Nonlinear Schrödinger Equations , 2008 .

[11]  Matthias Ehrhardt.  Nonlinear models in mathematical finance : new research trends in option pricing , 2008 .

[12]  E. Comparini.  Regularization Procedures of Singular Free Boundary Problems in Rotational Bingham Flows , 1997 .

[13]  F. Black,et al.  The Pricing of Options and Corporate Liabilities , 1973, Journal of Political Economy.

[14]  Eduardo S. Schwartz,et al.  The Valuation of American Put Options , 1977 .

[15]  Radha Panini,et al.  Option pricing with Mellin transnforms , 2004, Math. Comput. Model..

[16]  Brian J. McCartin,et al.  Accurate and efficient pricing of vanilla stock options via the Crandall-Douglas scheme , 2003, Appl. Math. Comput..

[17]  Matthias Ehrhardt,et al.  Discrete artificial boundary conditions , 2002 .

[18]  M. Broadie,et al.  American Option Valuation: New Bounds, Approximations, and a Comparison of Existing Methods , 1996 .

[19]  NUMERICAL METHOD FOR THREE-POINT VECTOR DIFFERENCE SCHEMES ON INFINITE INTERVAL , 2008 .

[20]  Muddun Bhuruth,et al.  A fast high-order finite difference algorithm for pricing American options , 2008 .

[21]  R. Geske,et al.  On Valuing American Call Options with the Black-Scholes European Formula , 1984 .

[22]  A. Erdélyi,et al.  Tables of integral transforms , 1955 .

[23]  Houde Han,et al.  A class of artificial boundary conditions for heat equation in unbounded domains , 2002 .

[24]  C. Cryer.  The Solution of a Quadratic Programming Problem Using Systematic Overrelaxation , 1971 .

[25]  Matthias Ehrhardt,et al.  Discrete transparent boundary conditions for the Schrödinger equation , 2001 .

[26]  Dawn Hunter,et al.  Penalty and front-fixing methods for the numerical solution of American option problems , 2002 .

[27]  Muddun Bhuruth,et al.  Numerical pricing of options using high-order compact finite difference schemes , 2008 .

[28]  Curt Randall,et al.  Pricing Financial Instruments: The Finite Difference Method , 2000 .

[29]  P. Wilmott,et al.  The Mathematics of Financial Derivatives: Preface , 1995 .

[30]  Xiaonan Wu,et al.  A Fast Numerical Method for the Black-Scholes Equation of American Options , 2003, SIAM J. Numer. Anal..

[31]  권용훈.  Finite Difference Bisection Algorithm for Free Boundaries of American Call and Put Options , 2003 .

[32]  Juan Carlos Cortés,et al.  A new direct method for solving the Black-Scholes equation , 2005, Appl. Math. Lett..

[33]  Marti G. Subrahmanyam,et al.  Pricing and Hedging American Options: A Recursive Integration Method , 1995 .

[34]  Thomas A. Manteuffel,et al.  The numerical solution of second-order boundary value problems on nonuniform meshes , 1986 .

[35]  Frank Cuypers.  Tools for Computational Finance , 2003 .

[36]  RAUL KANGRO,et al.  Far Field Boundary Conditions for Black-Scholes Equations , 2000, SIAM J. Numer. Anal..

[37]  Xiaonan Wu,et al.  Convergence of difference scheme for heat equation in unbounded domains using artificial boundary conditions , 2004 .

[38]  Walter Allegretto,et al.  Finite Element Error Estimates for a Nonlocal Problem in American Option Valuation , 2001, SIAM J. Numer. Anal..

[39]  Daniel Sevcovic,et al.  Analysis of the free boundary for the pricing of an American call option , 2001, European Journal of Applied Mathematics.