Option Pricing by the Legendre Wavelets Method

This paper presents the numerical solution of the Black-Scholes partial differential equation (PDE) for the evaluation of European call and put options. The proposed method is based on the finite difference and Legendre wavelets aproximation scheme. We derive a matrix structure for the Legendre wavelets integral operator which has been widely used so far. Moreover, in order to use the payoff function, another operational matrix is derived. By the proposed combined method, the solving Black-Scholes PDE problem reduces to those of solving a Sylvester equation. The proposed algorithms show that in compared to literature methods, the proposed method is easy to be implemented and have high execution speed. Furthermore, we prove that the obtained Sylvester equation has a unique solution. In addition, the effect of the finite difference space step size to the computational accuracy is studied. For having suitable solution, the numerical solutions show that there is no need to select very small step size. Also only a small number of basis functions in the Legendre wavelets series is needed. The numerical results demonstrate efficiency and capability of the proposed method.