+h(bfn-+b2f, + 上式称为显式阿当姆斯( Adams公式,其余式为 R(xn, h)=h (al+1)…(u+k-1 k1,(k+1) k! h y+)(2(+1)…(a+k-1) k+1,( k! Bh+(k+(s) ∈(x n-k2n-1 其中的系数值b(i=1,2,k)和误差常数B为(课本P216) k b 2 b 3 B 23/2 1/2 5/12 323/12-16/125/12 3/8 455/24-59/2437/24-9/24251/720
上式称为显式阿当姆斯(Adams)公式，其余式为 y y h b f b f b f n n n n k n k = + + + + − − − − 1 1 1 2 2 ( ) 1 ( 1) 0 1 1 ( 1) 0 1 ( 1) 1 ( 1) ( 1) ( , ) ( ) ! ( 1) ( 1) ( ) ! ( ) ( , ) k k n k k k k n k n u u u k R x h h h y du k u u u k h y du k Bh y x x + + + + + − − + + − = + + − = = 其中的系数值bi (i=1,2,…,k)和误差常数B为（课本P216） k b1 b2 b3 b4 B 2 3/2 -1/2 5/12 3 23/12 -16/12 5/12 3/8 4 55/24 -59/24 37/24 -9/24 251/720
511利用数值积分的方法 考虑方程y=f(x,y)在区间[xn=1,xn]上的积分 yt=」.f(x,y)x y(,)=y(m-)+ f(,y)dx 如果通过等距节点xn,x212…xn-k上的已知值fn,f21…,fnk 作 Lagrange插值,代入(*)式可得隐式阿当姆斯公式 yn=yn-1+h(bo/n +b,f n-1 …… 其余式为R(xn,h)=Bh+2y+2(),5∈(x2k,xn) 其中系数值与误差常数见课本P217表74。 系数绝对值较小,因而隐式公式比显式公式的舍入误差小
5.1.1利用数值积分的方法 1 1 ( , ) n n n n x x x x y dx f x y dx − − = 考虑方程 y f x y = ( , ) 在区间 [ , ] x x n n −1 上的积分 即 1 1 ( ) ( ) ( , ) (*) n n x n n x y x y x f x y dx − = + − 如果通过等距节点 上的已知值 作Lagrange插值，代入（*）式可得隐式阿当姆斯公式 1 , , , n n n k x x x − − 1 , , , n n n k f f f − − 其余式为 ( ) * * * n n n n k n k 1 0 1 1 y y h b f b f b f = + + + + − − − * 2 ( 2) ( , ) ( ), ( , ) k k R x h B h y x x n n k n + + = − 其中系数值与误差常数见课本P217表7.4。 系数绝对值较小，因而隐式公式比显式公式的舍入误差小
512基于泰勒展式的待定系数法 对于微分方程y=f(x,y),考虑泰勒展式 V(X v(x, +h)=y(n)+ h+ 2 得 y(xn+h)≈y(xn)+ 即得欧拉公式(k=1时的阿当姆斯公式,属线性单步法) n+1 +nh=y thfo n5·n
5.1.2 基于泰勒展式的待定系数法 对于微分方程 得 即得欧拉公式（k=1时的阿当姆斯公式，属线性单步法）： y f x y = ( , ), 考虑泰勒展式 2 ( ) ( ) ( ) ( ) 1! 2! n n n n y x y x y x h y x h h + = + + + ( ) ( ) ( ) 1! n n n y x y x h y x h + + 1 ( , ) 1 ( , ) ! n n n n n n n y y h f x y f x y y h + = + = +
确定线性多步法计算公式的一般方法—待定系数法 考虑线性多步法的局部截断误差公式 R(n,h)=a y(na-a,y n+k ay(xn)-∑ay(xm)+b∑Bf(xm) j=0 ∑ay(xn)-b∑.f(x,+) 要求该计算公式是p阶的。 (x+h)=yx)+j(你气 将上式中y(xn+以及y(xn+在xn处展开 2/1(x)+y"(xn)+… 3! y(x, +jh)=y(xn)+jhy(rn) gh y x.)+
确定线性多步法计算公式的一般方法——待定系数法 1 0 0 0 0 ( , ) ( ) ( ) ( ) ( ) ( ) ( ) n k n k k n k k k k n k j n j j n j j j k k j n j j n j j j R x h y x y y x y x h f x y x h f x + + − + + + = = + + = = = − = − − + = − 考虑线性多步法的局部截断误差公式 要求该计算公式是p阶的。 将上式中y(xn+j)以及y’(xn+j)在xn处展开 2 3 ( ) ( ) ( ) ( ) ( ) ( ) ( ) 2! 3! n n n n n jh jh y x jh y x jhy x y x y x + = + + + + 2 ( ) ( ) ( ) ( ) ( ) 2 ! n n n n jh y x jh y x jhy x y x + = + + +
代入上述误差公式,合并同类项得 RQ(xn,h)=cnyxn)+chy(x)+…+Cnhy(x)+ 其中 C=00+c1+…+ak C1=a1+2a2…+kak-(Bb+B1+…+B) c,=(a1+2…+ka) (A1+2B2+…+k4B (q=2,3,…) 要使该计算公式是p阶的,只要令C=C1=…=Cp=0,Cn1≠0 据此求出a,B(=0,1,…k)即可
( ) 0 1 ( , ) ( ) ( ) ( ) p p R x h c y x c hy x c h y x n n n p n = + + + + 其中 ( ) 0 0 1 1 1 2 0 1 1 1 1 2 1 2 2 ( ) 1 1 2 ( 2 ) ! ( 1)! ( 2,3, ) k k k q q q q q k k c c k c k k q q q − − = + + + = + + − + + + = + + − + + + − = 要使该计算公式是p阶的，只要令 0 1 1 0, 0. p p c c c c = = = = + 代入上述误差公式，合并同类项得 据此求出 j j , ( 0,1, , ) j k = 即可