通常仿真时多采用四阶龙格—库塔法，其原因就是这种计算公式的截断误差较小。（)