保证 IVP 解的准确性的一种方法是使用步长 h 和 h/2 解决问题两次，并在对应于较大步长的网格点上比较答案。 但是对于较小的步长，这需要大量的计算，并且如果确定一致性不够好，则必须重复。 Fehlberg 方法是尝试解决此问题的一种方法。 它有一个程序来确定是否使用了正确的步长 h。 在每个步骤中，都会对解决方案进行两种不同的近似处理并进行比较。 如果两个答案非常一致，则接受近似值。 如果两个答案不符合指定的准确度，则减小步长。 如果答案同意比所需的有效数字更多，则增加步长。

这是我在sourceforge上找到的比较好的一个关于卫星轨道动力学的一个源码包，简单易懂，用c++写的，作者是韩冬，向他致敬！