变步长梯形求积法计算定积分.doc

《变步长梯形求积法计算定积分.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、变步长梯形求积法计算定积分1．原理：变步长求积法的思想是利用若干小梯形的面积代替原方程的积分，当精度达不到要求时，可以通过增加点数对已有的区间再次划分，达到所需精度时即可；其中由于新的式子中有原来n点中的部分项主要公式：T2n=Tn/2+(h/2)*Σf(xk+0.5);2．源程序如下：#include"math.h"#include"iostream.h"doublef(doublex){doubles;s=log(x*x);return(s);}doubleffts(doublea,doubleb,doubleeps){intn,k;doublefa,fb,h,t

2、1,p,s,x,t;fa=f(a);fb=f(b);n=1;h=b-a;t1=h*(fa+fb)/2;p=eps+1;while(p>=eps){s=0;for(k=0;k<=n-1;k++){x=a+(k+0.5)*h;s=s+f(x);}t=t1/2+h*s/2;p=fabs(t1-t);cout<<"步长n为:"<

3、<<"需要求解的积分式为f(x)=log(x^2)"<>a;cout<<"输入边界值b="<<'t';cin>>b;cout<<"输入误差限"<<'t';cin>>eps;result=ffts(a,b,eps);cout<<"经过变步长梯形求积法得方程结果为:"<