suanfa

代码说明：

   数值解与理论解对比可知，四阶龙格-库塔法的精度已经很高，用它来解一般常微分方程已经足够了。(Numerical comparison shows that the theoretical solutions, Runge- Kutta method has high accuracy, and use it to solve ordinary differential equations general enough.)