IMA Journal of Numerical Analysis Advance Access published online on July 1, 2009
IMA Journal of Numerical Analysis, doi:10.1093/imanum/drp016
On convergence of a fitted finite-volume method for the valuation of options on assets with stochastic volatilities

Department of Applied Mathematics, National Sun Yat-sen University, No. 70, Lianhai Road, Gushan District, Kaohsiung City 804, Taiwan
Department of Mathematics and Statistics, The University of Western Australia, 35 Stirling Highway, Crawley WA6009, Australia
Email: hungch{at}math.nsysu.edu.tw
Received on 3 July 2008. Revised on 23 March 2009.
| Abstract |
|---|
In this paper we present a convergence analysis for a fitted finite-volume discretization method for the two-dimensional Black–Scholes equation arising in the Hull–White model for pricing European options with stochastic volatility. We first formulate the finite-volume method as a nonconforming Petrov–Galerkin finite–element method with each basis function of the trial space being determined by a set of two-point boundary-value problems defined on element edges. We then show that the bilinear form of the finite-element method is coercive and continuous. Finally we establish an upper bound of order
(h) on the discretization error of method, where h denotes the mesh parameter of the discretization.
Key Words: stochastic volatility; European option pricing; finite-volume method; stability and convergence