我正在尝试实现一个程序，即解决 ODE 的数值方法，称为Block BDF，如本文所述：https ://waset.org/journals/waset/v38/v38-49.pdf

由于它是可变步长，我需要计算 LTE（本地截断），即$LTE=y_{n+2}^{(5)}-y_{n+2}^{(4)}$. 文章中有一个公式$y_{n+2}^{(5)}$但我找不到一个$y_{n+2}^{(4)}$; 在哪里$y_{(n+2)}^{(k)}$意味着$k-$阶近似为$y_{(n+2)}$.

有人可以帮我找到四阶近似的公式吗？（由于我只是尝试实现，因此仅有关公式的信息可能会有所帮助，但鼓励进行解释）。

非常感谢！

PS。更通用的方法是 BDF：http ://en.wikipedia.org/wiki/Backward_differentiation_formula但此方法适用于 2 块（适用于$y(n+1)$和$y(n+2)$在同一步骤）。