We consider a system of coupled free boundary problems for pricing American put options with regime switching. To solve this system, we first fix the optimal exercise boundary for each regime resulting in multi-variable fixed domains. We further eliminate the first order derivatives associated with the regime switching model by taking derivatives to obtain a system of coupled partial differential equations which we called the asset-delta-gamma-speed option equations. The fourth-order compact finite difference scheme and Gauss-Seidel iterative method are then employed in each regime for solving the system of the equations. In particular, the third order Hermite interpolation technique is used for estimating the coupled asset and delta options in the set of equations. The numerical method is finally tested with several examples. Our results show that the scheme provides an accurate solution with the convergent rate in space of 2.44 and the rate in time of 1.86, which is accurate and fast in computation as compared with other existing numerical methods.
↧