想要使用四阶龙格库塔法生成模拟数据，但是R里面没有数值解的包，找了一番，发现deSolve里面就有rk4函数，但是现在问题来了，我不知道rk4函数如何应用于常微分方程组...
1.我要做的是以下论文的3.1部分也就是模拟数据的生成，如果不用Matlab，能否用deSolve包的rk4函数来生成呢？
2.局部线性估计和局部二次估计又该用哪个包呢

谢谢各位大神指教