复化梯形公式求积分.doc

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

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

2、月日实验内容设计分析复化数值积分:将区间[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;程序源代码:#include#includeusingna

3、mespacestd;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.0;}cout<<"对应所求的梯形求积分为"<

4、;doublex=mid,p=2*sin(2*x)*sin(2*x)*(1+sin(x)*sin(x))-2*cos(2*x)*pow(1+sin(x)*sin(x),2);doublef2d=p/pow(1+sin(x)*sin(x),2);for(i=0;i

5、>m1;EchelonIntegral(m1);

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

7、小区间上采取低阶数值积分公式，然后把这些小区间上的数值积分结果加起来作为函数在整个区间上的近似，类似于分段差值。在此次实验中，我学会了怎样以专业的知识，从专业的角度来考虑所遇见的问题或者与专业相关的问题，并以自己的看法，给予实现。同时，在此次实验中，也有一些不近人意的地方，也需要我不断改进和练习。比如误差怎样取值，保证其最小。