Pseudospectral methods for pricing options

Models with two or more risk sources have been widely applied in option pricing in order to capture volatility smiles and skews. However, the computational cost of implementing these models can be large—especially for American-style options. This paper illustrates how numerical techniques called ‘pseudospectral’ methods can be used to solve the partial differential and partial integro-differential equations that apply to these multifactor models. The method offers significant advantages over finite-difference and Monte Carlo simulation schemes in terms of accuracy and computational cost.