Before revision:  The high order analytical approximation is achieved by the method of extrapolation. As such, it enables us to employ fourth-order spatial discretization and an adaptive time integration with Dirichlet boundary conditions for obtaining the numerical solution of the asset option, option Greeks, and the optimal exercise boundary for each regime .
After revision:  We further compute the boundary values, asset option, option Greeks, and the optimal exercise boundary for each regime .