Solution of option pricing equations using orthogonal polynomial expansion
We study both analytic and numerical solutions of option pricing equations using systems of orthogonal polynomials. Using a Galerkin-based method, we solve the parabolic partial differential equation for the Black-Scholes model using Hermite polynomials and for the Heston model using Hermite and Laguerre polynomials. We compare the obtained solutions to existing semi-closed pricing formulas. Special attention is paid to the solution of the Heston model at the boundary with vanishing volatility.