In this paper we uniformly approximate the trajectories of the Cox–Ingersoll–Ross (CIR) process. At a sequence of random times the approximate trajectories will be even exact. In between, the approximation will be uniformly close to the exact trajectory. From a conceptual point of view, the proposed method gives a better quality of approximation in a path-wise sense than standard, or even exact, simulation of the CIR dynamics at some deterministic time grid.

The Doss-Sussmann (DS) approach is used for uniform simulation of the Cox-Ingersoll-Ross (CIR) process. The DS formalism allows us to express trajectories of the CIR process through solutions of some ordinary differential equation (ODE) depending on realizations of a Wiener process involved. By simulating the first-passage times of the increments of the Wiener process to the boundary of an interval and solving the ODE, we uniformly approximate the trajectories of the CIR process. In this respect special attention is payed to simulation of trajectories near 0. From a conceptual point of view the proposed method gives a better quality of approximation (from a pathwise point of view) than standard, even exact, simulation of the stochastic differential equation at some deterministic time grid.