复化梯形公求积分.doc

《复化梯形公求积分.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、计算机科学与工程学院武汉工程大学计算机科学与工程学院《计算方法》实验报告实验项目复化梯形公式求积分实验类别综合实验实验目的及要求实验目的：学会用复化梯形公式求函数的积分，并应用该算法于实际问题。实验要求：要求能随意输入被积函数，进行算法设计，打印出误差限例题：求被积函数在0

2、完整、体现收获40分说明：评阅教师：6计算机科学与工程学院日期：年月日实验内容设计分析复化数值积分:将区间[a,b]n等分，取等距节点由定积分的区间可加性，有在每一个小区间上利用梯形积分公式有一般记称做n+1点复化梯形积分公式。数学公式：算法描述：Step1:输入a,b和正整数n；Step2:置h=(b-a)/n;Step3:F=f(a)+f(b);l=0;Step4:对j=1,2,…,n循环执行5；Step5:置x=a+jh;l+=f(x);Step6:置T=h(F+2l)/2Step7:输出T;程序

3、源代码:#include#includeusingnamespacestd;6计算机科学与工程学院doublef(doublex)//求函数的值；{return1/(1+pow(sin(x),2.0));}voidEchelonIntegral(intn)//梯形积分{doubley=0;doubleh=(1.0-0.0)/n;for(inti=0;i<=n;i++){doublea=0.0+i*h,b=0.0+(i+1)*h;y=y+h*(f(a)+f(b))/2.

4、0;}cout<<"对应所求的梯形求积分为"<

5、}/*voidParabolicIntegral(intn)//抛物线积分{doubley=0;doublet=(1.0-0.0)/n;for(inti=0;i

6、pow(sin(x),2))的值****************"<>m1;EchelonIntegral(m1);cout<>answer1;while(answer1=='y'){cout<<"请输入把0到1的范围几等分？"<<"t";cin>>m1;EchelonIntegral(m1);//3.直线求积

7、分；cout<>answer1;}cout<

8、，在此次实验中，也有一些不近人意的地方，也需要我不断改进和练习。比如误差怎样取值，保证其最小。6