数值分析第五版全答案chap4资料

合集下载

数值分析第五版答案(全)

数值分析第五版答案(全)

第一章 绪论1.设0x >,x 的相对误差为δ,求ln x 的误差。

解:近似值*x 的相对误差为*****r e x xe x x δ-=== 而ln x 的误差为()1ln *ln *ln **e x x x e x =-≈进而有(ln *)x εδ≈2.设x 的相对误差为2%,求n x 的相对误差。

解:设()nf x x =,则函数的条件数为'()||()p xf x C f x = 又1'()n f x nx-=, 1||n p x nx C n n-⋅∴== 又((*))(*)r p r x n C x εε≈⋅且(*)r e x 为2((*))0.02n r x n ε∴≈3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:*1 1.1021x =,*20.031x =, *3385.6x =, *456.430x =,*57 1.0.x =⨯解:*1 1.1021x =是五位有效数字; *20.031x =是二位有效数字; *3385.6x =是四位有效数字; *456.430x =是五位有效数字; *57 1.0.x =⨯是二位有效数字。

4.利用公式(2.3)求下列各近似值的误差限:(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x .其中****1234,,,x x x x 均为第3题所给的数。

解:*41*32*13*34*151()1021()1021()1021()1021()102x x x x x εεεεε-----=⨯=⨯=⨯=⨯=⨯***124***1244333(1)()()()()1111010102221.0510x x x x x x εεεε----++=++=⨯+⨯+⨯=⨯ ***123*********123231132143(2)()()()()1111.10210.031100.031385.610 1.1021385.6102220.215x x x x x x x x x x x x εεεε---=++=⨯⨯⨯+⨯⨯⨯+⨯⨯⨯≈**24****24422*4335(3)(/)()()110.0311056.430102256.43056.43010x x x x x x xεεε---+≈⨯⨯+⨯⨯=⨯=5计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少? 解:球体体积为343V R π=则何种函数的条件数为23'4343p R V R R C V R ππ===(*)(*)3(*)r p r r V C R R εεε∴≈=又(*)1r V ε=%1故度量半径R 时允许的相对误差限为εr (V ∗)=13∗1%=13006.设028Y =,按递推公式1n n Y Y -= (n=1,2,…)计算到100Y 27.982≈(5位有效数字),试问计算100Y 将有多大误差?解:1n n Y Y -=-10099Y Y ∴=9998Y Y =9897Y Y =……10Y Y =依次代入后,有1000100Y Y =-即1000Y Y =27.982≈, 100027.982Y Y ∴=-*310001()()(27.982)102Y Y εεε-∴=+=⨯100Y ∴的误差限为31102-⨯。

数值分析课程第五版课后习题答案

数值分析课程第五版课后习题答案

*
[解] = (0.031 × 385.6) 1 × 10 − 4 + (1.1021 × 385.6) 1 × 10 −3 + (1.1021 × 0.031) 1 × 10 −3 ; 2 2 2 −3 −3 −3 = 0.59768 × 10 + 212.48488 × 10 + 0.01708255 × 10 = 213.09964255 × 10 −3 = 0.21309964255
ε * (R* ) 1 1 1 从而 ε * ( R * ) = 1% × R * ,故 ε r* ( R * ) = 。 = 1% × = * 3 300 3 R
6 、设 Y0 = 28 ,按递推公式 Yn = Yn −1 − 1 783 (n = 1,2, ) 计算到 Y100 ,若取 100
783 ≈ 27.982 (五位有效数字, )试问计算 Y100 将有多大误差? [解]令 Yn 表示 Yn 的近似值, e * (Yn ) = Yn − Yn ,则 e * (Y0 ) = 0 ,并且由 1 1 × 27.982 , Yn = Yn −1 − × 783 可知, 100 100 1 × (27.982 − 783 ) ,即 Yn − Yn = Yn −1 − Yn −1 − 100 1 2 从 e * (Yn ) = e * (Yn −1 ) − × (27.982 − 783 ) = e * (Yn − 2 ) − × (27.982 − 783 ) = , 100 100 Yn = Yn −1 − 而 e * (Y100 ) = e * (Y0 ) − (27.982 − 783 ) = 783 − 27.982 ,
而 783 − 27.982 ≤
1 1 × 10 −3 ,所以 ε * (Y100 ) = × 10 −3 。 2 2

(完整版)数值分析第五版答案(全)(最新整理)

(完整版)数值分析第五版答案(全)(最新整理)

第一章 绪论1.设,的相对误差为,求的误差。

0x >x δln x 解:近似值的相对误差为*x *****r e x x e x x δ-===而的误差为ln x ()1ln *ln *ln **e x x x e x =-≈进而有(ln *)x εδ≈2.设的相对误差为2%,求的相对误差。

x n x 解:设,则函数的条件数为()n f x x ='()||()p xf x C f x =又, 1'()n f x nx -= 1||n p x nx C n n-⋅∴==又((*))(*)r p r x n C x εε≈⋅ 且为2(*)r e x ((*))0.02n r x nε∴≈3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:,, , ,*1 1.1021x =*20.031x =*3385.6x =*456.430x =*57 1.0.x =⨯解:是五位有效数字;*1 1.1021x =是二位有效数字;*20.031x =是四位有效数字;*3385.6x =是五位有效数字;*456.430x =是二位有效数字。

*57 1.0.x =⨯4.利用公式(2.3)求下列各近似值的误差限:(1) ,(2) ,(3) .***124x x x ++***123x x x **24/x x 其中均为第3题所给的数。

****1234,,,x x x x 解:*41*32*13*34*151()1021()1021()1021()1021()102x x x x x εεεεε-----=⨯=⨯=⨯=⨯=⨯***124***1244333(1)()()()()1111010102221.0510x x x x x x εεεε----++=++=⨯+⨯+⨯=⨯***123*********123231132143(2)()()()()1111.10210.031100.031385.610 1.1021385.6102220.215x x x x x x x x x x x x εεεε---=++=⨯⨯⨯+⨯⨯⨯+⨯⨯⨯≈**24****24422*4335(3)(/)()()110.0311056.430102256.43056.43010x x x x x x x εεε---+≈⨯⨯+⨯⨯=⨯=5计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少?解:球体体积为343V R π=则何种函数的条件数为23'4343p R V R R C V R ππ===A A (*)(*)3(*)r p r r V C R R εεε∴≈=A 又%1(*)1r V ε=故度量半径R 时允许的相对误差限为εr (V ∗)=13∗1%=13006.设,按递推公式 (n=1,2,…)028Y =1n n Y Y -=-计算到(5位有效数字),试问计算将有多大误差?100Y 27.982≈100Y解: 1n n Y Y -=10099Y Y ∴=9998Y Y =9897Y Y =-……10Y Y =-依次代入后,有1000100Y Y =-即,1000Y Y =-, 27.982≈100027.982Y Y ∴=-*310001()()(27.982)102Y Y εεε-∴=+=⨯的误差限为。

数值分析课程第五版课后习题答案(李庆扬等)

数值分析课程第五版课后习题答案(李庆扬等)

数值分析课程第五版课后习题答案(李庆扬等)数值分析课程第五版课后习题答案(李庆扬等)第一章:数值分析导论1. 解答:数值分析是一门研究如何使用计算机来解决数学问题的学科。

它包括了从数学理论到计算实现的一系列技术。

数值分析的目标是通过近似的方式求解数学问题,其结果可能不是完全精确的,但是能够满足工程或科学应用的要求。

2. 解答:数值分析在实际应用中起着重要的作用。

它可以用于求解复杂的数学方程、计算机模拟及建模、数据的统计分析等等。

数值分析是科学计算和工程计算的基础,对许多领域都有着广泛的应用,如物理学、经济学、生物学等。

3. 解答:数值方法指的是使用数值计算的方式来求解数学问题。

与解析方法相比,数值方法一般更加灵活和高效,可以处理一些复杂的数学问题。

数值方法主要包括了数值逼近、插值、数值积分、数值微分、线性方程组的求解、非线性方程的求根等。

4. 解答:计算误差是指数值计算结果与精确解之间的差异。

在数值计算中,由于计算机的有限精度以及数值计算方法本身的近似性等因素,都会导致计算误差的产生。

计算误差可以分为截断误差和舍入误差两种。

第二章:数值误差分析1. 解答:绝对误差是指实际值与精确值之间的差异。

例如,对于一个计算出的数值近似解x和精确解x_0,其绝对误差为| x - x_0 |。

绝对误差可以衡量数值近似解的精确程度,通常被用作评估数值计算方法的好坏。

2. 解答:相对误差是指绝对误差与精确解之间的比值。

对于一个计算出的数值近似解x和精确解x_0,其相对误差为| (x - x_0) / x_0 |。

相对误差可以衡量数值近似解相对于精确解的精确度,常用于评估数值计算方法的收敛速度。

3. 解答:舍入误差是由于计算机的有限精度而引起的误差。

计算机中使用的浮点数系统只能表示有限的小数位数,因此在进行数值计算过程中,舍入误差不可避免地会产生。

舍入误差会导致计算结果与精确结果之间存在差异。

4. 解答:误差限度是指对于给定的数值计算问题,所能容忍的误差范围。

李庆扬-数值分析第五版第4章习题答案(20130714)

李庆扬-数值分析第五版第4章习题答案(20130714)
设 f (x) 1, x, x2 ,有
4h
A1
A0
A1
0 hA1 hA1
16 3
h3
h2 A1
h2 A1
4h A1 A0 A1
解出
8 A1 3 h
A0
4 3
h
A1
8 3
h
令 f (x) x3 ,
有 8h (h)3 4h 03 8h h3 8 h4 8 h4 0 h x3dx 0
有 h (h)4 4h 04 h h4 2 h5 0 h x4dx 1 x5 h 2 h5
3
3
33
h
5 h 5
因此,具有 3 次代数精度。
2h
2) 2h f (x)dx A1 f (h) A0 f (0) A1 f (h)
由于需要确定 3 个 未知量,因此,需要给定 3 个方程。
3
3
3
33
h
令 f (x) x4 ,
有 8h (h)4 4h 04 8h h4 16 h5 0 h x4dx 1 x5 2h 64 h5
3
3
3
3
h
5 2h 5
因此,具有 3 次代数精度。
1
3) 1 f (x)dx [ f (1) 2 f (x1) 3 f (x2)] / 3
T (k) m
4m 4m 1
T (k 1) m1
1 4m 1
T (k) m1
,
k
1, 2,3...
9、什么是高斯型求积公式?它的求积节点是如何确定的?它的代数精确度是多少?为何称 它是具有最高代数精确度的求积公式? 答: 如果求积公式
an
n
f (x)(x) d x Ak f (xk )

