数值分析第四版习及答案
最新应用数值分析第四版第一章课后作业答案
第一章1、 在下列各对数中，x 是精确值 a 的近似值。
3.14,7/100)4(143.0,7/1)2(0031.0,1000/)3(1.3,)1(========x a x a x a x a ππ试估计x 的绝对误差和相对误差。
解：（1）0132.00416.01.3≈=≈-=-=a ee x a e r π （2）0011.00143.0143.07/1≈=≈-=-=a ee x a e r （3）0127.000004.00031.01000/≈=≈-=-=aee x a e r π （4）001.00143.03.147/100≈=≈-=-=aee x a e r2. 已知四个数：x 1=26.3，x 2=0.0250， x 3= 134.25，x 4=0.001。
试估计各近似数的有效位数和误差限，并估计运算μ1= x 1 x 2 x 3和μ1= x 3 x 4 /x 1的相对误差限。
解：x 1=26.3 ｎ=3 δx 1=0.05 δr x 1=δx 1/∣x 1∣=0.19011×10-2x 2=0.0250 ｎ=3 δx 2=0.00005 δr x 2=δx 2/∣x 2∣=0.2×10-2x 3= 134.25 ｎ=5 δx 3=0.005 δr x 3=δx 3/∣x 3∣=0.372×10-4x 4=0.001 ｎ=1 δx 4=0.0005 δr x 4=δx 4/∣x 4∣=0.5由公式：e r （μ）= e （μ）/∣μ∣≦1/∣μ∣Σni=1∣∂f/∂x i ∣δx ie r （μ1）≦1/∣μ1∣［x 2 x 3δx 1+ x 1 x 3δx 2 +x 1x 2δx 3］ =0.34468/88.269275 =0.0039049e r （μ2）≦1/∣μ2∣［x 3 x 4/ x 21δx 1+ x 4/ x 1δx 3 + x 3/ x 1δx 4］ =0.5019373、设精确数ａ>0，x 是ａ的近似值，x 的相对误差限是0.2，求㏑x 的相对误差限。
应用数值分析(第四版)课后习题答案第9章-推荐下载
6 2 1 4.用反幂法求矩阵 A 2 3 1
1 1 1
y (0) (1,1,1)T 。
解：y=[1,1,1]';z=y;d=0; A=[6,2,1;2,3,1;1,1,1]; for k=1:100 AA=A-6*eye(3); y=AA\z; [c,i]=max(abs(y)); if y(i)<0,c=-c;end z=y/c; if abs(c-d)<0.0001,break; end d=c end d=6+1/c
ai1 x1 aii xi ain xn xi
aij x j
i 1
i j
xj xi
n
aij
i1
i j
n
aij x j
i 1
i j
2 3 2 3.用幂法求矩阵 A 10 3 4 的强特征值和特征向量，迭代初值取 y (0) (1,1,1)T 。
最接近 6 的特征值为 6+1/c=7.2880，特征向量为 (1.0000 0.5229 0.2422)T 。 5.设 A R nn 非奇异，A 的正交分解为 A=QR，作逆序相乘 A1=RQ，试证明
（1） 若 A 对称则 A1 也对称； （2） 若 A 是上 Hessenberg 阵，则 A1 也是上 Hessenberg 阵。
最接近 6 的特征值和特征向量，迭代初值取
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电通，力1根保过据护管生高线产中0不工资仅艺料可高试以中卷解资配决料置吊试技顶卷术层要是配求指置，机不对组规电在范气进高设行中备继资进电料行保试空护卷载高问与中题带资2负料2，荷试而下卷且高总可中体保资配障料置各试时类卷，管调需路控要习试在题验最到；大位对限。设度在备内管进来路行确敷调保设整机过使组程其高1在中正资，常料要工试加况卷强下安看与全22过，22度并22工且22作尽22下可护都能1关可地于以缩管正小路常故高工障中作高资；中料对资试于料卷继试连电卷接保破管护坏口进范处行围理整，高核或中对者资定对料值某试，些卷审异弯核常扁与高度校中固对资定图料盒纸试位，卷置编工.写况保复进护杂行层设自防备动腐与处跨装理接置，地高尤线中其弯资要曲料避半试免径卷错标调误高试高等方中，案资要，料求编试技5写、卷术重电保交要气护底设设装。备备置管4高调、动线中试电作敷资高气，设料中课并技3试资件且、术卷料中拒管试试调绝路包验卷试动敷含方技作设线案术，技槽以来术、及避管系免架统不等启必多动要项方高方案中式；资，对料为整试解套卷决启突高动然中过停语程机文中。电高因气中此课资，件料电中试力管卷高壁电中薄气资、设料接备试口进卷不行保严调护等试装问工置题作调，并试合且技理进术利行，用过要管关求线运电敷行力设高保技中护术资装。料置线试做缆卷到敷技准设术确原指灵则导活：。。在对对分于于线调差盒试动处过保，程护当中装不高置同中高电资中压料资回试料路卷试交技卷叉术调时问试，题技应，术采作是用为指金调发属试电隔人机板员一进，变行需压隔要器开在组处事在理前发；掌生同握内一图部线纸故槽资障内料时，、，强设需电备要回制进路造行须厂外同家部时出电切具源断高高习中中题资资电料料源试试，卷卷线试切缆验除敷报从设告而完与采毕相用，关高要技中进术资行资料检料试查，卷和并主检且要测了保处解护理现装。场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
应用数值分析（第四版）课后习题答案第9章
应⽤数值分析（第四版）课后习题答案第9章第九章习题解答1.已知矩阵=???=4114114114,30103212321A A 试⽤格希哥林圆盘确定A 的特征值的界。
解：,24)2(,33)1(≤-≤-λλ2.设T x x x x ),...,,(321=是矩阵A 属于特征值λ的特征向量，若i x x =∞，试证明特征值的估计式∑≠=≤-n i j j ij ii aa 1λ.解：,x Ax λ=∞∞∞∞≤==x A x x Ax i λλ由 i x x =∞ 得 i n in i ii i x x a x a x a λ=++++ 11j n j i i ij i ii x ax a ∑≠==-1)(λj n j i i ij j n j i i ij i ii x a x ax a ∑∑≠=≠=≤=-11λ∑∑≠=≠=≤≤-nj i i ij i j n j i i ijii a x x a a 11λ3.⽤幂法求矩阵=1634310232A 的强特征值和特征向量，迭代初值取T y )1,1,1()0(=。
解：y=[1,1,1]';z=y;d=0;A=[2,3,2;10,3,4;3,6,1];for k=1:100y=A*z;[c,i]=max(abs(y));if y(i)<0,c=-c;endz=y/cif abs(c-d)<0.0001,break; endd=cend11.0000=c ,0.7500) 1.0000 0.5000(z 10.9999 =c ,0.7500) 1.0000 0.5000(z 11.0003 =c ,0.7500) 1.0000 0.5000(z 10.9989=c ,0.7500) 1.0000 0.5000(z 11.0040 =c ,0.7498) 1.0000 0.5000(z 10.9859=c ,0.7506) 1.0000 0.5001(z 11.04981 =c ,0.7478) 1.0000 0.4995(z 10.8316 =c ,0.7574) 1.0000 0.5020(z 11.5839 =c ,) 0.7260 1.0000 0.4928 (z 9.4706 =c ,0.8261) 1.0000 0.5280(z 17 = c ,0.5882) 1.0000 0.4118(z 11T (11)10T (10)9T (9)8T (8)7T (7)6T (6)5T (5)4T (4)3T (3)2T (2)1T (1)===========强特征值为11，特征向量为T 0.7500)1.0000 0.5000(。
华中科技大学出版社—数值分析第四版—课后习题及答案
14. 由于 x1 , x 2 , , x n 是 f ( x ) 的 n 个互异的零点，所以 f ( x) a 0 ( x x1 )( x x 2 ) ( x x n )
a 0 ( x xi ) a 0 ( x x j ) ( x xi ),
i 1 i 1 i j n n
4 7 h 3 时，取得最大值 max | l 2 ( x ) |
10 7 7 x 0 x x3 27 . k x , x , , x n 处进行 n 次拉格朗日插值，则有 6. i) 对 f ( x) x , (k 0,1, , n) 在 0 1 x k Pn ( x ) Rn ( x ) l j ( x) x k j
。
14．
1000000000 999999998 x1 1.000000, x2 1.000000 999999999 999999999 方程组的真解为 ，
x 1.00, x2 1.00 ， 而无论用方程一还是方程二代入消元均解得 1 结果十分可 靠。 s b sin ca a sin cb ab cos cc a b c tan c c s ab sin c a b c 15．
可 得
计
算
( f1 ) ln(1
( f 2 ) ln(1
x x 1
2
) )
1 ( x x 2 1) 60 104 3 103 2 x x 1 ，
2
x x 1
2
x x 1
2
1 1 104 8.33 107 60 2
。
(Y100 ) 100
应用数值分析(第四版)课后习题答案第1章
应用数值分析(第四版)课后习题答案第1章第一章习题解答1.在下列各对数中，某是精确值ａ的近似值（1）ａ=π，某=3.1（2）ａ=1/7，某=0.143（3）ａ=π/1000，某=0.0031（4）ａ=100/7，某=14.3试估计某的绝对误差和相对误差。
解：（1）e=∣3.1-π∣≈0.0416，δr=e/∣某∣≈0.0143（2）e=∣0.143-1/7∣≈0.0143δr=e/∣某∣≈0.1（3）e=∣0.0031-π/1000∣≈0.0279δr=e/∣某∣≈0.9(4)e=∣14.3-100/7∣≈0.0143δr=e/∣某∣≈0.0012.已知四个数：某1=26.3，某2=0.0250，某3=134.25，某4=0.001。
试估计各近似数的有效位数和误差限，并估计运算μ1=某1某2某3和μ1=某3某4/某1的相对误差限。
-2解：某1=26.3ｎ=3δ某1=0.05δr某1=δ某1/∣某1∣=0.19011某10-2某2=0.0250ｎ=3δ某2=0.00005δr某2=δ某2/∣某2∣=0.2某10-4某3=134.25ｎ=5δ某3=0.005δr某3=δ某3/∣某3∣=0.372某10某4=0.001ｎ=1δ某4=0.0005δr某4=δ某4/∣某4∣=0.5n由公式：er（μ）=e（μ）/∣μ∣≦1/∣μ∣Σi=1∣f/某i∣δ某ier（μ1）≦1/∣μ1∣［某2某3δ某1+某1某3δ某2+某1某2δ某3］=0.34468/88.269275=0.00390492er（μ2）≦1/∣μ2∣［-某3某4/某1δ某1+某4/某1δ某3+某3/某1δ某4］=0.497073.设精确数ａ>0，某是ａ的近似值，某的相对误差限是0.2，求㏑某的相对误差限。
n解：δr≦Σi=1∣f/某i∣δ某i=1/㏑某·1/某·δ某=δr某/㏑某=0.2/㏑某即δr≦0.2/㏑某4.长方体的长宽高分别为50cm，20cm和10cm，试求测量误差满足什么条件时其表面积的2误差不超过1cm。
《数值分析》所有参考答案
等价三角方程组
， ，
11.设计算机具有4位字长。分别用Gauss消去法和列主元Gauss消去法解下列方程组，并比较所得的结果。
解：Gauss消去法
回代
列主元Gauss消去
15.用列主元三角分解法求解方程组。其中
A= ,
解：
等价三角方程组
回代得
， ， ，
16.已知 ，求 ， ， 。
解：
, ,
17.设 。证明
，(II)
，
当 时
当 时
迭代格式(II)对任意 均收敛
3) ，
构造迭代格式 (III)
，
当 时
当 时
迭代格式(III)对任意 均收敛
4)
取格式(III)
， ， ，
4.用简单迭代格式求方程 的所有实根，精确至有3位有效数。
解：
当 时， ，
1 2
当 时
，
，
， ，
1)
迭代格式 ，
,
当 时， ，
任取 迭代格式收敛于
是中的一种向量范数。
解：
当 时存在 使得
，
,
所给 为 上的一个范数
18.设 。证明
(1) ;
(2) ;
(3) 。
解：(1)
(2)
(3)
19.设
A=
求 ， ， 及 ， 。
解： ,
Newton迭代格式
,
20.设 为 上任意两种矩阵(算子)范数，证明存在常数
， 使得
对一切 均成立。
解：由向量范数的等价性知道存在正常数 使得
,
=0.187622
[23.015625 , 23.015625+0.187622]
数值分析第四版课后答案答案第八章
第八章 常微分方程初值问题数值解法1、解：欧拉法公式为221(,)(100),0,1,2+=+=++=n n n n n n n y y hf x y y h x y n代00y =入上式，计算结果为 123(0.1)0.0,(0.2)0.0010,(0.3)0.00501≈=≈=≈=y y y y y y2、解：改进的欧拉法为1112[(,)(,(,))]n n n n n n n n y y h f x y f x y hf x y ++=+++将2(,)=+-f x y x x y 代入上式，得2111111221n n n n n n h hh x x x x y h y +++)+[(-)(+)+(+)]=(-+ 同理，梯形法公式为211122[(1)(1)]-+++++=++++h h n nn n n n h h y y x x x x 将00,0.1y h ==代入上二式，，计算结果见表9—5表 9—5可见梯形方法比改进的欧拉法精确。
3、证明：梯形公式为111[(,)(,)]2n n n n n n hy y f x y f x y +++=++代(,)f x y y =-入上式，得11[]2++=+--n n n n hy y y y解得21110222()()()222n n n n h h h y y y y h h h++----===⋯=+++ 因为01y =，故2()2nn h y h-=+ 对0x∀>，以h 为步长经n 步运算可求得()y x 的近似值n y ，故,,xx nh n h==代入上式有2()2x hn hy h-=+22220000222lim lim()lim(1)lim[(1)]222x x h h xx h h h h hn h h h h h h h y e h h h+-+→→→→-==-=-=+++4、解：令2()xt y x e dt =⎰，则有初值问题2',(0)0x y e y ==对上述问题应用欧拉法，取h=0.5,计算公式为210.5,0,1,2,3n x n n y y e n +=+=由0(0)0,y y ==得1234(0.5)0.5,(1.0) 1.142012708(1.5) 2.501153623,(2.0)7.245021541≈=≈=≈=≈=y y y y y y y y5、解： 四阶经典龙格-库塔方法计算公式见式（9.7）。
数值分析第四版课后习题答案
第一章习题解答1、 在下列各对数中，x 是精确值 a 的近似值。
3.14,7/100)4(143.0,7/1)2(0031.0,1000/)3(1.3,)1(========x a x a x a x a ππ试估计x 的绝对误差和相对误差。
解：（1）0132.00416.01.3≈=≈−=−=aee x a e r π （2）0011.00143.0143.07/1≈=≈−=−=a ee x a e r （3）0127.000004.00031.01000/≈=≈−=−=aee x a e r π （4）001.00143.03.147/100≈=≈−=−=aee x a e r2、已知四个数：001.0,25.134,0250.0,3.264321====x x x x 。
试估计各近似数的有效位数和误差限，并估计运算3211x x x =μ和1431/x x x =μ的相对误差限。
解：21111121101901.0,1021,3,10263.06.23−−⨯≈=⨯==⨯==x x x x n x r δδδ22214212102.0,1021,3,10250.00250.0−−−⨯≈=⨯==⨯==x x x x n x r δδδ 43332333103724.0,1021,5,1013425.025.134−−⨯≈=⨯==⨯==x x x x n x r δδδ 5.0,1021,1,101.0001.04443424==⨯==⨯==−−x x x x n x r δδδ 由相对误差限公式：i r ini n in ni i ir x x fx x f x x x f x x f u δδδ∂∂=∂∂=∑∑==1111),,(),,()(所以有：232123113211103938.0)(1)(−⨯≈++=x x x x x x x x x r δδδμμδ4971.0)(1)(4133141214311≈++−=x x x x x x x x x x r δδδμμδ 3、设精确数ａ>0，x 是ａ的近似值，x 的相对误差限是0.2，求㏑x 的相对误差限。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Yn
Yn1
1 100
783
( n=1,2,…)
计算到 Y100 .若取 783 ≈27.982(五位有效数字),试问计算Y100 将有多大误差?
7. 求方程 x2 56x 1 0 的两个根,使它至少具有四位有效数字( 783 ≈27.982).
8.
当 N 充分大时,怎样求
N
1
1 x2
dx
24.
将
f
(x)
sin
1 2
x 在 1,1 上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方逼
近多项式并画出误差图形,再计算均方误差.
25. 把 f (x) arccos x 在 1,1 上展成切比雪夫级数.
26. 用最小二乘法求一个形如 y a bx2 的经验公式,使它与下列数据拟合,并求均方误差.
第四版 数值分析习题
第一章 绪 论
1. 设 x>0,x 的相对误差为δ ,求ln x 的误差.
2. 设 x 的相对误差为 2％,求 xn 的相对误差.
3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出 它们是几位有效数字:
x1* 1.1021, x2* 0.031, x3* 385.6, x4* 56.430, x5* 71.0.
19
25
31
38
44
xi
19.0
32.3
49.0
73.3
97.8
yi
27. 观测物体的直线运动,得出以下数据:
x2 C 0,1 的最佳平方逼近,并比较其结果.
22. f (x) x 在 1,1 上,求在 1 span 1, x2, x4 上的最佳平方逼近.
sin(n 1) arccos x
un (x)
23.
1 x2
是第二类切比雪夫多项式,证明它有递推关系
un1 x 2xun x un1 x .
5.
选取常数
a
,使
max
0 x 1
x3
ax
达到极小,又问这个解是否唯一?
6. 求 f (x) sin x 在 0, / 2 上的最佳一次逼近多项式,并估计误差.
7. 求 f (x) ex 在 0,1上的最佳一次逼近多项式.
8. 如何选取 r ,使 p(x) x2 r 在1,1 上与零偏差最小? r 是否唯一?
0.53
yj
0.5000
0.5477
0.6245
0.6708
0.7280
试求三次样条插值 S(x) 并满足条件
i) S(0.25) 1.0000, S(0.53) 0.6868; ii) S(0.25) S(0.53) 0.
25. 若 f (x) C2 a,b, S(x) 是三次样条函数,证明
马克劳林级数部分和误差做比较. 2. 求证:
(a)当 m f (x) M 时, m Bn ( f , x) M . (b)当 f (x) x 时, Bn ( f , x) x .
3. 在次数不超过 6 的多项式中,求 f (x) sin 4x 在0, 2 的最佳一致逼近多项式.
4. 假设 f (x) 在 a,b 上连续,求 f (x) 的零次最佳一致逼近多项式.
0.4
0.5
0.6
0.7
lnx
-0.916291
-0.693147
-0.510826
-0.357765
0.8 -0.223144
4. 给出 cos x,0°≤x ≤90°的函数表,步长 h =1′=(1/60)°,若函数表具有 5 位有效数字, 研究用线性插值求 cos x 近似值时的总误差界.
5.
16. f (x) x7 x4 3x 1,求 f 20, 21,L , 27 及 f 20, 21,L , 28 .
17. 证明两点三次埃尔米特插值余项是
R3(x) f (4) ()(x xk )2 (x xk1)2 / 4!, (xk , xk1)
并由此求出分段三次埃尔米特插值的误差限.
b f (x)2dx bS(x)2dx b f (x) S(x)2dx 2 b S(x) f (x) S(x)dx
i) a
a
a
a
;
ii) 若 f (xi ) S(xi )(i 0,1,L , n) , 式 中 xi 为 插 值 节 点 , 且 a x0 x1 L xn b , 则
a, b, c.证明面积的误差 s 满足
s a b c . s abc
第二章 插值法
1. 根据(2.2)定义的范德蒙行列式,令
1 x0 x02 L x0n
LL LLL
Vn (x) Vn (x0 , x1,L , xn1, x) 1
xn1
x2 n1
L
xn n1
1 x x2 L xn
证明Vn (x) 是 n 次多项式,它的根是 x0,L , xn1 ,且
Vn (x) Vn1(x0, x1,L , xn1)(x x0 )L (x xn1) .
2. 当 x= 1 , -1 , 2 时, f(x)= 0 , -3 , 4 ,求 f(x)的二次插值多项式. 3. 给出 f(x)=ln x 的数值表用线性插值及二次插值计算 ln 0.54 的近似值.
x
Fn*(x) Hn 也是奇(偶)函数.
17.
求 a 、b 使
2 0
ax
b
sin
x2dx
为最小.并与
1
题及
6
题的一次逼近多项式误差作比较.
18. f (x) 、 g(x) C1 a,b,定义
b
b
(a)( f , g) a f (x)g(x)dx;(b)( f , g) a f (x)g(x)dx f (a)g(a);
设 f (x) C2 a,b 且 f (a) f (b) 0 ,求证 axb
f (x) 1 (b a)2
8
a xb
f (x) .
8. 在 4 x 4 上给出 f (x) ex 的等距节点函数表,若用二次插值求 ex 的近似值,要使截
断误差不超过106 ,问使用函数表的步长 h 应取多少?
改用另一等价公式
ln(x x2 1) ln(x x2 1)
计算,求对数时误差有多大?
x11010 x2 1010 ;
14. 试用消元法解方程组 x1 x2 2.
假定只用三位数计算,问结果是否可靠?
s 1 ab sin c,
0c
15. 已知三角形面积 2
其中 c 为弧度,
2 ,且测量 a ,b ,c 的误差分别为
18. 求一个次数不高于 4 次的多项式 P(x) ,使它满足 P(0) P(k 1) 并由此求出分段三次
埃尔米特插值的误差限.
19. 试求出一个最高次数不高于 4 次的函数多项式 P(x) ,以便使它能够满足以下边界条件
P(0) P(0) 0 , P(1) P(1) 1, P(2) 1.
问它们是否构成内积?
1 x6 dx
19. 用许瓦兹不等式(4.5)估计 0 1 x 的上界,并用积分中值定理估计同一积分的上下界,
并比较其结果.
20.
选择 a ,使下列积分取得最小值:
1 (x ax2 )2dx,
1
1 x ax2 dx
1
.
21. 设空间 span1, x, 2 span x100, x101 ,分别在 1 、2 上求出一个元素,使得其为
11. 证明 ( fk gk ) fk gk gk1fk .
n1
n1
fk gk fn gn f0g0 gk1fk .
12. 证明 k0
k 0
n1
2 y j yn y0.
13. 证明 j0
14. 若 f (x) a0 a1x L an1xn1 anxn 有 n 个不同实根 x1, x2,L , xn ,证明
n
xkj
j1 f (x j )
0,0k n2; an1 ,k n1.
15. 证明 n 阶均差有下列性质:
i) 若 F(x) cf (x) ,则 F x0, x1,L , xn cf x0, x1,L , xn ;
ii) 若 F(x) f (x) g(x) ,则 F x0, x1,L , xn f x0, x1,L , xn g x0, x1,L , xn .
20. 设 f (x) C a,b,把 a,b 分为 n 等分,试构造一个台阶形的零次分段插值函数 n (x)
并证明当 n 时, n (x) 在a,b 上一致收敛到 f (x) .
21. 设 f (x) 1/(1 x2 ) ,在 5 x 5 上取 n 10 ,按等距节点求分段线性插值函数 Ih (x) ,
?
9. 正方形的边长大约为 100 ㎝,应怎样测量才能使其面积误差不超过 1 ㎝ 2 ?
S 1 gt2 10. 设 2 假定 g 是准确的,而对 t 的测量有±0.1 秒的误差,证明当 t 增加时 S 的绝对
误差增加,而相对误差却减小.
11. 序列{yn} 满足递推关系 yn 10yn1 1 (n=1,2,…),若 y0 2 1.41 (三位有效数字),
设 xk
x0
kh
,k=0,1,2,3,求
max
x0 xx3
l2 (x)
.
6. 设 x j 为互异节点(j=0,1,…,n),求证:
n
xkj l j (x) xk (k 0,1,L , n);
i) j0
n
(xj x)k l j (x) k 1, 2,L , n).
ii) j0
max max 7.
9. 设 f (x) x4 3x3 1,在 0,1 上求三次最佳逼近多项式.