使用程序包中的rcCIR和rcOU函数估计CIR和Vasicek随机利率模型的利率。
程序包中CIR和Vasicek的形式是：

和 ；
而论文中使用的一般是
和 。
程序包把括号外的回归速度参数代入到括号里面了，于是我在进行运行时也这么做了。

假设计算三年内的利率变化，划分为1000份，Dt就是单位变化，近似连续变化，此时3/1000记为0.003；x0为初始利率，theta为程序包中对应的值。此时回归速度=5，远期利率均值为0.03，波动率为0.005.
这时候有点不对劲，于是把远期利率均值改为0.05.
远期利率为0.05，但是最大的随机数还是这么点，按理说应该趋向0.05才对，结果最大值这么小。（实际上使用rcCIR也是这样）

但是论文中实际的结果是这样。
究竟是哪里出现了问题，是我的思路错了？
如果不能用这个程序，求推荐程序能模拟随机利率的；另外这样的利率变化图具体应该怎么画。