Merton model; wavelet-Galerkin method; integro-differential equation; spline wavelets; Crank-Nicolson scheme; sparse matrix; option pricing
This paper examines the pricing of two-asset European options under the Merton model represented by a nonstationary integro-differential equation with two state variables. For its numerical solution, the wavelet-Galerkin method combined with the Crank-Nicolson scheme is used. A drawback of most classical methods is the full structure of discretization matrices. In comparison, the wavelet method enables the approximation of discretization matrices with sparse matrices. Sparsity is essential for the efficient application of iterative methods in solving the resulting systems and the efficient computation of the matrices arising from the discretization of integral terms. To illustrate the efficiency of the method, we provide the results of numerical experiments concerning a European option on the maximum of two assets.