Simulate Zero Coupon Bond Using Chen's Model

Define Chen's model for the evolution of interest rates:
Visualize some trajectories using the stochastic RungeKutta method, correcting for negative interest rates that result from the use of the approximation:
Evaluate zerocoupon bond maturing at time 1 using the Monte Carlo method. The present value of the bond is given by the mean of the process :