数值分析课程第五版课后习题答案

数值分析课程第五版课后习题答案

= (N + 1) − N 1 + N (N + 1)
=
N2
1, + N +1
∫N +1
因此
1
dx = α − β = arctan
1

N 1+ x2
N2 + N +1
9、正方形的边长大约为 100cm,应怎样测量才能使其面积误差不超过 1 cm2 ?
[ 解 ] 由 ε * ((l * )2 ) = [(l * )2 ]′ ε * (l * ) = 2l *ε * (l * ) 可 知 , 若 要 求 ε * ((l * )2 ) = 1 , 则
∆s = 2
2
2
s 所以
1 ab sin c

2
= ∆c + ∆b + ∆c ≤ ∆c + ∆b + ∆c c b tan c c b c
第二章 插值法(40-42)
1、根据(2.2)定义的范德蒙行列式,令
1
Vn
(x0
,
x1 ,,
xn−1 ,
x)
=
1
x0 xn−1
x02
x
2 n−1
2
2
2
= 0.59768 ×10−3 + 212.48488 ×10−3 + 0.01708255 ×10−3
= 213.09964255 ×10−3 = 0.21309964255
(3)
x
* 2
/
x4*

∑ e* (x2*
/
x
* 4
)
=
n k =1

∂f ∂xk

