Simulate Zero Coupon Bond Using Chen's Model

Define Chen's model for the evolution of interest rates:

Visualize few trajectories using the stochastic Runge-Kutta method, making sure to correct the artifact of getting negative interest rates resulting from the use of an approximate simulation method:

Evaluate zero-coupon bond maturing at time 1 using Monte Carlo. It is given by the mean of process .