Pricing Interest Rate Derivatives in a Multifactor HJM Model with Time

We investigate the partial differential equation (PDE) for pricing interest derivatives in the multi-factor Cheyette Model, which involves time-dependent volatility functions with a special structure. The high dimensional parabolic PDE that results is solved numerically via a modified sparse grid approach, that turns out to be accurate and efficient. In addition we study the corresponding Monte Carlo simulation, which is fast since the distribution of the state variables can be calculated explicitly. The results obtained from both methodologies are compared to the known analytical solutions for bonds and caplets. When there is no analytical solution, both European and Bermudan swaptions have been evaluated using the sparse grid PDE approach that is shown to outperform the Monte Carlo simulation.