对于常微分方程组，很难求得解析解，因此工程上常用数值解法，常微分方程组的R-K解法-For ordinary differential equations, analytical solution is difficult to obtain, so the numerical method used on the project,

RK solution of ordinary differential equations