常微分方程组和高阶微分方程的数值解法.PDF

R-K方法的稳定区域 将R-K方法用到实验方程y =y 上得  h 2   K1      K1 hyn , K 2 h yn yn h   2   2   h 2 h 3   K2      K 3 h yn  yn h    2   2 4   h 3 h 4  2         K h y K y h h  4  n 3  n    2 4    1 y n1 y n  K 1 2K 2 2K 3 K 4  6  h 2 h 3 h 4        1    y n  h  2 6 24    因此误差方程是：  h 2 h 3 h 4          1 h n1 n    2  6  24    R K 于是  公式的绝对稳定区域为 2 3 4 h h h       1 h