数值分析第五版答案(全)

数值分析第五版答案(全)

第一章 绪论1.设0x >,x 的相对误差为δ,求ln x 的误差。

解:近似值*x 的相对误差为*****r e x xe x x δ-=== 而ln x 的误差为()1ln *ln *ln **e x x x e x =-≈进而有(ln *)x εδ≈2.设x 的相对误差为2%,求n x 的相对误差。

解:设()nf x x =,则函数的条件数为'()||()p xf x C f x = 又1'()n f x nx-=, 1||n p x nx C n n-⋅∴== 又((*))(*)r p r x n C x εε≈⋅且(*)r e x 为2((*))0.02n r x n ε∴≈3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:*1 1.1021x =,*20.031x =, *3385.6x =, *456.430x =,*57 1.0.x =⨯解:*1 1.1021x =是五位有效数字; *20.031x =是二位有效数字; *3385.6x =是四位有效数字; *456.430x =是五位有效数字; *57 1.0.x =⨯是二位有效数字。

4.利用公式(2.3)求下列各近似值的误差限:(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x .其中****1234,,,x x x x 均为第3题所给的数。

解:*41*32*13*34*151()1021()1021()1021()1021()102x x x x x εεεεε-----=⨯=⨯=⨯=⨯=⨯***124***1244333(1)()()()()1111010102221.0510x x x x x x εεεε----++=++=⨯+⨯+⨯=⨯ ***123*********123231132143(2)()()()()1111.10210.031100.031385.610 1.1021385.6102220.215x x x x x x x x x x x x εεεε---=++=⨯⨯⨯+⨯⨯⨯+⨯⨯⨯≈**24****24422*4335(3)(/)()()110.0311056.430102256.43056.43010x x x x x x xεεε---+≈⨯⨯+⨯⨯=⨯=5计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少? 解:球体体积为343V R π=则何种函数的条件数为23'4343p R V R R C V R ππ===(*)(*)3(*)r p r r V C R R εεε∴≈=又(*)1r V ε=%1故度量半径R 时允许的相对误差限为εr (V ∗)=13∗1%=13006.设028Y =,按递推公式1n n Y Y -= (n=1,2,…)计算到100Y 27.982≈(5位有效数字),试问计算100Y 将有多大误差?解:1n n Y Y -=-10099Y Y ∴=9998Y Y =9897Y Y =……10Y Y =依次代入后,有1000100Y Y =-即1000Y Y =27.982≈, 100027.982Y Y ∴=-*310001()()(27.982)102Y Y εεε-∴=+=⨯100Y ∴的误差限为31102-⨯。

数值分析第五版答案(全)

数值分析第五版答案(全)

第一章 绪论1.设0x >,x 的相对误差为δ,求ln x 的误差。

解:近似值*x 的相对误差为*****r e x xe x x δ-=== 而ln x 的误差为()1ln *ln *ln **e x x x e x =-≈进而有(ln *)x εδ≈2.设x 的相对误差为2%,求n x 的相对误差。

解:设()nf x x =,则函数的条件数为'()||()p xf x C f x = 又1'()n f x nx-=, 1||n p x nx C n n-⋅∴== 又((*))(*)r p r x n C x εε≈⋅且(*)r e x 为2]((*))0.02n r x n ε∴≈3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:*1 1.1021x =,*20.031x =, *3385.6x =, *456.430x =,*57 1.0.x =⨯解:*1 1.1021x =是五位有效数字; *20.031x =是二位有效数字; *3385.6x =是四位有效数字; *456.430x =是五位有效数字; *57 1.0.x =⨯是二位有效数字。

