上式称为辛普森求积公式，也称三点公式或抛物线公式b-a+[f(a) +4f("S=I,(f)记为)+ f(b)l624.5辛普森公式的余项为3.5R(S) = R(I,)=/ R(x)dx32.5b-a2b-a)* f((n)1.521800.5o0.500.51.5上页辛普森公式具有3次代数精度下页返园
上页 下页 返回 -0.5 0 0.5 1 1.5 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 上式称为辛普森求积公式，也称三点公式或抛物线公式 记为 ( ) 2 S I f 辛普森公式的余项为 ( ) ( ) 2 R S R I b a R2 (x)dx ) ( ) 2 ( 180 4 (4) f b a b a 辛普森公式具有3次代数精度 ) ( )] 2 [ ( ) 4 ( 6 f b a b f a f b a
(-1)"-kC(n)"II(t-j)dt10n.k!(n-k)!0≤j≤n3.柯特斯公式及其余项jtkb-a取n= 4,则 x, =a+kh,k =0,1,...,4,h=47C(4)柯特斯系数为(t -1)(t -2)(t -3)(t - 4)dt =904.4!32C(4)t(t -2)(t - 3)(t - 4)dt9004.3!V121Cht(t -1)(t -3)(t -4)dt904.2!.2!321C(4)t(t -1)(t -2)(t - 4)dtC3Jo904.3!7t(t -1)(t - 2)(t - 3)dt上页904.4Jo下页返回
上页 下页 返回 3.柯特斯公式及其余项 4 4, , 0,1, ,4 , b a n xk a kh k h 取 则 柯特斯系数为 C (t 1)(t 2)(t 3)(t 4)dt 4 4! 1 4 0 (4) 0 90 7 C t(t 2)(t 3)(t 4)dt 4 3! 1 4 0 (4) 1 90 32 C t(t 1)(t 3)(t 4)dt 4 2! 2! 1 4 0 (4) 2 90 12 C t(t 1)(t 2)(t 4)dt 4 3! 1 4 0 (4) 3 90 32 C t(t 1)(t 2)(t 3)dt 4 4! 1 4 0 (4) 4 90 7 n j k j n n k n k t j t n k n k C 0 0 ( ) ( )d !( )! ( 1)
求积公式为I(f) =(b-a)Ec(" f(x)K=3212327=(b-a)[ f(x)+2.-7--f(x)+f(x2)+f(xg)+f(x))9090909090b-a[7f(x)+32f(x)+12f(x2)+32f(x,)+7f(x4))90上式称为柯特斯求积公式，也称五点公式记为 C=I(f)柯特斯公式的余项为2(b-a)")° f(6) (n)R(C)= R(I) = R (x)dx945上页柯特斯公式具有5次代数精度下页返回
上页 下页 返回 求积公式为 ( ) 4 I f 4 0 (4) ( ) ( ) k k xk b a C f ( )] 9 0 7 ( ) 9 0 3 2 ( ) 9 0 1 2 ( ) 9 0 3 2 ( ) 9 0 7 ( )[ 0 1 2 3 x4 b a f x f x f x f x f [7 ( ) 3 2 ( ) 1 2 ( ) 3 2 ( ) 7 ( )] 9 0 0 1 2 3 x4 f x f x f x f x f b a 上式称为柯特斯求积公式，也称五点公式 记为 ( ) 4 C I f 柯特斯公式的余项为 ( ) ( ) 4 R C R I b a R4 (x)dx ) ( ) 4 ( 945 2( ) 6 (6) f b a b a 柯特斯公式具有5次代数精度
21etdr的近似值，并例5试分别使用梯形公式和Simpson公式计算积分1估计截断误差[edx~2,l(e+et)=2. 1835.解用梯形公式计算：2f(z)et,f(a)et，f"(a)(2+)e,max（）/m(1)8.1548(21)#max1f(r)10.6796.估计截断误差为R121616221-e+4e+电量）2.0263用Sinpson公式计算：eid61612413612tmax4(z)1mft(1)198.43F0()(4++内利十35135(2-1)5上页Rf(*)(x) = 0.06890估计裁断误差为max<28801≤x≤2下页返园
上页 下页 返回 max ( ) 0.06890 2880 ( 2 1 ) ( 4 ) 1 2 5 2 R f x x 例 5
四、牛顿一柯特斯公式的稳定性（舍入误差考察柯特斯系数(-1)"-kC(n)II(t- j)dt1JOn.k!(n-k)!0≤jSnj+k只与积分区间a,b的节点x,的划分有关与函数f(x)无关其值可以精确给定因此用牛顿一柯特斯公式计算积分的舍入误差主要由函数值f(x)的计算引起只需讨论(x)的舍入误差对公式的影响假设f(x)为精确值而以f(x)作为f(x）的近似值计算值）8=f(x）-f(x）为误差记 i,=(b-a)Ec(")j(x)为I,的近似值计算值上页k=0而理论值为 I,=(b-a)Zcl"f(x)下页返圆k=0
上页 下页 返回 四、牛顿-柯特斯公式的稳定性(舍入误差) t j d t n k n k C n j k j n n k n k 0 0 ( ) ( ) !( )! ( 1) 考察柯特斯系数 只与积分区间[a,b]的节点xj 的划分有关,与函数f (x)无关 因此用牛顿-柯特斯公式计算积分的舍入误差主要由函数值 ( )的计算引起, xk f 其值可以精确给定 只需讨论f (xk )的舍入误差对公式的影响 假设 ( )为精确值,而以 ( )作为 ( )的近似值(计算值) k k xk f x f x f k f (xk ) f (xk )为误差 n I n k k n b a Ck f x 0 ( ) 记 ( ) ( ) 为 的近似值(计算值) n I 而理论值为 n I n k k n b a Ck f x 0 ( ) ( ) ( )