吉林大学工程数学计算方法(第三章习题答案)
计算方法第三章习题答案
计算方法第三章习题答案计算方法第三章习题答案计算方法是一门涵盖了数值计算和计算机编程的学科，它在现代科学和工程中扮演着重要的角色。
第三章是计算方法课程中的重要章节，主要涉及到数值计算中的误差分析和插值方法。
本文将为大家提供第三章习题的详细答案，帮助读者更好地理解和应用这些概念。
1. 误差分析误差分析是计算方法中非常重要的一部分，它帮助我们理解和评估数值计算中的误差来源。
以下是一些常见的误差类型：- 绝对误差：绝对误差是指数值计算结果与真实值之间的差异。
它可以通过计算两者之差来得到。
- 相对误差：相对误差是指绝对误差与真实值之间的比值。
通常以百分比的形式表示。
- 截断误差：截断误差是由于在计算过程中舍入或截断数字而引入的误差。
它通常是由于计算机的有限精度导致的。
- 舍入误差：舍入误差是由于将无限位数的小数截断为有限位数而引入的误差。
它通常是由于计算机的有限精度或计算方法的近似性质导致的。
2. 插值方法插值方法是一种用于通过已知数据点来估计未知数据点的技术。
以下是一些常见的插值方法：- 线性插值：线性插值是一种简单的插值方法，它假设两个已知数据点之间的未知数据点的取值在直线上。
通过已知数据点的斜率和截距，我们可以计算出未知数据点的值。
- 拉格朗日插值：拉格朗日插值是一种使用多项式来逼近已知数据点的方法。
它通过构造一个满足通过已知数据点的多项式来估计未知数据点的值。
- 牛顿插值：牛顿插值是一种使用差商来逼近已知数据点的方法。
它通过构造一个满足通过已知数据点的差商多项式来估计未知数据点的值。
3. 习题答案以下是一些第三章习题的答案，供大家参考：- 习题1：已知函数f(x)在区间[a, b]上连续，且在[a, b]上的导数存在且连续，证明存在一点c∈(a, b)，使得f(b) - f(a) = (b - a)f'(c)。
这是拉格朗日中值定理的一个特例，根据定理的条件，我们可以得到上述结论。
- 习题2：已知函数f(x)在区间[a, b]上连续，且在(a, b)内可导，证明存在一点c∈(a, b)，使得f'(c) = (f(b) - f(a))/(b - a)。
计算方法的课后答案
《计算方法》习题答案第一章 数值计算中的误差1．什么是计算方法？（狭义解释）答：计算方法就是将所求的的数学问题简化为一系列的算术运算和逻辑运算，以便在计算机上编程上机，求出问题的数值解，并对算法的收敛性、稳定性和误差进行分析、计算。
2．一个实际问题利用计算机解决所采取的五个步骤是什么？答：一个实际问题当利用计算机来解决时，应采取以下五个步骤： 实际问题→建立数学模型→构造数值算法→编程上机→获得近似结果 4．利用秦九韶算法计算多项式4)(53-+-=x x x x P 在3-=x 处的值，并编程获得解。
解：400)(2345-+⋅+-⋅+=x x x x x x P ，从而所以，多项式4)(53-+-=x x x x P 在3-=x 处的值223)3(-=-P 。
5．叙述误差的种类及来源。
答：误差的种类及来源有如下四个方面：（1）模型误差：数学模型是对实际问题进行抽象，忽略一些次要因素简化得到的，它是原始问题的近似，即使数学模型能求出准确解，也与实际问题的真解不同，我们把数学模型与实际问题之间存在的误差称为模型误差。
（2）观测误差：在建模和具体运算过程中所用的一些原始数据往往都是通过观测、实验得来的，由于仪器的精密性，实验手段的局限性，周围环境的变化以及人们的工作态度和能力等因素，而使数据必然带有误差，这种误差称为观测误差。
（3）截断误差：理论上的精确值往往要求用无限次的运算才能得到，而实际运算时只能用有限次运算的结果来近似，这样引起的误差称为截断误差（或方法误差）。
（4）舍入误差：在数值计算过程中还会用到一些无穷小数，而计算机受机器字长的限制，它所能表示的数据只能是一定的有限数位，需要把数据按四舍五入成一定位数的近似的有理数来代替。
这样引起的误差称为舍入误差。
6．掌握绝对误差（限）和相对误差（限）的定义公式。
答：设*x 是某个量的精确值，x 是其近似值，则称差x x e -=*为近似值x 的绝对误差（简称误差）。
计算方法各习题及参考答案
第二章 数值分析2.1 已知多项式432()1p x x x x x =-+-+通过下列点：试构造一多项式()q x 通过下列点：答案：54313()()()3122q x p x r x x x x x =-=-++-+． 2.2 观测得到二次多项式2()p x 的值：表中2()p x 的某一个函数值有错误，试找出并校正它．答案：函数值表中2(1)p -错误，应有2(1)0p -=．2.3 利用差分的性质证明22212(1)(21)/6n n n n +++=++ ．2.4 当用等距节点的分段二次插值多项式在区间[1,1]-近似函数xe 时，使用多少个节点能够保证误差不超过61102-⨯． 答案：需要143个插值节点．2.5 设被插值函数4()[,]f x C a b ∈，()3()h H x 是()f x 关于等距节点01n a x x x b =<<<= 的分段三次艾尔米特插值多项式，步长b ah n-=．试估计()3||()()||h f x H x ∞-．答案：()443||()()||384h M f x H x h ∞-≤．第三章 函数逼近3.1 求()sin ,[0,0.1]f x x x =∈在空间2{1,,}span x x Φ=上最佳平方逼近多项式，并给出平方误差．答案：()sin f x x =的二次最佳平方逼近多项式为-522sin ()0.832 440 710 1.000 999 10.024 985 1x p x x x ≈=-⨯+-，二次最佳平方逼近的平方误差为0.122-1220(sin )())0.989 310 710x p x dx δ=-=⨯⎰．3.2 确定参数,a b c 和，使得积分2121(,,)[I a b c ax bx c -=++-⎰取最小值．答案：810, 0, 33a b c ππ=-== 3.3 求多项式432()251f x x x x =+++在[1,1]-上的３次最佳一致逼近多项式()p x ．答案：()f x 的最佳一致逼近多项式为323()74p x x x =++． 3.4 用幂级数缩合方法，求() (11)x f x e x =-≤≤上的３次近似多项式6,3()p x ，并估计6,3||()()||f x p x ∞-．答案：236,3()0.994 574 650.997 395 830.542 968 750.177 083 33p x x x x =+++， 6,3||()()||0.006 572 327 7f x p x ∞-≤3.5 求() (11)xf x e x =-≤≤上的关于权函数()x ρ=的三次最佳平方逼近多项式3()S x ，并估计误差32||()()||f x S x -和3||()()||f x S x ∞-．答案：233()0.994 5710.997 3080.542 9910.177 347S x x x x =+++，32||()()||0.006 894 83f x S x -=，3||()()||0.006 442 575f x S x ∞-≤．第四章 数值积分与数值微分4.1 用梯形公式、辛浦生公式和柯特斯公式分别计算积分1(1,2,3,4)n x dx n =⎰，并与精确值比较．答案：计算结果如下表所示4.2 确定下列求积公式中的待定参数，使得求积公式的代数精度尽量高，并指明所确定的求积公式具有的代数精度． （１）101()()(0)()hh f x dx A f h A f A f h --≈-++⎰（２）11211()[(1)2()3()]3f x dx f f x f x -≈-++⎰ （３）20()[(0)()][(0)()]2h h f x dx f f h h f f h α''≈++-⎰答案：（１）具有三次代数精确度（２）具有二次代数精确度（３）具有三次代数精确度．4.3 设10h x x =-，确定求积公式12300101()()[()()][()()][]x x x x f x dx h Af x Bf x h Cf x Df x R f ''-=++++⎰中的待定参数,,,A B C D ，使得该求积公式的代数精确度尽量高，并给出余项表达式．答案：3711,,,20203020A B C D ====-，(4)6()[]1440f R f h η=，其中01(,)x x η∈．4.4 设2()P x 是以0,,2h h 为插值点的()f x 的二次插值多项式，用2()P x 导出计算积分30()hI f x dx =⎰的数值积分公式h I ，并用台劳展开法证明：453(0)()8h I I h f O h '''-=+． 答案：3203()[(0)3(2)]4h h I p x dx h f f h ==+⎰．4.5 给定积分10sin xI dx x =⎰（１）运用复化梯形公式计算上述积分值，使其截断误差不超过31102-⨯． （２）取同样的求积节点，改用复化辛浦生公式计算时，截断误差是多少？（３）要求的截断误差不超过610-，若用复化辛浦生公式，应取多少个节点处的函数值？ 答案：（１）只需7.5n ≥，取９个节点，0.946I ≈（２）4(4)46111|[]||()|()0.271102880288045n b a R f h f η--=-≤=⨯ （３）取７个节点处的函数值．4.6 用变步长的复化梯形公式和变步长的复化辛浦生公式计算积分10sin xI dx x =⎰．要求用事后误差估计法时，截断误不超过31102-⨯和61102-⨯． 答案：使用复化梯形公式时，80.946I T ≈=满足精度要求；使用复化辛浦生公式时，40.946 083I s ≈=满足精度要求．4.7（１）利用埃尔米特插值公式推导带有导数值的求积公式2()()[()()][()()][]212ba b a b a f x dx f a f b f b f a R f --''=+--+⎰，其中余项为 5(4)()[](), (,)4!30b a R f f a b ηη-=∈． （２）利用上述公式推导带修正项的复化梯形求积公式020()[()()]12Nx N N x h f x dx T f x f x ''≈--⎰，其中 0121[()2()2()2()()]2N N N hT f x f x f x f x f x -=+++++ ，而 00, (0,1,2,,), i N x x ih i N Nh x x =+==- ．4.8 用龙贝格方法计算椭圆2214x y +=的周长，使结果具有五位有效数字． 答案：49.6884l I =≈．4.9确定高斯型求积公式0011()()()x dx A f x A f x ≈+⎰的节点0x ，1x 及系数0A ，1A ．答案：00.289 949x =，10.821 162x =，00.277 556A =，10.389 111A =．4.10 验证高斯型求积公式00110()()()x e f x dx A f x A f x +∞-≈+⎰的系数及节点分别为0001 2 2A A x x ===-=+第五章 解线性方程组的直接法5.1 用按列选主元的高斯-若当消去法求矩阵A 的逆矩阵，其中111210110A -⎛⎫ ⎪= ⎪ ⎪-⎝⎭． 答案： 1110331203321133A -⎛⎫ ⎪ ⎪⎪=- ⎪ ⎪⎪-- ⎪⎝⎭5.2 用矩阵的直接三角分解法解方程组1234102050101312431701037x x x x ⎛⎫⎛⎫⎛⎫⎪ ⎪ ⎪⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭ 答案： 42x =，32x =，21x =，11x =．5.3 用平方根法(Cholesky 分解法)求解方程组12341161 4.25 2.750.51 2.75 3.5 1.25x x x -⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪-=- ⎪⎪ ⎪ ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭答案： 12x =，21x =，31x =-．5.4 用追赶法求解三对角方程组123421113121112210x x x x ⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭ 答案：42x =，31x =-，21x =，10x =．第六章 解线性代数方程组的迭代法６.１ 对方程1212123879897x x x x x x x -+=⎧⎪-+=⎨⎪--=⎩作简单调整，使得用高斯－赛得尔迭代法求解时对任意初始向量都收敛，并取初始向量(0)[0 0 0]T x =，用该方法求近似解(1)k x+，使(1)()3||||10k k x x +-∞-≤． 答案：近似解为(4)[1.0000 1.0000 1.0000]Tx =．６.２ 讨论松弛因子 1.25ω=时，用SOR 方法求解方程组121232343163420412x x x x x x x +=⎧⎪+-=⎨⎪-+=-⎩ 的收敛性．若收敛，则取(0)[0 0 0]T x=迭代求解，使(1)()41||||102k k x x +-∞-<⨯． 答案：方程组的近似解为*1 1.50001x =，*2 3.33333x =，*3 2.16667x =-．６.３ 给定线性方程组Ax b =，其中111221112211122A ⎛⎫ ⎪ ⎪⎪= ⎪ ⎪ ⎪ ⎪⎝⎭，证明用雅可比迭代法解此方程组发散，而高斯－赛得尔迭代法收敛．６.４ 设有方程组112233302021212x b x b x b -⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪= ⎪⎪ ⎪ ⎪⎪ ⎪-⎝⎭⎝⎭⎝⎭，讨论用雅可比方法和高斯－赛得尔方法解此方程组的收敛性．如果收敛，比较哪种方法收敛较快．答案：雅可比方法收敛，高斯－赛得尔方法收敛,且较快．６.５ 设矩阵A 非奇异．求证：方程组Ax b =的解总能通过高斯－赛得尔方法得到．６.６ 设()ij n nA a ⨯=为对称正定矩阵，对角阵1122(,,,)nn D diag a a a = ．求证：高斯－赛得尔方法求解方程组1122D AD x b --=时对任意初始向量都收敛．第七章 非线性方程求根例７.４ 对方程230xx e -=确定迭代函数()x ϕ及区间[,]a b ，使对0[,]x a b ∀∈，迭代过程1(), 0,1,2,k x x k ϕ+== 均收敛，并求解．要求51||10k k x x -+-<． 答案：若取2()x x ϕ=，则在[1,0]-中满足收敛性条件，因此迭代法121, 0,1,2,k x k x k +== 在(1,0)-中有惟一解．取00.5x =-，*70.458960903x x ≈=-．取2()x x ϕ=，在[0,1上满足收敛性条件，迭代序列121, 0,1,2,k x k x k +== 在[0,1]中有惟一解．取00.5x =，*140.910001967x x ≈=- 在[3,4]上，将原方程改写为23xe x =，取对数得2ln(3)()x x x ϕ==．满足收敛性条件，则迭代序列21ln(3), 0,1,2,k k x x k +== 在[3,4]中有惟一解．取0 3.5x =， *16 3.733067511x x ≈=．例７.６ 对于迭代函数2()(3)x x c x ϕ=+-，试讨论：（１）当c 为何值时，1()k k x x ϕ+=产生的序列{}k x（２）c 取何值时收敛最快？（３）取1,2c =-()x ϕ51||10k k x x -+-<．答案：（１）(c ∈时迭代收敛．（２）c =时收敛最快．（３）分别取1, 2c =--，并取0 1.5x =，计算结果如下表７.７所示表７.７例７.13 设不动点迭代1()k x x ϕ+=的迭代函数()x ϕ具有二阶连续导数，*x 是()x ϕ的不动点，且*()1x ϕ'≠，证明Steffensen 迭代式21(), (), 0,1,2,()2k k k k k k k k k k k y x z x k y x x x z y xϕϕ+==⎧⎪=-⎨=-⎪-+⎩二阶收敛于*x ．例7.15 设2()()()()()x x p x f x q x f x ϕ=--，试确定函数()p x 和()q x ，使求解()0f x =且以()x ϕ为迭代函数的迭代法至少三阶收敛．答案：1()()p x f x ='，31()()2[()]f x q x f x ''=' 例７.19 设()f x 在[,]a b 上有高阶导数，*(,)x a b ∈是()0f x =的(2)m m ≥重根，且牛顿法收敛，证明牛顿迭代序列{}k x 有下列极限关系：111lim2k kk k k k x x m x x x -→∞-+-=-+．第八章 矩阵特征值8.1 用乘幂法求矩阵A 的按模最大的特征值与对应的特征向量，已知5500 5.51031A -⎛⎫ ⎪=- ⎪ ⎪-⎝⎭，要求(1)()611||10k k λλ+--<，这里()1k λ表示1λ的第k 次近似值．答案：15λ≈，对应的特征向量为[5,0,0]T-；25λ≈-，对应的特征向量为[5,10,5]T --． 8.2 用反幂法求矩阵110242012A -⎛⎫⎪=-- ⎪ ⎪-⎝⎭的按模最小的特征值．知A 的按模较大的特征值的近似值为15λ=，用5p =的原点平移法计算1λ及其对应的特征向量．答案：(1) A 的按模最小的特征值为30.2384428λ≈(2) 1 5.1248854λ≈，对应的特征向量为(8)[0.242 4310, 1 ,0.320 011 7]T U =--．8.3 设方阵A 的特征值都是实数，且满足121, ||||n n λλλλλ>≥≥> ，为求1λ而作原点平移，试证：当平移量21()2n p λλ=+时，幂法收敛最快． 8.4 用二分法求三对角对称方阵1221221221A ⎛⎫ ⎪ ⎪= ⎪ ⎪ ⎪⎝⎭的最小特征值，使它至少具有２位有效数字．答案：取5 2.234375λ≈-即有２位有效数字．8.5 用平面旋转变换和反射变换将向量[2 3 0 5]T x =变为与1[1 0 0 0]Te =平行的向量．答案：203/2/00001010/0T ⎛⎫⎪- ⎪=⎪--⎝0.324 442 8400.486 664 26200.811 107 1040.486 664 2620.812 176 04800.298 039 92200100.811 107 1040.298 039 92200.530 266 798H --⎛⎫⎪--⎪= ⎪ ⎪⎪--⎝⎭8.6 若532644445A -⎛⎫⎪=- ⎪ ⎪-⎝⎭，试把A 化为相似的上Hessenberg 阵，然后用QR 方法求A 的全部特征值．第九章 微分方程初值问题的数值解法9.1 用反复迭代（反复校正）的欧拉预估－校正法求解初值问题0, 0<0.2(0)1y y x y '+=≤⎧⎨=⎩，要求取步长0.1h =，每步迭代误差不超过510-． 答案： [4]11(0.1)0.904 762y y y ≈==，[4]22(0.2)0.818 594y y y ≈==9.2 用二阶中点格式和二阶休恩格式求初值问题2, 0<0.4(0)1dy x y x dx y ⎧=+≤⎪⎨⎪=⎩的数值解(取步长0.2h =，运算过程中保留五位小数)．答案：用二阶中点格式，取初值01y =计算得0n =时，1211.000 00, 1.200 00, (0.2)=1.240 00K K y y ==≈ 1n =时，1221.737 60, 2.298 72, (0.4)=1.699 74K K y y ==≈用二阶休恩格式，取初值01y =计算得0n =时，1211.000 00, 1.266 67, (0.2)=1.240 00K K y y ==≈ 1n =时，1221.737 60, 2.499 18, (0.4)=1.701 76K K y y ==≈9.3 用如下四步四阶阿达姆斯显格式1123(5559379)/24n n n n n n y y h f f f f +---=+-+-求初值问题, (0)1y x y y '=+=在[0,0.5]上的数值解．取步长0.1h =，小数点后保留８位．答案：4(0.4)0.583 640 216y y ≈=，5(0.5) 1.797 421 984y y ≈=． 9.4 为使二阶中点公式1(,(,))22n n n n n n h hy y hf x y f x y +=+++，求解初值问题 , (0)y y y aλλ'=-⎧⎨=⎩为实常数绝对稳定，试求步长h 的大小应受到的限制条件． 答案：2h λ≤．9.5 用如下反复迭代的欧拉预估－校正格式(0)1(1)()111(,)[(,)(,)]2 0,1,2,; 0,1,2,nn n n k k n n n n n n y y hf x y h y y f x y f x y k n +++++⎧=+⎪⎪=++⎨⎪⎪==⎩，求解初值问题sin(), 01(0)1x y e xy x y '⎧=<≤⎨=⎩时，如何选择步长h ，使上述格式关于k 的迭代收敛． 答案：2h e<时上述格式关于k 的迭代是收敛的．9.6 求系数,,,a b c d ，使求解初值问题0(,), ()y f x y y x a '==的如下隐式二步法221()n n n n n y ay h bf cf df +++=+++的误差阶尽可能高，并指出其阶数．答案：系数为142,,33a b d c ====，此时方法的局部截断误差阶最高，为五阶5()O h ．9.7 试用欧拉预估－校正法求解初值问题, (0)=1, 0<0.2()/, (0)2dyxy z y dxx dz x y z z dx⎧=-⎪⎪≤⎨⎪=+=⎪⎩，取步长0.1h =，小数点后至少保留六位．答案：由初值00(0)1, (0)2y y z z ====可计算得110.800 000z 2.050 000y =⎧⎨=⎩ ， 11(0.1)0.801 500(0.1) 2.046 951y y z z ≈=⎧⎨≈=⎩ 220.604 820z 2.090 992y =⎧⎨=⎩ ， 22(0.2)0.604 659(0.2) 2.088 216y y z z ≈=⎧⎨≈=⎩。
工程数学-复变函数与积分变换吉林大学数学学院习题详解
《工程数学-复变函数与积分变换》课后习题详解 吉林大学数学学院 （主编：王忠仁 张静）高等教育出版社 习题一（P12）对任何z ，22z z =是否成立如果是，就给出证明。
如果不是，对哪些z 值才成立解：设z x iy =+，则2222z x y xyi =-+，222z x y =+；若22z z =成立，则有22222x y xyi x y -+=+，即222220x y x yxy ⎧-=+⎨=⎩，解得0y =，即z x =。
所以，对任何z ，22z z =不成立，只对z 为实数时才成立。
求下列各式的值：（1）5)i ； （2）6(1)i +； （3； （4）13(1)i -。
解：（162ii eπ-=，所以555556661)223232())2i i i i e e e i i πππ--⨯-⎛⎫====-=- ⎪⎝⎭（2）因为41ii e π+=，所以63663442(1)288i i i e e e i πππ⨯⎫+====-⎪⎭（3）因为1cos sin i ππ-=+，所以()1622cos sin cossin66k k k w i i ππππππ++==+=+，其中0,1,2,3,4,5k =；即01cossin6622w i i ππ=+=+，1cos sin 22w i i ππ=+=，2551cossin 662w i i ππ=+=+，3771cos sin 662w i i ππ=+=-，433cossin 22w i i ππ=+=-，511111cos sin 662w i i ππ=+=-。
（4）因为1cos()sin()44i i ππ⎤-=-+-⎥⎦，所以11362244(1)2cos sin 33k k k w i i ππππ⎡⎤-+-+⎢⎥=-=+⎢⎥⎢⎥⎣⎦，其中0,1,2k =；即1602cos()sin()1212w i ππ⎡⎤=-+-⎢⎥⎣⎦，161772cos sin1212w i ππ⎡⎤=+⎢⎥⎣⎦，162552cos sin 44w i ππ⎡⎤=+⎢⎥⎣⎦。
计算方法 课后习题答案
l0
(x ( x0
x1)(x x2 ) x1)(x0 x2 )
(x
1)( x 6
2)
l1
(x ( x1
x0 )(x x2 ) x0 )(x1 x2 )
(x
1)( x 2
2)
l2
(x ( x2
x0 )(x x1) x0 )(x2 x1)
L3 x 的最高次项系数是 6，试确定 y1 。
解： l0 (x)
x x1 x0 x1
x x2 x0 x2
x x3 x0 x3
x 0.5 0 0.5
x 1 0 1
x2 02
= x3
7 2
x2
7 2
x 1
l1 ( x)
x x0 x1 x0
1 8(x 0) 3(x 0)(x 1) 114 (x 0)(x 1)(x 2) 11 x3 45 x2 1 x 1
4 42
由求解结果可知： L3 (x) N3(x)
说明插值问题的解存在且唯一。
6. 已知由数据 (0, 0), (0.5, y1), (1,3)和(2, 2) 构造出的 Lagrange 插值多项式
2 02 12 4 23 4 04 14 2 3
1 x2 3x 2 x 4 3x x2 6x 8 23 x x2 5x 4 1 x x2 3x 2
8
4
8
11 x3 45 x2 1 x 1
y0
吉林大学工程数学计算方法第三章习题答案word资料9页
第三章习题答案1.分别用梯形公式、Simpson 公式、Cotes 公式计算积分1,I =⎰并估计误差。
解：1）用梯形公式有： 事实上， 2）Simpson 公式事实上，()()()110.50.510.5410.000030462S E f f f f -⎡+⎤⎛⎫=-++= ⎪⎢⎥⎝⎭⎣⎦⎰3）由Cotes 公式有： 事实上，()0.0000003C E f =2．证明Simpson 公式()2.8具有三次代数精度。
证明：而当()4f x x =时左侧：()()45515b b f x dx x dx b a a a ==-⎰⎰ 右侧：左侧不等于右侧。
所以Simpson 具有三次代数精度. 3.分别用复化梯形公式和复化公式Simpson 计算下列积分.(1)21,804x dx n x =+⎰，（3）,4n =⎰，6,sin 4602=-⎰n d ϕϕπ解：（1）用复化梯形公式有：10188b a h n --===，()()[]12345672128888888102(0.0311280.0615380.0905660.117650.142350.164380.18361)0.20.111416n h T f a f f f f f f f f ⎡⎤⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫=++++++++⎢⎥ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎣⎦=+⨯+++++++=由复化Simpson 公式有： 解：删去 解（3）：,4n =⎰由复化梯形公式有： 由复化Simpson 公式有：（4）解：6,sin 4602=-⎰n d ϕϕπ由复化梯形公式： 由复化Simpson 公式：4．给定求积节点012113,,,424x x x ===试推出计算积分()10f x dx ⎰的插值型求积公式，并写出它的截断误差。
解：考虑到对称性，有20A A =，于是有求积公式由于原式含有3个节点，故它至少有2阶精度。
计算方法课后习题集规范标准答案
习 题 一3.已知函数y =4, 6.25,9x x x ===处的函数值，试通过一个二次插值函解：0120124, 6.25,9;2, 2.5,3y x x x y y y =======由题意 （1） 采用Lagrange插值多项式220()()j j j y L x l x y ==≈=∑27020112012010*********()|()()()()()()()()()()()()(7 6.25)(79)(74)(79)(74)(7 6.25)2 2.532.255 2.25 2.75 2.7552.6484848x y L x x x x x x x x x x x x x y y y x x x x x x x x x x x x ==≈------=++------------=⨯+⨯+⨯⨯-⨯⨯= 其误差为(3)25(3)25(3)2[4,9]2()(7)(74)(7 6.25)(79)3!3()83max |()|40.0117281|(7)|(4.5)(0.01172)0.008796f R f x x f x R ξ--=---==<∴<=又则（2）采用Newton插值多项式2()y N x =≈ 根据题意作差商表：224(7)2(74)()(74)(7 6.25) 2.64848489495N =+⨯-+-⨯-⨯-≈4. 设()()0,1,...,k f x x k n ==，试列出()f x 关于互异节点()0,1,...,i x i n =的Lagrange 插值多项式。
注意到：若1n +个节点()0,1,...,i x i n =互异，则对任意次数n ≤的多项式()f x ，它关于节点()0,1,...,i x i n =满足条件(),0,1,...,i i P x y i n ==的插值多项式()P x 就是它本身。
可见，当k n ≤时幂函数()(0,1,...,)kf x x k n ==关于1n +个节点()0,1,...,i x i n =的插值多项式就是它本身，故依Lagrange 公式有()00(),0,1,...,nn n k kk i j j j j j i j ii jx x x l x x x k n x x ===≠-=≡=-∑∑∏特别地，当0k =时，有()0001nn n ij j j i j ii jx x l x x x ===≠-=≡-∑∑∏而当1k =时有()000nnn ij j j j j i j ii jx x x l x x x x x ===≠⎛⎫- ⎪=≡ ⎪- ⎪⎝⎭∑∑∏ 5.依据下列函数表分别建立次数不超过3的Lagrange 插值多项式和Newton 插值多项式，并验证插值多项式的唯一性。
《工程数学-复变函数与积分变换》吉林大学数学学院 习题详解
《工程数学-复变函数与积分变换》课后习题详解大学数学学院 （主编：王忠仁 静）高等教育 习题一（P12）1.1 对任何z ，22z z =是否成立？如果是，就给出证明。
如果不是，对哪些z 值才成立？解：设z x iy =+，则2222z x y xyi =-+，222z x y =+；若22z z =成立，则有22222x y xyi x y -+=+，即222220x y x yxy ⎧-=+⎨=⎩，解得0y =，即z x =。
所以，对任何z ，22z z =不成立，只对z 为实数时才成立。
1.2 求下列各式的值：（1）5)i ； （2）6(1)i +； （3； （4）13(1)i -。
解：（162ii eπ-=，所以555556661)223232())2i i i i e e e i i πππ--⨯-⎛⎫====-=- ⎪⎝⎭（2）因为41ii e π+=，所以63663442(1)288i i i e e e i πππ⨯⎫+====-⎪⎭（3）因为1cos sin i ππ-=+，所以()1622cos sin cossin66k k k w i i ππππππ++==+=+，其中0,1,2,3,4,5k =；即01cossin6622w i i ππ=+=+，1cos sin 22w i i ππ=+=，2551cossin 662w i i ππ=+=+，3771cos sin 662w i i ππ=+=-，433cossin 22w i i ππ=+=-，511111cos sin 662w i i ππ=+=-。
（4）因为1cos()sin()44i i ππ⎤-=-+-⎥⎦，所以11362244(1)2cos sin 33k k k w i i ππππ⎡⎤-+-+⎢⎥=-=+⎢⎥⎢⎥⎣⎦，其中0,1,2k =；即1602cos()sin()1212w i ππ⎡⎤=-+-⎢⎥⎣⎦，161772cos sin1212w i ππ⎡⎤=+⎢⎥⎣⎦，162552cos sin 44w i ππ⎡⎤=+⎢⎥⎣⎦。
工程数学(线性代数与概率统计)第三章典型例题分析
第三章例1 设A 为n 阶方阵, 若存在正整数k 和向量, 使, 且.证明: 向量组线性无关.证明: （利用线性无关定义证明） 假设有常数, 使得1120k k A A λαλαλα-++= （1）将（1）两边左乘, 可得122120k k k k A A A λαλαλα--++=由已知条件, 可知上式从第二项全等于零, 所以, 又由条件, 所以. 类似地, 将（1）两边左乘, 可得； 类似地可证得，所以向量组线性无关.例2 设向量组线性相关, 向量组线性无关, 问:（1）能否由线性表示？ 证明你的结论; （2）能否由线性表示？ 证明你的结论. 解: （1）能由线性表示.证明：由于向量组线性无关, 那么其部分组也线性无关。
又由已知条件有线性相关, 故能由线性表示. (2) 4α不能由123,,ααα线性表示.证明:假设4α能由123,,ααα线性表示,即存在不全为零的常数123,,λλλ,使得4112233ααλαλαλ=++由(1)的结论,我们可以设12233k k ααα=+,代入上式,可得421223133()()k k αλλαλλα=+++即4α可由23,αα线性表示,从而234,,ααα线性相关,与已知条件矛盾.因此假设不成立, 4α不能由123,,ααα线性表示.例3 设两向量组()()()123(1)1,2,3,3,0,1,9,6,7TTTααα=-==- ()()()123(2)0,1,1,,2,1,,1,0TTTa b βββ===已知两向量组的秩相等,且3β能由123,,ααα线性表示,求a,b. 解:令123123(,,),(,,)A B αααβββ==由于矩阵A 已知, 可以先对A 进行初等变换求秩.12231313913913925206061206123331701020000r r A r r r r ⎛⎫⎛⎫⎛⎫-+ ⎪ ⎪ ⎪=--+-- ⎪ ⎪ ⎪+ ⎪ ⎪ ⎪--⎝⎭⎝⎭⎝⎭因此()2r A =,且12,αα为(1)的一个极大无关组.由已知条件两向量组的秩相等,所以()2r B =,从而0B =,即0121011a bB a b ==-= 所以a b =.又由条件3β能由123,,ααα线性表示而12,αα为(1)的一个极大无关组.所以3β能由12,αα线性表示,则1230ααβ=,即123132012100310b b ααβ⎛⎫⎪==-= ⎪⎪-⎝⎭,解得 5b =,所以有5a b ==.例4求向量组()11,1,1,3,T α=-()21,3,5,1Tα=-,()32,6,10,Ta α=-,()44,1,6,10Tα=-,()53,2,1,Tc α=-的秩和一个极大无关组.解:对以12345,,,,ααααα为列构成的矩阵A,做初等变换112431124313612024311510610612243110046291124311243024310243100077000110028110203A a c a c Ba c a c ----⎡⎤⎡⎤⎢⎥⎢⎥---⎢⎥⎢⎥=⎢⎥⎢⎥--⎢⎥⎢⎥+--⎣⎦⎣⎦----⎡⎤⎡⎤⎢⎥⎢⎥--⎢⎥⎢⎥=⎢⎥⎢⎥--⎢⎥⎢⎥-----⎣⎦⎣⎦当a=2且c=3时, ,B 中第1.2.4列线性无关, 此时向量组的秩为3, 是一个极大无关组；当时, , B 中第1.2.3.4列线性无关, 此时向量组的秩为4, 是一个极大无关组；当, , B 中第1.2、4、5列线性无关此时向量组的秩为4, 是一个极大无关组.例5设向量组（1）的秩为3；向量组（2）的秩为4,证明:向量组的秩为4.证明: （要证明的秩为4, 可通过证明线性无关来得到想要的结论） 由向量组（2）的秩为4, 可知线性无关, 又由向量组（1）的秩为3, 可知线性相关, 从而可由线性表示, 即存在不全为零的常数, 使得, 不妨设, 将代入, 可得14112422343345()()()0k k l k k l k k l k αααα-+-+-+= 由于线性无关, 所以1412421234343400000k k l k k l k k k k k k l k -=⎧⎪-=⎪⇒====⎨-=⎪⎪=⎩故线性无关, 从而该向量组的秩为4.例6 设向量组的秩为r, , , , , 证明向量组12,,,m βββ的秩为r证明:(由推论等价的向量组有相同的秩, 此题只需证明两个向量组等价即可)由已知可由线性表示, 且有下式成立1212(1)()m m m βββααα+++=-+++从而,于是有, 即也可由, 故向量组与向量组等价, 从而他们的秩相等, 从而向量组的秩为r.。
2019-2020学年第一学期期末考试《计算方法》大作业答案
吉林大学网络教育学院2019-2020学年第一学期期末考试《计算方法》大作业答案学生姓名专业层次年级学号学习中心成绩年月日作业完成要求：大作业要求学生手写，提供手写文档的清晰扫描图片，并将图片添加到word文档内，最终wod文档上传平台，不允许学生提交其他格式文件（如JPG，RAR等非word文档格式），如有雷同、抄袭成绩按不及格处理。
一、解线性方程（每小题8分，共80分）1、用矩阵的LU分解算法求解线性方程组X1+2X2+3X3= 02X1+2X2+8X3= -4-3X1-10X2-2X3= -11答：2、用矩阵的Doolittle分解算法求解线性方程组X1+2X2+3X3= 12X1– X2+9X3= 0-3X1+ 4X2+9X3= 1答：3、用矩阵的Doolittle分解算法求解线性方程组2X1+X2+X3= 46X1+4X2+5X3=154X1+3X2+6X3= 13答：4、用高斯消去法求解线性方程组2X1- X2+3X3= 24X1+2X2+5X3= 4-3X1+4X2-3X3= -3答：5、用无回代过程消元法求解线性方程组2X1- X2+3X3= 24X1+2X2+5X3= 4-3X1+4X2-3X3= -3答：6、用主元素消元法求解线性方程组2X1- X2+3X3= 24X1+2X2+5X3= 4-3X1+4X2-3X3= -3答：7、用高斯消去法求解线性方程组1231231232344272266x x x x x x x x x -+=++=-++=答：8、利用Doolittle 分解法解方程组Ax=b ，即解方程组12341231521917334319174262113x x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥--⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎣⎦ 答：9、利用Doolittle 分解法解方程组Ax=b ，即解方程组123421111443306776081011112x x x x ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦ 答：10、用高斯消元法解方程组1237811351341231x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-=-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦答案：二、计算（每小题10分，共20分）1、已知节点x1,x2及节点处函数值f(x1)，f(x2)，构造线性插值多项式p1(x). 答：2、设f(xi)=i(i=0,1,2),构造二次式p2(x)，使满足： p2(xi)=f(xi)(i=0,1,2)答：。
吉林大学工程数学计算方法第三章习题答案共7页word资料
第三章习题答案11. 分别用梯形公式、Simpson公式、Cotes公式计算积分I xdx,并估计误差。
0.5解：1)用梯形公式有：事实上，2) Simps on 公式事实上，E S(f) —[T Xdx —(0.5)+4f l f(1)]= 0.00003043) 由Cotes公式有：事实上，E C f = 0.00000032. 证明Simpson公式2.8具有三次代数精度。
证明: 而当f x =x4时b b 4 i 5 5左侧：f x dx x4dx b-a5'a 'a 5右侧：左侧不等于右侧。
所以Simps on具有三次代数精度.3•分别用复化梯形公式和复化公式Simps on计算下列积分.1 9 M r _____________________(1) dx, n=8 , (3) . xdx,n=4 , 仏4-sin2'd n = 604 x2 1 0解：(1)用复化梯形公式有:b - a 1_ 0 1h = ■n 8 8T n f a 210 2 (0.031128 0.061538 0.090566 0.11765 0.14235 0.16438 0.18361) 0.21 - 0.111416由复化Simpson公式有:解：删去解(3)：存dx,® 由复化梯形公式有: 由复化Simps on公式有:(4)解：\4_sin 2 d , n = 6 由复化梯形公式： 由复化 Simps on 公式：113 14•给定求积节点x o ,x ),X 2,试推出计算积分.0f x dx 的插值型求积公式， 并写出它的截断误差。
解:考虑到对称性，有 氏=氏，于是有求积公式由于原式含有 3个节点，故它至少有 2阶精度。
考虑到其对称性，可以猜想到它可能有 3阶精度。
事实上，对 f =X 原式左右两端相等：此外，容易验证原式对 f =X 4不准确，故所构造出的求积公式有 3阶精度。
计算方法 课后习题答案
计算方法课后习题答案计算方法课后习题答案计算方法是一门重要的学科，它涉及到数值计算、算法设计和数据处理等方面的内容。
在学习计算方法的过程中，课后习题是不可或缺的一部分。
通过解答习题，我们可以巩固所学的知识，提高自己的计算能力。
下面是一些计算方法课后习题的答案，希望对大家的学习有所帮助。
1. 矩阵的转置矩阵的转置是将矩阵的行和列互换得到的新矩阵。
对于一个m×n的矩阵A，它的转置记作A^T。
转置后的矩阵A^T的行数和列数分别为原矩阵A的列数和行数。
例如，对于一个3×2的矩阵A，它的转置A^T是一个2×3的矩阵。
2. 矩阵的加法和减法矩阵的加法和减法是对应位置上的元素进行相加或相减得到的新矩阵。
对于两个相同大小的矩阵A和B，它们的和记作A+B，差记作A-B。
加法和减法的运算规则是相同位置上的元素进行相应的运算。
3. 矩阵的乘法矩阵的乘法是指将两个矩阵相乘得到一个新矩阵的运算。
对于两个矩阵A和B，它们的乘积记作AB。
矩阵乘法的运算规则是矩阵A的行与矩阵B的列进行相乘，并将结果相加得到新矩阵的对应位置上的元素。
4. 矩阵的逆矩阵的逆是指对于一个可逆矩阵A，存在一个矩阵B，使得AB=BA=I，其中I是单位矩阵。
如果一个矩阵A存在逆矩阵，则称其为可逆矩阵或非奇异矩阵。
求解矩阵的逆可以使用伴随矩阵和行列式的方法。
5. 线性方程组的求解线性方程组是指由一组线性方程组成的方程组。
求解线性方程组的方法有很多，包括高斯消元法、LU分解法、迭代法等。
其中，高斯消元法是一种常用的求解线性方程组的方法，它通过消元和回代的过程，将线性方程组转化为上三角形矩阵或对角矩阵，从而求解出方程组的解。
6. 数值积分的方法数值积分是指通过数值计算的方法来求解定积分的近似值。
常用的数值积分方法包括梯形法则、辛普森法则和龙贝格法则等。
这些方法都是基于将定积分转化为离散求和的形式，通过计算离散点上的函数值来估计定积分的近似值。
精选吉林大学工程数学计算方法第三章习题答案资料
第三章习题答案1.分别用梯形公式、Simpson公式、Cotes公式计算积分1,I=⎰并估计误差。
解：1）用梯形公式有：()()110.51[10.5]10.42678242f f⎛-≈+=+≈⎝⎭⎰()()()333333220.512.6042107.36571012124Tb aE f fηηη-----⎛⎫''=-=--=⨯≤⨯⎪⎝⎭事实上，()()()()()()110.430964410.50.510.4267767210.50.510.00418772Tf x II f fE f f f===-≈+=⎡⎤⎣⎦-∴=-+=⎡⎤⎣⎦⎰⎰2）Simpson公式()110.53111410.43093 64212f f f⎛-⎡⎤⎛⎫⎛⎫≈++=+=⎪ ⎪⎢⎥⎝⎭⎝⎭⎣⎦⎝⎭⎰[]()()44744211111522 1.1837710180218028Sb a b aE f fηη--⎛⎫--⎪⎛⎫--⎛⎫=-=--≤⨯⎪ ⎪⎪⎝⎭⎝⎭⎪⎝⎭()312()''()48T f fb aE h=?--事实上，()()()110.50.510.5410.000030462SE f f f f-⎡+⎤⎛⎫=-++=⎪⎢⎥⎝⎭⎣⎦⎰3）由Cotes公式有：()() ()111537270.5321232719084814.9497525.2982210.3923029.9332670.43096180f f f f f-⎡⎤⎛⎫⎛⎫⎛⎫≈++++⎪ ⎪ ⎪⎢⎥⎝⎭⎝⎭⎝⎭⎣⎦=++++=⎰17)180+++()6116211294522 2.697410945464C E f η--⎛⎫⨯ ⎪⎛⎫=-⨯-≤⨯ ⎪ ⎪⎝⎭⎪⎝⎭()7(6)945*4()()82Cf b a E f h =?-- 事实上，()0.0000003C E f =2．证明Simpson 公式()2.8具有三次代数精度。
工程数学(线性代数与概率统计)答案(3章)
工程数学（线性代数与概率统计）习题三1、2、3、略4、)1,0,1()1,1,0()0,1,1(21-=-=-αα)2,1,0()0,4,3()1,1,0(2)0,1,1(323321=-+=-+ααα5、)523(61)(5)(2)(3321321αααααααααα-+=→+=++-6、设存在一组数r k k k ,,,21 使得 0)()()()(02212121212112211=++++++++=+++++++==+++r r r r r r r r k k k k k k k k k k k k αααααααααβββ因r ααα ,,21线性无关，有⎪⎪⎩⎪⎪⎨⎧==++=+++000221rr r k k k k k k 即021====r k k k ，所以r βββ ,,21线性无关。
7、设存在一组数4321,,,k k k k 使得044332211=+++ββββk k k k 有0)()()()(443332221141=+++++++ααααk k k k k k k k 因0000000043322141=k k k k k k k k ，且不全为0，所以4321,,,ββββ线性相关。
8、讨论向量组相关性。
（本题的特点是向量组的个数等于向量的维数， 其判断法是求向量组成的行列式值是否为0）（1）052520111631520111321===ααα，相关 （2）02102011321≠==ααα，无关 9、由向量组组成的行列式为 1221011131321111321-==t tααα（1）如果,5,41=→=-t t 行列式等于0，向量组线性相关， （2）如果,5,41≠→≠-t t 行列式不等于0，向量组线性无关， （3）当5=t 时，向量组相关，设22113αααk k += 即⎩⎨⎧=-=⇒⎪⎪⎪⎭⎫ ⎝⎛+⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛213211115312121k k k k 10、用矩阵的秩判别向量组的相关性（方法是求由向量组构成的矩阵的秩r 与向量组个数关系） （1）()⎪⎪⎪⎪⎪⎭⎫⎝⎛---→⎪⎪⎪⎪⎪⎭⎫⎝⎛----−−→−⎪⎪⎪⎪⎪⎭⎫⎝⎛---==--01502601401051562641401041242031111323213321c c c c A ααα所以 2)(=A R ，相关。
工程数学测试题及答案3
第三章 复变函数的积分一、选择题：1．设c 为从原点沿x y =2至i +1的弧段，则=+⎰cdz iy x )(2( )（A ）i 6561- （B ）i 6561+- （C ）i 6561-- （D ）i 6561+ 2．设c 为不经过点1与1-的正向简单闭曲线，则dz z z zc⎰+-2)1)(1(为( ) （A ）2i π （B ）2iπ- （C ）0 （D ）(A)(B)(C)都有可能 3．设1:1=z c 为负向，3:2=z c 正向，则=⎰+=dz z zc c c 212sin ( ) （A ） i π2- （B ）0 （C ）i π2 （D ）i π4 4．设c 为正向圆周2=z ，则=-⎰dz z zc 2)1(cos ( ) （A ）1sin - （B ）1sin （C ）1sin 2i π- （D ）1sin 2i π5．设c 为正向圆周21=z ，则=--⎰dz z z z c23)1(21cos( )（A ）)1sin 1cos 3(2-i π （B ）0 （C ）1cos 6i π （D ）1sin 2i π-6．设ξξξξd ze zf ⎰=-=4)(，其中4≠z ，则='）i f π(( ) （A ）i π2- （B ）1- （C ）i π2 （D ）17．设)(z f 在单连通域B 内处处解析且不为零，c 为B 内任何一条简单闭曲线，则积分dz z f z f z f z f c⎰+'+'')()()(2)( ( )（A ）于i π2 （B ）等于i π2- （C ）等于0 （D ）不能确定8．设c 是从0到i 21π+的直线段，则积分=⎰cz dz ze （ ）（A ）21eπ-(B) 21eπ-- (C)i e21π+(D) i e21π-9．设c 为正向圆周0222=-+x y x ，则=-⎰dz z z c1)4sin(2π( ) （A ）i π22 （B ）i π2 （C ）0 （D ）i π22- 10．设c 为正向圆周i a i z ≠=-,1，则=-⎰c dz i a zz 2)(cos ( ) （A ）ie π2 （B ）eiπ2 （C ）0 （D ）i i cos 11．设)(z f 在区域D 内解析，c 为D 内任一条正向简单闭曲线，它的内部全属于D ．如果)(z f 在c 上的值为2，那么对c 内任一点0z ,)(0z f ( )（A ）等于0 （B ）等于1 （C ）等于2 （D ）不能确定 12．下列命题中，不正确的是( ) （A ）积分⎰=--ra z dz az 1的值与半径)0(>r r 的大小无关 （B ）2)(22≤+⎰cdz iy x ,其中c 为连接i -到i 的线段 （C ）若在区域D 内有)()(z g z f ='，则在D 内)(z g '存在且解析 （D ）若)(z f 在10<<z 内解析，且沿任何圆周)10(:<<=r r z c 的积分等于零，则)(z f 在0=z 处解析13．设c 为任意实常数，那么由调和函数22y x u -=确定的解析函数iv u z f +=)(是 ( )(A)c iz +2（B ） ic iz +2（C ）c z +2（D ）ic z +214．下列命题中，正确的是( )（A ）设21,v v 在区域D 内均为u 的共轭调和函数，则必有21v v = （B ）解析函数的实部是虚部的共轭调和函数 （C ）若iv u z f +=)(在区域D 内解析，则xu∂∂为D 内的调和函数 （D ）以调和函数为实部与虚部的函数是解析函数15．设),(y x v 在区域D 内为),(y x u 的共轭调和函数，则下列函数中为D 内解析函数的是( )（A ）),(),(y x iu y x v + （B ）),(),(y x iu y x v -（C ）),(),(y x iv y x u - （D ）xv i x u ∂∂-∂∂二、填空题1．设c 为沿原点0=z 到点i z +=1的直线段，则=⎰cdz z 22．设c 为正向圆周14=-z ，则=-+-⎰c dz z z z 22)4(233．设⎰=-=2)2sin()(ξξξξπd zz f ,其中2≠z ，则=')3(f 4．设c 为正向圆周3=z ，则=+⎰cdz zzz 5．设c 为负向圆周4=z ，则=-⎰c zdz i z e 5)(π6．解析函数在圆心处的值等于它在圆周上的 7．设)(z f 在单连通域B 内连续，且对于B 内任何一条简单闭曲线c 都有0)(=⎰cdz z f ，那么)(z f 在B 内8．调和函数xy y x =),(ϕ的共轭调和函数为9．若函数23),(axy x y x u +=为某一解析函数的虚部，则常数=a10．设),(y x u 的共轭调和函数为),(y x v ，那么),(y x v 的共轭调和函数为 三、计算积分 1.⎰=+-R z dz z z z)2)(1(62,其中1,0≠>R R 且2≠R ; 2.⎰=++22422z z z dz． 四、设)(z f 在单连通域B 内解析，且满足)(1)(1B x z f ∈<-.试证１．在B 内处处有0)(≠z f ； ２．对于B 内任意一条闭曲线c ，都有0)()(=''⎰cdz z f z f 五、设)(z f 在圆域R a z <-内解析，若)0()()(max R r r M z f ra z <<==-，则),2,1()(!)()( =≤n r r M n a fnn .六、求积分⎰=1z zdz z e ，从而证明πθθπθ=⎰0cos )cos(sin d e . 七、设)(z f 在复平面上处处解析且有界，对于任意给定的两个复数b a ,，试求极限⎰=+∞→--R z R dz b z a z z f ))(()(lim并由此推证)()(b f a f =（刘维尔Liouville 定理）.八、设)(z f 在)1(><R R z 内解析，且2)0(,1)0(='=f f ，试计算积分⎰=+122)()1(z dz z z f z 并由此得出⎰πθθθ202)(2cos d e f i 之值.九、设iv u z f +=)(是z 的解析函数，证明222222222))(1()(4))(1ln())(1ln(z f z f yz f xz f +'=∂+∂+∂+∂.十、若)(22y x u u +=，试求解析函数iv u z f +=)(.答案第三章 复变函数的积分一、1．（D） 2．（D） 3．（B） 4．（C）５．（B）６．（A ） ７．（C ） ８．（A ） ９．（A ） 10．（C ） 11．（C ） 12．（D ） 13．（D ） 14．（C ） 15．（B ）二、1．2 2．i π10 3．0 4．i π6 5．12iπ 6．平均值 7．解析 8．C x y +-)(21229．3- 10．),(y x u - 三、1．当10<<R 时，0； 当21<<R 时，i π8； 当+∞<<R 2时，0.2．0. 六、i π2. 七、0.八、,8)()1(122i dz z z f z z π=+⎰=π=θθ⎰πθ2)(2cos 202d e f i . 十、321ln 2)(ic c z c z f ++=（321,,c c c 为任意实常数）.（注：可编辑下载，若有不当之处，请指正，谢谢!）。
计算方法第三章习题答案
第三章习题解答1.试讨论a 取什么值时，下列线性方程组有解，并求出解 。
123123123123212312311(1)1(2)1ax x x ax x x x ax x x ax x a x x ax x x ax a⎧++=++=⎧⎪⎪++=++=⎨⎨⎪⎪++=++=⎩⎩ 解：(1)111111111a A a a ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦ 经初等行变换化为1001/(2)0101/(2)0011/(2)a a a +⎡⎤⎢⎥+⎢⎥⎢⎥+⎣⎦当2a ≠-时，方程组有解，解为111(,,).222Tx a a a =+++ (2)21111111a A a a a a ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦ 经初等行变换化为2100(1)/(2)0101/(2)001(21)/(2)a a a a a a -++⎡⎤⎢⎥+⎢⎥⎢⎥+++⎣⎦当2a ≠-时，方程组有解，解为21121(,,).222Ta a a x a a a +++=-+++2.证明下列方程组Ax=b12341123421233234432432385x x x x b x x x x b x x x b x x x b+--=⎧⎪-+-=⎪⎨+-=⎪⎪-+-=⎩ 当（1）(10,4,16,3).T b =-时无解；（2）(2,3,1,3).T b =时有无穷多组解。
解：(1) r(A)=3≠r(A,b)=4 当(10,4,16,3).T b =-时无解；(2) r(A)=3,r(A,b)=3 当(2,3,1,3).T b =时有无穷多组解。
3.用列主元高斯消元法求解Ax=b2233(1)477,12457A b ⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦ 1231(2)234,13462A b ⎡⎤⎡⎤⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(1)x=(2,-2,1)T (2)x=(0,-7,5)T4.证明上（下）三角方阵的逆矩阵任是上（下）三角方阵。
2020年奥鹏吉大网络教育《计算方法》大作业解答
2020年奥鹏吉⼤⽹络教育《计算⽅法》⼤作业解答2020年奥鹏吉⼤⽹络教育《计算⽅法》⼤作业解答(说明：前⾯是题⽬，后⾯⼏页是答案完整解答部分，注意的顺序。
)⼀、解线性⽅程⽤矩阵的LU分解算法求解线性⽅程组⽤矩阵的Doolittle分解算法求解线性⽅程组⽤矩阵的Doolittle分解算法求解线性⽅程组⽤⾼斯消去法求解线性⽅程组⽤⾼斯消去法求解线性⽅程组⽤主元素消元法求解线性⽅程组⽤⾼斯消去法求解线性⽅程组利⽤Doolittle分解法解⽅程组Ax=b，即解⽅程组1、⽤矩阵的LU分解算法求解线性⽅程组X1+2X2+3X3 = 02X1+2X2+8X3 = -4-3X1-10X2-2X3 = -112、⽤矩阵的Doolittle分解算法求解线性⽅程组X1+2X2+3X3 = 12X1– X2+9X3 = 0-3X1+ 4X2+9X3 = 13、⽤矩阵的Doolittle分解算法求解线性⽅程组2X1+X2+X3 = 46X1+4X2+5X3 =154X1+3X2+6X3 = 134、⽤⾼斯消去法求解线性⽅程组2X 1- X 2+3X 3 = 24X 1+2X 2+5X 3 = 4-3X 1+4X 2-3X 3 = -35、⽤⽆回代过程消元法求解线性⽅程组2X 1- X 2+3X 3 = 24X 1+2X 2+5X 3 = 4-3X 1+4X 2-3X 3 = -36、⽤主元素消元法求解线性⽅程组2X 1- X 2+3X 3 = 24X 1+2X 2+5X 3 = 4-3X 1+4X 2-3X 3 = -37、⽤⾼斯消去法求解线性⽅程组1231231232344272266xx x x x x x x x -+=++=-++=8、利⽤Doolittle 分解法解⽅程组Ax=b ，即解⽅程组 12341231521917334319174262113x x x x -?-=?--??--9、利⽤Doolittle 分解法解⽅程组Ax=b ，即解⽅程组 123421111443306776081011112x x x x =-????10、⽤⾼斯消元法解⽅程组1237811351341231x x x --=-??????⼆、计算1、已知节点x1,x2及节点处函数值f(x1)，f(x2)，构造线性插值多项式p1(x).2、设f(xi)=i(i=0,1,2),构造⼆次式p2(x)，使满⾜： p2(xi)=f(xi)(i=0,1,2)以下为解答部分：⼀、解线性⽅程1、⽤矩阵的LU 分解算法求解线性⽅程组X 1+2X 2+3X 3 = 02X 1+2X 2+8X 3 = -4-3X1-10X2-2X3 = -11解答：2、⽤矩阵的Doolittle分解算法求解线性⽅程组 X1+2X2+3X3 = 12X1– X2+9X3 = 0-3X1+ 4X2+9X3 = 1解答：3、⽤矩阵的Doolittle分解算法求解线性⽅程组 2X1+X2+X3 = 46X1+4X2+5X3 =154X1+3X2+6X3 = 13解答：4、⽤⾼斯消去法求解线性⽅程组 2X1- X2+3X3 = 24X1+2X2+5X3 = 4-3X1+4X2-3X3 = -3解：⽅程组的扩⼤矩阵为：5、⽤⽆回代过程消元法求解线性⽅程组 2X1- X2+3X3 = 2 4X1+2X2+5X3 = 4-3X1+4X2-3X3 = -3解：⽅程组的扩⼤矩阵为：接下⼀页：。
吉林大学2021年《计算方法》 答案
计算方法一计算题 (共10题，总分值100分 )1. 用尤拉法解初值问题取步长h=0.1计算。
（10 分）因为y'=x2+100y2即f(x，y)=x2+100y2，因为欧拉法公式为y n+1=y n+hf(x n，y n)取h=0.1,x0=0,y(x0)=y0=0f(x0,y0)=0所以y1=y0+0.1f(x0，y0)=0f(x1,y1)=f(0.1,0)=0.01y2=y1+0.1f(x1，y1)=0.001f(x2，y2)=f(0.2，0.001)=0.0401y3=y2+hf(x2,y2)=0.005012. 已知函数表：用Simpson公式求的近似值。
（10 分）3. 基于迭代原理证明（10 分）证明：令x=√(2+√(2+......x²=2+√(2+√(2+......x²-2=√(2+√(2+......x²-2=xx²-x-2=0(x+1)(x-2)=0x1=-1（不符合题义，舍去）x2=24. 求用高斯－塞德尔迭代求解线性代数方程组的两次迭代解（取初始向量X(0)=0）。
（10 分）5. 利用Doolittle分解法解方程组Ax=b，即解方程组（10 分）6. 试证明Euler显格式是一阶方法。
（10 分）7. 用高斯消元法解方程组（10 分）8. 下列矩阵矩阵能否分解为LU（其中L为单位下三角阵，U为上三角阵）?若能分解是否唯一？（10 分）A为A9.（10 分）解：根据插值公式构造三次式 3222(1)(2)(3)()(01)(02)(03)(0)(1)(3)(0)(1)(2)(7)26(0)(1)(3)(30)(31)(32)x x x p x x x x x x x x x x ---=+----------+------ 10. 用高斯消去法求解线性方程组 2X 1- X 2+3X 3 = 24X 1+2X 2+5X 3 = 4-3X 1+4X 2-3X 3 = -3 （10 分）。
奥鹏吉大17秋学期《计算方法》在线作业一含答案100分.doc
A.一^+毛=2B. -2牛+1.5也=3.5
己知A =
q 2■ ■
,则化为A为对角阵的平面旋转变换角。=（）°
A.
B.-
7
4
C.
D.-
I
2.
A.A
B.BC.CD.D源自正确答案：4.所谓（）插值，就是将被插值函数逐段多项式化。
A.牛顿
B.拉格朗日
C.三次样条
D.分段
正确答案：
5.题面如下图所示，正确的是（）
玉+2沔+毛=0
.用列主元消元法解线性方程组2毛+2些+3忍=3作第一次消元后得到的第3个方程（）.
使用方法：为保护知识产权，本文档中的正确答案：被设置
奥鹏吉大17秋学期《计算方法》在线作业一含答案100分
一、单选题（共15道试题，共60分。）
1.数值3.1416的有效位数为（）
A.3
B.4
C.5
D.6
正确答案：C
2.常用的折线函数是简单（）次样条函数
A.零
B.—
C.二
D.三
正确答案：
3.题面如下，正确的是（）
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第三章习题答案1.分别用梯形公式、Simpson公式、Cotes公式计算积分1,I=⎰并估计误差。
解：1）用梯形公式有：()()110.51[10.5]10.42678242f f⎛-≈+=+≈⎝⎭⎰()()()333333220.512.6042107.36571012124Tb aE f fηηη-----⎛⎫''=-=--=⨯≤⨯⎪⎝⎭事实上，()()()()()()110.430964410.50.510.4267767210.50.510.00418772Tf x II f fE f f f===-≈+=⎡⎤⎣⎦-∴=-+=⎡⎤⎣⎦⎰⎰2）Simpson公式()110.53111410.43093 642122f f f⎛-⎡⎤⎛⎫⎛⎫≈++=+=⎪ ⎪⎢⎥⎝⎭⎝⎭⎣⎦⎝⎭⎰[]()()44744211111522 1.1837710180218028Sb a b aE f fηη--⎛⎫--⎪⎛⎫--⎛⎫=-=--≤⨯⎪ ⎪⎪⎝⎭⎝⎭⎪⎝⎭3122()''()48T f fb aE事实上，()()()10.510.50.510.5410.000030462SE f f f f-⎡+⎤⎛⎫=-++=⎪⎢⎥⎝⎭⎣⎦⎰3）由Cotes公式有：()() ()111537270.5321232719084814.9497525.2982210.3923029.9332670.43096180f f f f f-⎡⎤⎛⎫⎛⎫⎛⎫≈++++⎪ ⎪ ⎪⎢⎥⎝⎭⎝⎭⎝⎭⎣⎦=++++=⎰15732127)18088()6116211294522 2.697410945464C E f η--⎛⎫⨯ ⎪⎛⎫=-⨯-≤⨯ ⎪ ⎪⎝⎭⎪⎝⎭7(6)945*42()()82Cf b aEf事实上，()0.0000003C E f =2．证明Simpson 公式()2.8具有三次代数精度。
证明：()()()()()()333344444224446243baa b a b f x x f a a f f b bxb a f x dx b a a b b a f a f f b ++⎛⎫⎛⎫==== ⎪ ⎪⎝⎭⎝⎭-===-⎡+⎤-⎛⎫=++=⎪⎢⎥⎝⎭⎣⎦⎰b a令，则，，左边右边故该公式的代数精度是。
而当()4f x x =时左侧：()()45515b b f x dx x dx b a a a ==-⎰⎰ 右侧：()()()44455432234446268552232a b b a a b b a f a f f b a b b a a b a b a b ab ⎡⎤+-⎡+⎤-⎛⎫++=+⨯+⎢⎥ ⎪⎢⎥⎝⎭⎣⎦⎢⎥⎣⎦---+-=左侧不等于右侧。
所以Simpson 具有三次代数精度.3.分别用复化梯形公式和复化公式Simpson 计算下列积分.(1)21,804x dx n x =+⎰，（3）,4n =⎰，6,sin 4602=-⎰n d ϕϕπ解：（1）用复化梯形公式有：10188b a h n --===，()()[]12345672128888888102(0.0311280.0615380.0905660.117650.142350.164380.18361)0.20.111416n h T f a f f f f f f f f ⎡⎤⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫=++++++++⎢⎥ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎣⎦=+⨯+++++++=由复化Simpson 公式有：()()()()811123135702()146444488881020.0615380.117650.1643840.0311280.0905660.412350.183510.2240.11157S f f f f f f f f f ⎡⎤⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫=⨯+⨯+++++++⎢⎥⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎣⎦=+⨯+++⨯++++⎡⎤⎣⎦=()()12112,100x e dx n x--=⎰解：删去 解（3）：,4n =⎰由复化梯形公式有：()()()()()()()()4912,412192357213217.2277b a h n T f f f f f --====⨯++++=++⨯≈由复化Simpson 公式有：()()()()()()()()414192543762132417.32203S f f f f f =⨯++++=⨯++≈（4）解：6,sin 4602=-⎰n d ϕϕπ由复化梯形公式：0356219.1]36362)0([36)]()(2)([25,4,3,2,1,,3660651516=⎪⎭⎫⎝⎛+⎪⎭⎫ ⎝⎛+=++==+==-=-=∑∑==πππϕϕππf k f f b f f a f h T k kh a n a b h k k k k由复化Simpson 公式：035763886.1,035834878.13672365,4,3,2,1,0,2,,3231456215216664==⎪⎭⎫⎝⎛+==+=⎪⎪⎭⎫ ⎝⎛=+=∑∑=+=+S k f H k hf h H H T S k kk k k πππϕϕϕ4．给定求积节点012113,,,424x x x ===试推出计算积分()10f x dx ⎰的插值型求积公式，并写出它的截断误差。
解：()10120113424f x dx A f A fA f ⎛⎫⎛⎫⎛⎫≈++ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎰10011013224,11133424413144111332424x x A dx x x A dx ⎛⎫⎛⎫-- ⎪⎪⎝⎭⎝⎭==⎛⎫⎛⎫-- ⎪⎪⎝⎭⎝⎭⎛⎫⎛⎫-- ⎪⎪⎝⎭⎝⎭==-⎛⎫⎛⎫-- ⎪⎪⎝⎭⎝⎭⎰⎰考虑到对称性，有20A A =，于是有求积公式()121311[]34432f x dx f f f ⎛⎫⎛⎫⎛⎫≈+- ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎰由于原式含有3个节点，故它至少有2阶精度。
考虑到其对称性，可以猜想到它可能有3阶精度。
事实上，对3f x =原式左右两端相等：333130213111[]344324x dx ⎛⎫⎛⎫⎛⎫+-== ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎰ 此外，容易验证原式对4f x =不准确，故所构造出的求积公式有3阶精度。
5．给定积分2sin I xdx π=⎰。
（1） 利用复化梯形公式计算上述积分值，使其截断误差不超过3110;2-⨯ （2） 取同样的求积节点，改用复化Simpson 公式计算时，截断误差是多少？ （3） 如果要求截断误差不超过610-，那么使用复化Simpson 公式计算时，应将积分区间分成多少等分？解：(1) 33''''22()()()()1296nTb a E f f f n n πηη-=-=- ()f x =sin x ,'''()cos ,()sin f x x f x x ==-∴3322()sin ,0,96962n TE f nn πππφη⎡⎤=≤∈⎢⎥⎣⎦当误差3()0.510nT E f -≤⨯时，n ≥25.6, 所以取n =26。
25h :h=[(0)()2()]52221T f f f x n k k ππ⇒=++∑=则12325{012[sin()sin()sin()...sin()]}25252525252πππππ=⨯++++++0.9465= （2）1S 4''''42E []()()()sin()n 180218022nh f f ππηη=-=-⨯⨯b-a 11S 44922E []()(26)()710n 18022n 18022n f n ππππ-≤⨯⨯=⇒⨯⨯=⨯则 1S 462(3)E []()10n 18022nf ππ-≤⨯⨯≤ 7.6=8n n ≥⇒则6．用Romberg 求积方法计算下列积分，使误差不超过510-。
（11x e dx -⎰；（2）20sin x xdx π⎰；（3）3⎰;(4)12041dx x +⎰解（1）：dx ex⎰-12π112111212422()[0,1]:(0)(1)]0.771743332()[01]:110.68439656,()0.728069946,2210.7135121533()[01]:131]0.705895578,()0.442a T f f b H f T T H S T T c H f f T T H =+=⎛⎫===+= ⎪⎝⎭=-=⎛⎫⎛⎫=+==+= ⎪ ⎪⎝⎭⎝⎭在上用梯形公式，二等分，四等分716982762,101052.3,713271669.0141144,713271674.0141144713272634.03134,714200166.0)(21711417571.0]87858381[412]10[)(713272026.0141144,713287034.0313457111323312242224844484122221142--<⨯=-=---==---==-==+==⎪⎭⎫⎝⎛+⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛==---==-=C R C C R S S C T T S H T T f f f f H d S S C T T S π八等分：，将计算可以停止。
解（2）：dx x x ⎰π20sin956833201.5)(21,9788642.6242:]20[)(018385352.7141144,579736267.63134934802201.4)(21,869604401.9]232[:]20[)(03134,0)(21,0)(2:]20[)(0)]2()0([22:]20[)(448341222212422242212111211-=+=-=⎪⎭⎫⎝⎛+==---=-=-=-=+=-=-=⎪⎭⎫⎝⎛+⎪⎭⎫ ⎝⎛==-==+====+=∑=H T T i f H ，d S S C T T S H T T f f H ，c T T S H T T f H ，b f f T ，a i πππππππππππππππ八等分将四等分将二等分将上用梯形公式得在283132311.6141144,284030929.63134202231497.6)(21,447629792.6484]20[)(266954014.6141144278695129.6141144,2975102.6313442822481688816708132331224222484-=---=-=-=-=+=-=⎪⎭⎫⎝⎛+=-=---=-=---=-=-=∑=S S C T T S H T T i f H ，e C C R S S C T T S i ππππ十六等分将581116266125455244844483163381623222163264323232643103215255124444243833482162281632161616321501614244123433210105.9Y -Z 283185304.6Y 141Y 144Z ,283185304.6X 141X 144Y 283185304.6R 141R 144,283185304.6141144283185292.6141144,283188551.63134278137899.6)(21,293289853.6163216:]20[)(283185209.6141_144283185288.6141144,283185356.6141144283184528.6141144,283237428.63134262985945.6)(21,323740394.68168]20[)(283266463.6141144,283202742.6141144--==<⨯=-=---=-=---=-=---=-=---=-=---=-=-=-=+=-=⎪⎭⎫⎝⎛+=-=--=-=---=-=---=-=---=-=-=-=+=-=⎪⎭⎫⎝⎛+=-=---=-=---=∑∑X C C R S S C T T S H T T i f H ，g X X Y R R X C C R S S C T T S H T T i f H ，f R R X C C R i i ππππππππ六十四等分将三十二等分将（3）解:⎰+321dx x x20762073.10141144,20722396.10313426636719.10][21,08893752.10438343:]3,0[)(20457443.10141144,20127249.10313444379685.10)[21,71622377.9]4943[23:]3,0[)(1517434.103134,17136992.11)(21,11249037.8233:]3,0[)(23024947.14)]0()3([23]30[)(22422248444830412222124222421211211=---==-==+==⎪⎭⎫ ⎝⎛+==---==-==+==⎪⎭⎫⎝⎛+⎪⎭⎫ ⎝⎛==-==+==⎪⎭⎫⎝⎛==+=∑=S S C T T S H T T k f H d S S C T T S H T T f f H c T T S H T T f H b f f T ，a k 八等分将四等分将二等分将上用梯形公式在计算可以停止三十二等分将十六等分将56111525512444424383348216228163216161632150161424412343324282248168708816813233110104.1,20759219.1014114420759219.10141144,20759219.1014114420759223.10141144,20759091.10313421126074.10)(21,20025127.10163323163:]30[)(20759364.10141144,20759393.1014114420759435.10141144,2075712.1031342222702.10)(21,1781732.108316383:]30[)(20766908.101411442--==<⨯=-=---==---==---==---==-==+==⎪⎭⎫⎝⎛+==---==---==---==-==+==⎪⎭⎫ ⎝⎛+==---=∑∑X Y X X Y R R X C C R S S C T T S H T T i f H ，f R R X C C R S S C T T S H T T i f H ，e C C R i i 解（4）:dx x ⎰+10214141585784.3141144141594094.3141144,141592502.33134148988495.3)(21,146800518.348141:]10[)(142117647.3141144,141568627.33134131176471.3)(21,162352941.3]4341[21:]10[)(133333333.33134,1.3)(21,2.321:]10[)(3)]0()1([21:]10[)(13233122422248444834122221242224212111211=---==---==-==+==⎪⎭⎫⎝⎛+==---==-==+==⎪⎭⎫⎝⎛+⎪⎭⎫ ⎝⎛==-==+==⎪⎭⎫⎝⎛==+=∑=C C R S S C T T S H T T i f H ，d S S C T T S H T T f f H ，c T T S H T T f H ，b f f T ，a i 八等分将四等分将二等分将上用梯形公式在7816880281684842234242121334411()[01]:1113.14289473,() 3.1409416138168241413.141592652, 3.14159266233414141413.141592639, 3.141592666414141416i e i H f T T H S T T C S S R C C X R R X R =⎛⎫=+==+= ⎪⎝⎭=-==-=--=-==-=-----=∑，十六等分65.881010,--⨯<算可以停止。