4.利用公式求下列各近似值的误差限:(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x .其中****1234,,,x x x x 均为第3题所给的数。

解: {*41*32*13*34*151()1021()1021()1021()1021()102x x x x x εεεεε-----=⨯=⨯=⨯=⨯=⨯***124***1244333(1)()()()()1111010102221.0510x x x x x x εεεε----++=++=⨯+⨯+⨯=⨯ ***123*********123231132143(2)()()()()1111.10210.031100.031385.610 1.1021385.6102220.215x x x x x x x x x x x x εεεε---=++=⨯⨯⨯+⨯⨯⨯+⨯⨯⨯≈**24****24422*4335(3)(/)()()110.0311056.430102256.43056.43010x x x x x x xεεε---+≈⨯⨯+⨯⨯=⨯=5计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少 解:球体体积为343V R π=则何种函数的条件数为23'4343p R V R R C V R ππ===(*)(*)3(*)r p r r V C R R εεε∴≈=又(*)1r V ε=%1:故度量半径R 时允许的相对误差限为εε(ε∗)=13∗1%=13006.设028Y =,按递推公式1n n Y Y -=(n=1,2,…)计算到100Y 27.982≈(5位有效数字),试问计算100Y 将有多大误差解:1n n Y Y -=10099Y Y ∴=9998Y Y =9897Y Y =……10Y Y =依次代入后,有1000100Y Y =- %即1000Y Y =27.982≈, 100027.982Y Y ∴=-*310001()()(27.982)102Y Y εεε-∴=+=⨯100Y ∴的误差限为31102-⨯。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

数值分析第五版全答案c h a p4第四章 数值积分与数值微分1.确定下列求积公式中的特定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度:10121012112120(1)()()(0)();(2)()()(0)();(3)()[(1)2()3()]/3;(4)()[(0)()]/2[(0)()];hh hhh f x dx A f h A f A f h f x dx A f h A f A f h f x dx f f x f x f x dx h f f h ah f f h -----≈-++≈-++≈-++''≈++-⎰⎰⎰⎰ 解:求解求积公式的代数精度时,应根据代数精度的定义,即求积公式对于次数不超过m 的多项式均能准确地成立,但对于m+1次多项式就不准确成立,进行验证性求解。

(1)若101(1)()()(0)()hh f x dx A f h A f A f h --≈-++⎰ 令()1f x =,则1012h A A A -=++令()f x x =,则110A h A h -=-+令2()f x x =,则3221123h h A h A -=+ 从而解得011431313A hA h A h -⎧=⎪⎪⎪=⎨⎪⎪=⎪⎩令3()f x x =,则3()0hhh h f x dx x dx --==⎰⎰ 101()(0)()0A f h A f A f h --++=故101()()(0)()hh f x dx A f h A f A f h --=-++⎰成立。

令4()f x x =,则4551012()52()(0)()3hh h h f x dx x dx h A f h A f A f h h ---==-++=⎰⎰故此时,101()()(0)()hh f x dx A f h A f A f h --≠-++⎰故101()()(0)()hh f x dx A f h A f A f h --≈-++⎰具有3次代数精度。

(2)若21012()()(0)()hh f x dx A f h A f A f h --≈-++⎰令()1f x =,则1014h A A A -=++令()f x x =,则110A h A h -=-+令2()f x x =,则32211163h h A h A -=+ 从而解得011438383A h A h A h -⎧=-⎪⎪⎪=⎨⎪⎪=⎪⎩令3()f x x =,则22322()0hh h h f x dx x dx --==⎰⎰101()(0)()0A f h A f A f h --++=故21012()()(0)()hh f x dx A f h A f A f h --=-++⎰成立。

令4()f x x =,则22452264()5h h h h f x dx x dx h --==⎰⎰ 510116()(0)()3A f h A f A f h h --++= 故此时,21012()()(0)()h h f x dx A f h A f A f h --≠-++⎰因此, 21012()()(0)()hh f x dx A f h A f A f h --≈-++⎰具有3次代数精度。

(3)若1121()[(1)2()3()]/3f x dx f f x f x -≈-++⎰ 令()1f x =,则1121()2[(1)2()3()]/3f x dx f f x f x -==-++⎰令()f x x =,则120123x x =-++令2()f x x =,则22122123x x =++从而解得120.28990.5266x x =-⎧⎨=⎩或120.68990.1266x x =⎧⎨=⎩ 令3()f x x =,则11311()0f x dx x dx --==⎰⎰ 12[(1)2()3()]/30f f x f x -++≠故1121()[(1)2()3()]/3f x dx f f x f x -=-++⎰不成立。

因此,原求积公式具有2次代数精度。

(4)若20()[(0)()]/2[(0)()]hf x dx h f f h ah f f h ''≈++-⎰ 令()1f x =,则0(),hf x dx h =⎰2[(0)()]/2[(0)()]h f f h ah f f h h ''++-=令()f x x =,则200221()21[(0)()]/2[(0)()]2hhf x dx xdx h h f f h ah f f h h ==''++-=⎰⎰令2()f x x =,则23002321()31[(0)()]/2[(0)()]22hhf x dx x dx h h f f h ah f f h h ah ==''++-=-⎰⎰故有33211232112h h ah a =-=令3()f x x =,则340024441()41111[(0)()]/2[(0)()]12244hh f x dx x dx h h f f h h f f h h h h ==''++-=-=⎰⎰ 令4()f x x =,则450025551()51111[(0)()]/2[(0)()]12236hh f x dx x dx h h f f h h f f h h h h ==''++-=-=⎰⎰ 故此时,201()[(0)()]/2[(0)()],12hf x dx h f f h h f f h ''≠++-⎰ 因此,201()[(0)()]/2[(0)()]12hf x dx h f f h h f f h ''≈++-⎰ 具有3次代数精度。

2.分别用梯形公式和辛普森公式计算下列积分:12012101(1),8;4(1)(2),10;(3),4;(4),6;x x dx n x e dx n x n n ϕ-=+-===⎰⎰⎰ 解:21(1)8,0,1,,()84x n a b h f x x=====+ 复化梯形公式为781[()2()()]0.111402k k h T f a f x f b ==++=∑ 复化辛普森公式为7781012[()4()2()()]0.111576k k k k h S f a f x f x f b +===+++=∑∑ 121(1)(2)10,0,1,,()10x e n a b h f x x --===== 复化梯形公式为9101[()2()()] 1.391482k k h T f a f x f b ==++=∑ 复化辛普森公式为99101012[()4()2()()] 1.454716k k k k h S f a f x f x f b +===+++=∑∑(3)4,1,9,2,()n a b h f x =====复化梯形公式为341[()2()()]17.227742k k h T f a f x f b ==++=∑ 复化辛普森公式为3341012[()4()2()()]17.322226(4)6,0,,,()636k k k k h S f a f x f x f b n a b h f x ππ+===+++======∑∑ 复化梯形公式为 561[()2()()] 1.035622k k h T f a f x f b ==++=∑ 复化辛普森公式为5561012[()4()2()()] 1.035776k k k k h S f a f x f x f b +===+++=∑∑ 3。

直接验证柯特斯教材公式(2。

4)具有5交代数精度。

证明:柯特斯公式为01234()[7()32()12()32()7()]90ba b a f x dx f x f x f x f x f x -=++++⎰ 令()1f x =,则01234()90[7()32()12()32()7()]90ba b a f x dx b a f x f x f x f x f x b a -=-++++=-⎰令()f x x =,则2222012341()()21[7()32()12()32()7()]()902b ba a f x dx xdxb a b a f x f x f x f x f x b a ==--++++=-⎰⎰ 令2()f x x =,则23333012341()()31[7()32()12()32()7()]()903bb a a f x dx x dx b a b a f x f x f x f x f x b a ==--++++=-⎰⎰ 令3()f x x =,则34444012341()()41[7()32()12()32()7()]()904bb a a f x dx x dx b a b a f x f x f x f x f x b a ==--++++=-⎰⎰ 令4()f x x =,则45555012341()()51[7()32()12()32()7()]()905bb a a f x dx x dx b a b a f x f x f x f x f x b a ==--++++=-⎰⎰ 令5()f x x =,则56666012341()()61[7()32()12()32()7()]()906bb a a f x dx x dx b a b a f x f x f x f x f x b a ==--++++=-⎰⎰ 令6()f x x =,则012340()[7()32()12()32()7()]90hb a f x dx f x f x f x f x f x -≠++++⎰ 因此,该柯特斯公式具有5次代数精度。

4。

用辛普森公式求积分10x e dx -⎰并估计误差。

解:辛普森公式为[()4()()]62b a a b S f a f f b -+=++ 此时,0,1,(),x a b f x e -===从而有1121(14)0.632336S e e --=++= 误差为4(4)04()()()1802110.00035,(0,1)1802b a b a R f f e ηη--=-≤⨯⨯=∈5。

推导下列三种矩形求积公式:223()()()()();2()()()()();2()()()()();224b a b a ba f f x dxb a f a b a f f x dx b a f b b a a b f f x dx b a f b a ηηη'=-+-'=---''+=-+-⎰⎰⎰证明:(1)()()()(),(,)f x f a f x a a b ηη'=+-∈两边同时在[,]a b 上积分,得()()()()()bba a f x dxb a f a f x a dx η'=-+-⎰⎰ 即2()()()()()2(2)()()()(),(,)ba f f x dxb a f a b a f x f b f b x a b ηηη'=-+-'=--∈⎰ 两边同时在[,]a b 上积分,得()()()()()bba a f x dxb a f a f b x dx η'=---⎰⎰ 即22()()()()()2()(3)()()()()(),(,)22222ba f f x dxb a f b b a a b a b a b f a b f x f f x x a b ηηη'=---''++++'=+-+-∈⎰ 两连边同时在[,]a b 上积分,得2()()()()()()()22222b b b a a a a b a b a b f a b f x dx b a f f x dx x dx η''++++'=-+-+-⎰⎰⎰ 即 3()()()()();224ba ab f f x dx b a f b a η''+=-+-⎰ 6。

相关文档
最新文档