Fast approximations of bond option prices under CKLS models

A new computational method for approximating prices of zero-coupon bonds and bond option prices under general Chan–Karolyi–Longstaff–Schwartz models is proposed. The pricing partial differential equations are discretized using second-order finite difference approximations and an exponential time integration scheme combined with best rational approximations based on the Carathéodory–Fejér procedure is employed for solving the resulting semi-discrete equations. The algorithm has a linear computational complexity and provides accurate bond and European bond option prices. We give several numerical results which illustrate the computational efficiency of the algorithm and uniform second-order convergence rates for the computed bond and bond option prices.