Comparison of numerical methods to price zero coupon bonds in a two-factor CIR model

In this paper we price a zero coupon bond under a Cox–Ingersoll–Ross (CIR) two-factor model using various numerical schemes. To the best of our knowledge, a closed-form or explicit price functional is not trivial and has been less studied. The use and comparison of several numerical methods to determine the bond price is one contribution of this paper. Ordinary differential equations (ODEs) , finite difference schemes and simulation are the three classes of numerical methods considered. These are compared on the basis of computational efficiency and accuracy, with the second aim of this paper being to identify the most efficient numerical method. The numerical ODE methods used to solve the system of ODEs arising as a result of the affine structure of the CIR model are more accurate and efficient than the other classes of methods considered, with the Runge–Kutta ODE method being the most efficient. The Alternating Direction Implicit (ADI) method is the most efficient of the finite difference scheme methods considered, while the simulation methods are shown to be inefficient. Our choice of considering these methods instead of the other known and apparently new numerical methods (eg Fast Fourier Transform (FFT) method, Cosine (COS) method, etc.) is motivated by their popularity in handling interest rate instruments.