计算方法实验六-----数值积分(共8页).doc
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_05.gif)
《计算方法实验六-----数值积分(共8页).doc》由会员分享,可在线阅读,更多相关《计算方法实验六-----数值积分(共8页).doc(8页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、精选优质文档-倾情为你奉上山西大学计算机与信息技术学院实验报告姓 名学 号专业班级课程名称计算方法实验日期成 绩指导老师批改日期实验名称实验六 数值积分一 实验目的: 利用复化梯形公式、复化辛普生公式和龙贝格数值积分公式计算的近似值。二 实验方法:(1) 将a,b区间n等分,记分点为,并在每个小区间上应用梯形公式 (2) 在每个小区间上,用辛普生公式 式中为的中点,即(3) 先用梯形公式计算,然后,将求积区间(a,b)逐次折半的方法,令区间长度 计算,式中。 于是,得到辛普生公式。 柯斯特求积公式。 最后,得龙贝格求积公式。 利用上述各公式计算,直到相邻两次的积分结果之差满足精度要求。三 实验
2、内容利用复化梯形公式、复化辛普生公式和龙贝格数值积分公式计算的近似值,要求误差为,将计算结果与精确值比较,并对计算结果进行分析(计算量、误差) 四 实验程序:复合梯形公式:#include #include #define esp 0.5e-7#define a 1 #define b 2 #define c 0#define d 1#define E 2.#define f1(x) (x*pow(E,x) #define f2(x) (4/(1+(x*x) void fun1()int i,n,k=0; double h,q,t,g; n=1; h=(double)(b-a)/2; t=h*
3、(f1(a)+f1(b); do k+; q=t; g=0; for (i=1;iesp); printf(函数1分了%d次:n,k);printf(积分结果为:); printf(%12.8lfn,t); void fun2()int i,n,k=0; double h,q,t,g; n=1; h=(double)(d-c)/2; t=h*(f2(c)+f2(d); do k+; q=t; g=0; for (i=1;iesp);printf(函数2分了%d次:n,k);printf(积分结果为:); printf(%12.8lfn,t); int main()printf(/*复合梯形公式
4、*/n); fun1();fun2();return 0;复合辛普生公式代码:#include #include #define esp 0.5e-7#define a 1 #define b 2 #define c 0#define d 1#define E 2.#define f1(x) (x*pow(E,x) #define f2(x) (4/(1+(x*x) void fun1() int i,n,k=0; double f1,f2,f3,h,s0,s; f1=f1(a)+f1(b); f2=f1(double)(b+a)/2); f3=0; s=(double)(b-a)/6)*(f
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算方法 实验 数值 积分
![提示](https://www.deliwenku.com/images/bang_tan.gif)
限制150内