We describe an improvement of Han and Wu's algorithm [H. Han, X.Wu, A fast numerical method for the Black-Scholes equation of American options, SIAM J. Numer. Anal. 41 (6) (2003Anal. 41 (6) ( ) 2081Anal. 41 (6) ( -2095] ] for American options. A high-order optimal compact scheme is used to discretis
Numerical pricing of options using high-order compact finite difference schemes
β Scribed by D.Y. Tangman; A. Gopaul; M. Bhuruth
- Publisher
- Elsevier Science
- Year
- 2008
- Tongue
- English
- Weight
- 262 KB
- Volume
- 218
- Category
- Article
- ISSN
- 0377-0427
No coin nor oath required. For personal study only.
β¦ Synopsis
We consider high-order compact (HOC) schemes for quasilinear parabolic partial differential equations to discretise the Black-Scholes PDE for the numerical pricing of European and American options. We show that for the heat equation with smooth initial conditions, the HOC schemes attain clear fourth-order convergence but fail if non-smooth payoff conditions are used. To restore the fourth-order convergence, we use a grid stretching that concentrates grid nodes at the strike price for European options. For an American option, an efficient procedure is also described to compute the option price, Greeks and the optimal exercise curve. Comparisons with a fourth-order non-compact scheme are also done. However, fourth-order convergence is not experienced with this strategy. To improve the convergence rate for American options, we discuss the use of a front-fixing transformation with the HOC scheme. We also show that the HOC scheme with grid stretching along the asset price dimension gives accurate numerical solutions for European options under stochastic volatility.
π SIMILAR VOLUMES
Direct numerical simulation (DNS) has become a powerful tool in studying fundamental phenomena of laminar-turbulent transition of high-speed boundary layers. Previous DNS studies of supersonic and hypersonic boundary layer transition have been limited to perfect-gas flow over flat-plate boundary lay