Two-Dimensional Fourier Cosine Series Expansion Method for Pricing Financial Options

Abstract

The COS method for pricing European and Bermudan options with one underlying asset was developed in [F. Fang and C. W. Oosterlee, SIAM J. Sci. Comput., 31 (2008), pp. 826--848] and [F. Fang and C. W. Oosterlee, Numer. Math., 114 (2009), pp. 27--62]. In this paper, we extend the method to higher dimensions, with a multidimensional asset price process. The algorithm can be applied to, for example, pricing two-color rainbow options but also to pricing under the popular Heston stochastic volatility model. For smooth density functions, the resulting method converges exponentially in the number of terms in the Fourier cosine series summations; otherwise we achieve algebraic convergence. The use of an FFT algorithm, for asset prices modeled by Lévy processes, makes the algorithm highly efficient. We perform extensive numerical experiments.