计算方法上机题doc
计算方法习题
《计算方法》练习题一练习题第1套参考答案一、填空题1.二=3.14159…的近似值3.1428，准确数位是( 10，)。
f "心2.满足f(a)二c, f(b)二d 的插值余项R(x)二( (x-a)(x-b) )。
23.设｛R(x)｝为勒让德多项式，则(P2(X),F2(X))二(- )。
54•乘幕法是求实方阵(按模最大)特征值与特征向量的迭代法。
5•欧拉法的绝对稳定实区间是( [-2,0])。
二、单选题1 .已知近似数a,b,的误差限；(a), ；(b)U ；(ab) = (C )。
A. E(a)g(b)B. g(a)+w(b)C. a s(a) +|^(b)D. a^(b) + b 客(a) 2•设f (x) =X2 x，则f [1,2,3] = ( A )oA.lB.2C.3D.43 13 .设A= | ，则化A为对角阵的平面旋转 e = ( C ).〔1 3 一兀A.—2兀B. 一3JIC. D.4兀64 .若双点弦法收敛，则双点弦法具有(B )敛速.A.线性B.超线性C.平方D.三次5 .改进欧拉法的局部截断误差阶是( C).2 3 4A. o(h)B. o(h )C. o(h )D. o(h )三、计算题| x1 X2 — 31.求矛盾方程组：x1 2x^ 4的最小二乘解。
X r _ X2 = 2「(X1,X2)=(X1 X2 -3)2 (x「2X2 -4)2(X1 -X2 -2)2,由一=0： =0 得: -X] :x23X1 *2X? = 9 2x1 +6x2 =9解得X1 =^,X2 =2 o7 14|2x 1 5x 2 3x 3 = 6 I3•用列主元消元法解方程组： <2X 1 +4X 2 +3X 3 =5。
4禺 +6x 2+2x 3= 42 53 64 62 44 6 2 4 2 4 35 T1 2 3 T 2 2 4 4 6 2 4i 2 2 4i i1 1回代得：^(-1,1,1)T计算得：x 1 =0.25。
计算方法复习题
计算方法复习题一、判断题1．四舍五入得到的最后一位数字是有效数字。
（ ）2．运算量是衡量一个算法好坏的唯一指标。
（ ）3．从计算方法近似解角度考虑，方程组都有解。
（ ）4．最小二乘拟合本质是解矛盾方程组。
（ ）5．高斯—塞德尔迭代法一定比雅可比迭代法收敛速度快。
（ ）6．数值积分中求积系数与被积函数f (x )有关。
（ ）7．同一组数据采用拉格朗日插值与牛顿插值的结果不同。
（ ）8．迭代法求非线性方程f (x )=0收敛的条件是|f ’(x)|<1。
（ ）9．常微分方程数值解中龙格库塔法的系数可由Taylor 公式展开求取。
（ ）10．线性方程组的迭代法不适合用于求解大型稀疏矩阵。
（ ）11．加减计算量是衡量一个算法好坏的最重要的指标。
（ ）12．计算方法应考虑各种误差的影响。
（ ）13．插值法是函数逼近的唯一方法。
（ ）14．求解同一个问题时，结果的有效数字位数越多说明的近似解精度越高。
（ ）15．高斯—塞德尔迭代法不一定比雅可比迭代法求解精度高。
（ ）16．所有插值法都是只要求构造的φ(x )与f (x )在给定点的函数值相等。
（ ）17．f(x)没有解析表达式，只有数表形式时，可以对f (x )进行积分。
（ ）18．线性方程组的直接解法适合用于求解小型稠密矩阵。
（ ）19．可以用代数精确度度量数值积分的精度。
（ ）20．计算方法中各种算法只考虑舍入误差。
（ ）21．计算方法考虑数学问题的近似解，信息量越少近似解越准确。
（ ）22．所有插值法只要求构造的φ(x )与f (x )在给定点的函数值相等。
（ ）23．线性方程组迭代收敛与矩阵A 的特征值有关。
（ ）24．可以用代数精确度度量数值积分的精度。
（ ）二、填空题1．微分方程离散化的方法有：数值积分、差商和_________________。
2．你学习或知道的线性方程组求解方法，除了简单迭代法（雅克比）外，还有____________等。
(完整版)统计计算方法复习考试题
二、选择题：
6、能产生等可能取值为 中一个数的MATLAB程序是（ ）
(A)ceil(5*rand)(B)ceil(4*rand)(C)floor(4*rand)(D)randperm(4)
7、在MATLAB中,表示负二项分布的概率密度函数的是（ ）
(A) binopdf (B)binocdf (C)nbinpdf (D)nbincdf
一、填空题：
1、若随机变量 的概率密度为 ,则 的方差为。
2、若 服从二项分布B(5000,0.001),则由泊松定理知 。
3、若 服从均值为5的指数分布，则 。
4、设 服从参数为2的泊松过程,则 。
5、设 的概率密度为 ,则其分布函数的逆函数为。
二、选择题：
6、能产生等可能取值为 中一个数的MATLAB程序是（ ）
（2） 某一时段的状态为0，定义为初始状态，即 ，所求概率为：
17、首先由C-K方程得两步转移矩阵为：
一、填空题：
1、若随机变量 的概率密度为 ,则 的方差为。
2、若 服从二项分布B(500,0.01),则由泊松定理知 。
3、若 服从失效率为0.05的指数分布，则 。
4、设 服从参数为0.5的泊松过程,则 。
MATLAB程序：
alpha=5;beta=3; U=rand;
X=(-log(U)/alpha)^(1/beta
因 ,计算得
，
由P值为0.8931，说明不能拒绝原假设，即不认为发生事故与星期几有关。
16、(1)一步转移概率可用频率近似地表示为：
所以一步转移矩阵为： ；
(A)ceil(5*rand)(B)floor(5*rand)(C)floor(6*rand)(D)randperm(5)
计算机二级上机考试(程序改错题)1
程序编程题目1、第一大类(交换类)共7道2. 程序修改（第3套）给定程序中函数fun 的功能是：通过某种方式实现两个变量的交换，规定不允许增加语句和表达式。
例如变量a中的值原为8，b中的值原为3，程序运行后a中值为3，b 中的值为8./**found**/t=*x; *x=y;/**found**/return(t);2. 程序修改（第60套）给定程序中函数fun的功能是：实现两个整数的交换。
例如给a和b分别输入60和65，输出为：a=65 b=60/**found**/void fun(int *a,int *b)/**found**/t=*b;*b=*a;*a=t;2. 程序修改（第88套）给定程序中函数fun的功能是：将主函数中两个变量的值进行交换。
例如，若变量a中的值为8，b中的值为3，则程序运行后，a中的值为3，b中的值为8。
/**found**/void fun(int *x, int *y)/**found**/t=*x,*x=*y,*y=t;2、第二大类(计算类)共22道2. 程序修改（第1套）给定程序中函数fun的功能是：求出以下分数序列的前n项之和。
2/1+3/2+5/3+8/5+13/8+21/13+……和值通过函数值返回main函数。
例如n=5,则应输出：。
/**found**/double fun(int n)/**found**/s=s+(double)a/b;2. 程序修改（第6套）给定程序中函数fun的功能是：用递归算法计算列中第n项的值。
从第一项起，斐波拉契数列为：1、1、2、3、5、8、13、21、……/**found**/switch(g)/**found**/case 1:case 2:return 1;2. 程序修改（第13套）给定程序中函数fun的功能是：求 s=aa…aa—…—aaa—aa—a（此处aa…aa表示n个a，a和n的值在1至9之间）。
例如，a=3,n=6,则以下表达式为：S=333333-33333-3333-333-33-3其值是296298。
计算方法复习题-试题卷
一计算题
1. 能不能用迭代法求解以下方程，如果不能时，试将方程改写成能用迭代法求解的形式。
2. 用矩阵的LU分解算法求解线性方程组
X1+2X2+3X3 = 0
2X1+2X2+8X3 = -4
-3X1-10X2-2X3 = -11
3. 用高斯消去法求解线性方程组
解：消元过程
4. 给定常微分初值问题试构造一个求解常微分初值问题的两步差分格式。
5. 用矩阵的Doolittle分解算法求解线性方程组
2X1+X2+X3 = 4
6X1+4X2+5X3 =15
4X1+3X2+6X3 = 13
6. 利用Doolittle分解法解方程组Ax=b，即解方程组
解：用公式
7. 用矩阵的Doolittle分解算法求解线性方程组 X1+2X2+3X3 = 1
2X1– X2+9X3 = 0
-3X1+ 4X2+9X3 = 1
解：
8. 用Doolittle分解法解方程组
解：方程组的系数矩阵为
根据分解公式得
9. 方程将其改写为
10. 用高斯消元法解方程组
解：方程组的扩大矩阵为
11. 方程将其改写为
解：注意到迭代公式的形式，
12. 用Doolittle三角分解法求解线性代数方程组：
解：由公式
13. 用高斯消去法求解线性方程组
2X1- X2+3X3 = 2
4X1+2X2+5X3 = 4
-3X1+4X2-3X3 = -3
解：方程组的扩大矩阵为
14. 给定方程
〔1〕分析该方程存在几个根；
〔2〕构造迭代公式，说明迭代公式是收敛的。
15. 用Euler方法求解
(取h=0.2)。
计算方法第二版课后练习题含答案
计算方法第二版课后练习题含答案前言本文将为大家提供计算方法第二版课后练习题的答案，旨在帮助读者更好地学习和掌握计算方法的知识。
本文全部内容均为作者整理，尽可能保证每一题的答案正确性。
读者可以借助本文的答案，检验自己的练习成果，加强对计算方法知识的理解和掌握程度。
同时，读者也应该注意切勿直接复制答案，本文的答案仅供参考，希望读者能够通过自己的思考和探索，获得更深层次的学习感悟。
第一章引论1.1 计算方法的基本概念和思想练习题 1写出计算方法的三要素，并分别简要解释。
答案计算方法的三要素为：模型、算法、误差分析。
•模型：计算方法所涉及的实际问题所对应的数学模型，是解决问题的基础；•算法：根据模型，构造相应的计算程序，即算法；•误差分析：计算结果与实际应用中所需的精度之间的差异，称为误差。
误差分析是对计算结果质量的保障。
1.2 算法的误差练习题 2写出二分法算法，并解释其误差。
答案算法：function binarySearch(a, target) {let low = 0;let high = a.length - 1;while (low <= high) {let midIndex = Math.floor((low + high) / 2);let midValue = a[midIndex];if (midValue === target) {return midIndex;} else if (midValue < target) {low = midIndex + 1;} else {high = midIndex - 1;}}return -1;}误差:二分法算法的误差上界为O(2−k)，其中k为迭代次数。
在二分法被成功应用时，k取决于与目标值x的距离，即 $k=\\log _{2}(\\frac{b-a}{\\epsilon})$，其中[a,b]是区间，$\\epsilon$ 是目标值的精度。
全国计算机等级考试一级上机Excel题库
全国计算机等级考试一级上机Excel题库第1题、请在“考试项目”菜单下选择“电子表格软件使用”菜单项，然后按照题目要求打开相应的子菜单，完成下面的内容，具体要求如下：注意：下面出现的所有文件都必须保存在考生文件夹[%USER%]下所有中英文状态的括号、小数位数必须与题面相符合。
，(1)打开工作簿文件table13.xls，将下列已知数据建立一抗洪救灾捐献统计表(存放在A1：D5的区域内)，将当前工作表Sheet1更名为“救灾统计表”。
，单位捐款(万元) 实物(件) 折合人民币(万元)第一部门1.95 89 2.45第二部门1.2 87 1.67第三部门0.95 52 1.30总计，(2)计算各项捐献的总计,分别填入“总计”行的各相应列中。
(结果的数字格式为常规样式) ，(3)选“单位”和“折合人民币”两列数据(不包含总计)，绘制部门捐款的三维饼图，要求有图例并显示各部门捐款总数的百分比，图表标题为“各部门捐款总数百分比图”。
嵌入在数据表格下方(存放在A8：E18 的区域内)。
，第2题、请在“考试项目”菜单下选择“电子表格软件使用”菜单项，然后按照题目要求打开相应的子菜单，完成下面的内容，具体要求如下：注意：下面出现的所有文件都必须保存在考生文件夹[%USER%]下所有中英文状态的括号、小数位数必须与题面相符合。
，（1）打开工作簿文件table.xls，请将下列两种类型的股票价格随时间变化的数据建成一个数据表存放在(A1:E7的区域内）,其数据表保存在sheet1工作表中。
股票种类时间盘高盘低收盘价A 10：30 114.2 113.2 113.5A 12：20 215.2 210.3 212.1A 14：30 116.5 112.2 112.3B 12：20 120.5 119.2 119.5B 14：30 222.0 221.0 221.5B 16：40 125.5 125.0 125.0 ，（2）对建立的数据表选择“盘高”、“盘低”、“收盘价”、“时间”数据建立盘高-盘低-收盘价簇状柱形图图表，图表标题为“股票价格走势图”，并将其嵌入到工作表的A9：F19区域中。
东南大学《数值分析》上机题
数值分析上机题1设2211NN j S j ==-∑，其精确值为1311221N N ⎛⎫-- ⎪+⎝⎭。
（1）编制按从大到小的顺序22211121311N S N =+++---，计算N S 的通用程序。
（2）编制按从小到大的顺序2221111(1)121N S N N =+++----，计算NS 的通用程序。
（3）按两种顺序分别计算210S ，410S ，610S ，并指出有效位数。
（编制程序时用单精度） （4）通过本上机题，你明白了什么？程序代码（matlab 编程）：clc cleara=single(1./([2:10^7].^2-1)); S1(1)=single(0); S1(2)=1/(2^2-1); for N=3:10^2 S1(N)=a(1); for i=2:N-1S1(N)=S1(N)+a(i); end endS2(1)=single(0); S2(2)=1/(2^2-1); for N=3:10^2 S2(N)=a(N-1);for i=linspace(N-2,1,N-2) S2(N)=S2(N)+a(i); end endS1表示按从大到小的顺序的S N S2表示按从小到大的顺序的S N通过本上机题，看出按两种不同的顺序计算的结果是不相同的，按从大到小的顺序计算的值与精确值有较大的误差，而按从小到大的顺序计算的值与精确值吻合。
从大到小的顺序计算得到的结果的有效位数少。
计算机在进行数值计算时会出现“大数吃小数”的现象，导致计算结果的精度有所降低，我们在计算机中进行同号数的加法时，采用绝对值较小者先加的算法，其结果的相对误差较小。
数值分析上机题220．（上机题）Newton 迭代法（1）给定初值0x 及容许误差ε，编制Newton 法解方程()0f x =根的通用程序。
（2）给定方程3()/30f x x x =-=，易知其有三个根1x *=，20x *=，3x *=。
《数值计算方法》试题集及答案(1-6)-2..
《计算方法》期中复习试题一、填空题：1、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得⎰≈31_________)(dx x f ,用三点式求得≈')1(f 。
答案：2.367，0.252、1)3(,2)2(,1)1(==-=f f f ，则过这三点的二次插值多项式中2x 的系数为 ，拉格朗日插值多项式为 。
答案：-1，)2)(1(21)3)(1(2)3)(2(21)(2--------=x x x x x x x L3、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字；4、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( )；答案)(1)(1n n n n n x f x f x x x '---=+5、对1)(3++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 )；6、计算方法主要研究( 截断 )误差和( 舍入 )误差；7、用二分法求非线性方程 f (x )=0在区间(a ,b )内的根时，二分n 次后的误差限为( 12+-n a b )；8、已知f (1)＝2，f (2)＝3，f (4)＝5.9，则二次Newton 插值多项式中x 2系数为( 0.15 )； 11、 两点式高斯型求积公式⎰1d )(xx f ≈(⎰++-≈1)]3213()3213([21d )(f f x x f )，代数精度为( 5 )；12、 为了使计算32)1(6)1(41310---+-+=x x x y 的乘除法次数尽量地少，应将该表达式改写为11,))64(3(10-=-++=x t t t t y ，为了减少舍入误差，应将表达式19992001-改写为199920012+ 。
13、 用二分法求方程01)(3=-+=x x x f 在区间[0,1]内的根,进行一步后根的所在区间为 0.5，1 ,进行两步后根的所在区间为 0.5，0.75 。
计算方法练习题集和答案解析
练习题与答案练习题一练习题二练习题三练习题四练习题五练习题六练习题七练习题八练习题答案练习题一一、是非题1.*x=–12.0326作为x的近似值一定具有6位有效数字，且其误差限£41021-⨯。
( )2.对两个不同数的近似数，误差越小，有效数位越多。
()3.一个近似数的有效数位愈多，其相对误差限愈小。
()4.用212x-近似表示c o s x产生舍入误差。
() 5. 3.14和 3.142作为π的近似值有效数字位数相同。
() 二、填空题1.为了使计算()()2334912111yx x x=+-+---的乘除法次数尽量少，应将该表达式改写为；2.*x=–0.003457是x舍入得到的近似值，它有位有效数字，误差限为，相对误差限为；3.误差的来源是；4.截断误差为；5.设计算法应遵循的原则是。
三、选择题1．*x=–0.026900作为x的近似值，它的有效数字位数为( ) 。
(A) 7； (B) 3；(C) 不能确定 (D) 5.2．舍入误差是( )产生的误差。
(A) 只取有限位数 (B) 模型准确值与用数值方法求得的准确值(C) 观察与测量 (D) 数学模型准确值与实际值3．用 1+x近似表示e x所产生的误差是( )误差。
(A). 模型 (B). 观测 (C). 截断 (D). 舍入4．用s *=21g t 2表示自由落体运动距离与时间的关系式 (g 为重力加速度)，s t 是在时间t 的实际距离，则s t s *是（ ）误差。
(A). 舍入 (B). 观测 (C). 模型 (D). 截断 5．1.41300作为2的近似值，有( )位有效数字。
(A) 3； (B) 4； (C) 5； (D) 6。
四、计算题1． 3.142，3.141，227分别作为π的近似值，各有几位有效数字？2． 设计算球体积允许的相对误差限为1%，问测量球直径的相对误差限最大为多少？3． 利用等价变换使下列表达式的计算结果比较精确：(1)1||,11211<<+-++x x x x ， (2) 1||1112<<+⎰+x dt t x x (3) 1||,1<<-x e x ， (4) 1)1ln(2>>-+x x x4．真空中自由落体运动距离s 与时间t 的关系式是s =21g t 2，g 为重力加速度。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第八题：一，题目从函数表 x 0.0 0.1 0.195 0.3 0.401 0.5 f(x) 0.39894 0.39695 0.39142 0.38138 0.36812 0.35206 出发，用下列方法计算f(0.15),f(0.31),f(0.47)的近似值： （1）分段线性插值；（2）分段二次插值;(3)全区间上拉格朗日插值(要求一次性输入整张函数表，并利用计算机选择在插值计算中所需的结点）。
二，基本方法（1）先选取两个结点x i 1-与x i 使∈x [x i 1-,x i ],然后在区间[x i 1-,x i ]上作线性插值，即得f(x)≈xx x x y x x x x yL i ii iii ii x 11111)(------+--=。
（2）选取距点x 最近的三个结点xi 1-，xi，xi 1+进行二次插值，即取f(x)≈)]([)(111,12∑∏+-=+≠-=--=i i k i kj i j jkjkxx x x y L x 。
（3）n=5,))...()()...())...()()...(()()(111100(x x x x xx x x x x x x x x x x yly L n k k k k kkn k k nk kknk knx x --------==+-+-==∑∑三，计算结果 (1) x 0.15 0.31 0.47 f(x)0.394039 0.380067 0.356927 (2)x 0.15 0.31 0.47 f(x)0.3944820.3802250.361399(3)x 0.15 0.31 0.47f(x) 0.394473 0.380219 0.357222 四，结果分析编程基本按照课本套公式，方法从（1）至（3），计算结果越来越精确。
五，源程序（1）#include<stdio.h>int main(void){int i,j;double a[6],b[6],x,y;printf("Enter x:");for(i=0;i<6;i++)scanf("%lf",&a[i]);printf("Enter f(x):");for(i=0;i<6;i++)scanf("%lf",&b[i]);for(j=0;j<3;j++){printf("Enter x:");scanf("%lf",&x);for(i=0;x>a[i];i++);y=b[i-1]*(x-a[i])/(a[i-1]-a[i])+b[i]*(x-a[i-1])/(a[i]-a[i-1]); printf("x=%lf , ",x);printf("y=%lf\n",y);}return 0;}（2）#include<stdio.h>int main(void){int i,j;double a[6],b[6],x,y;printf("Enter x:");for(i=0;i<6;i++)scanf("%lf",&a[i]);printf("Enter f(x):");for(i=0;i<6;i++)scanf("%lf",&b[i]);for(j=0;j<3;j++){printf("Enter x:");scanf("%lf",&x);for(i=0;x>a[i];i++);if((a[i+1]-x)<(x-a[i-2])){y=b[i-1]*(x-a[i])/(a[i-1]-a[i])*(x-a[i+1])/(a[i-1]-a[i+1]);y=y+b[i]*(x-a[i-1])/(a[i]-a[i-1])*(x-a[i+1])/(a[i]-a[i+1]);y=y+b[i+1]*(x-a[i-1])/(a[i+1]-a[i-1])*(x-a[i])/(a[i+1]-a[i]); printf("x=%lf , ",x); printf("y=%lf\n",y);}else{y=b[i-2]*(x-a[i-1])/(a[i-2]-a[i-1])*(x-a[i])/(a[i-2]-a[i]);y=y+b[i-1]*(x-a[i-2])/(a[i-1]-a[i-2])*(x-a[i])/(a[i-1]-a[i]);y=y+b[i]*(x-a[i-2])/(a[i]-a[i-2])*(x-a[i-1])/(a[i]-a[i-1]); printf("x=%lf , ",x); printf("y=%lf\n",y);}}return 0;}（3）#include<stdio.h>int main(void){int i,j,k;double a[6],b[6],x,y,p;printf("Enter x:");for(i=0;i<6;i++)scanf("%lf",&a[i]);printf("Enter f(x):");for(i=0;i<6;i++)scanf("%lf",&b[i]);for(j=0;j<3;j++){y=0;printf("Enter x:");scanf("%lf",&x);for(i=0;i<6;i++){p=1;for(k=0;k<6;k++){if(k!=i)p=p*(x-a[k])/(a[i]-a[k]);}y=y+p*b[i];}printf("x=%lf , ",x); printf("y=%lf\n",y);}return 0;}六，附图（1）（2）（3）第九题：一，题目已知直升机旋转机翼外形曲线轮廓线上的某些型值点及端点处的一阶导数值y`(x0)=1.86548,y`(x18)=-0.046115,试计算该曲线上横坐标为2,4,6,12,16,30,60,110,180,280,400,515处点的纵坐标（要求该曲线具有二阶光滑度）。
k 0 1 2 3 4 5 6二，基本方法根据课本图4-9流程图及追赶法解方程。
三，计算结果 x 246121630y7.834331 10.476880 12.358485 16.579653 19.093173 25.386597x 60110180280 400 515y 32.804369 36.647886 35.917141 29.368427 16.799141 0.524459四，结果分析利用追赶法求M[i]，比高斯消去法步骤更简单，虽然结果有微小偏差，尤其在s[515]处与课本结果偏差较大。
但在精度要求不高的情况下，追赶法可以减少工作量，提高工作效率。
五，源程序#include<stdio.h>x k 0.52 3.1 8.0 17.95 28.65 39.62 50.65 y k5.28794 9.4 13.84 20.2 24.9 28.44 31.1k13 14 15 16 17 18 x k 364.4 416.3 468 494 507 520 y k20.914.87.83.71.50.2k7 8 9 10 11 12 x k 78 104.6 156.6 208.6 260.7 312.5 y k35 36.5 36.6 34.6 31.0 26.34int main(void){int i,j;doublea[19],b[19],h[19],s[18],t[18],g[19],x,y,p,q,f[19],m[20],n[20];printf("输入左端点一阶导:");scanf("%lf",&p);printf("输入右端点一阶导:");scanf("%lf",&q);printf("Enter x:");for(i=0;i<19;i++)scanf("%lf",&a[i]);printf("Enter f(x):");for(i=0;i<19;i++)scanf("%lf",&b[i]);for(i=0;i<19;i++)h[i]=a[i]-a[i-1];for(i=1;i<18;i++)s[i]=h[i]/(h[i]+h[i+1]);for(i=1;i<18;i++)t[i]=1-s[i];for(i=1;i<18;i++)g[i]=6/(h[i]+h[i+1])*((b[i+1]-b[i])/h[i+1]-(b[i]-b[i-1])/h[i]);g[0]=6/h[1]*((b[1]-b[0])/h[1]-p);g[18]=6/h[18]*(q-(b[18]-b[17])/h [18]);f[1]=0.5;for(i=2;i<19;i++)f[i]=t[i-1]/(2-s[i-1]*f[i-1]);n[1]=g[0]/2;for(i=2;i<20;i++)n[i]=(g[i-1]-s[i-1]*n[i-1])/(2-s[i-1]*f[i-1]);m[19]=n[19];for(i=18;i>0;i--)m[i]=n[i]-f[i]*m[i+1];printf("输出M:\n");for(i=1;i<20;i++)printf("%lf ",m[i]);for(j=0;j<12;j++){printf("\n");printf("Enter x:");scanf("%lf",&x);for(i=0;x>a[i];i++);y=m[i]*(a[i]-x)*(a[i]-x)*(a[i]-x)/h[i]/6;y=y+m[i+1]*(x-a[i-1])*(x-a[i-1])*(x-a[i-1])/h[i]/6;y=y+(b[i-1]-m[i]*h[i]*h[i]/6)*(a[i]-x)/h[i];y=y+(b[i]-m[i+1]*h[i]*h[i]/6)*(x-a[i-1])/h[i];printf("f(x):%lf",y);}return 0;}六，附图第十二题一，题目某类疾病发生率y‰和年龄段x(每五年为一段，例如0-5岁为第一段，6-10岁为第二段)之间有形如e bx=的关系。