复化梯形公式matlab程序_复合梯形公式（compound trapezoidal formula）

复合梯形公式(compound trapezoidal formula)是1993年公布的数学名词。

将积分区间分为若干份, 在每一个“小区间”上用低阶梯形求积公式可得

复合梯形公式的收敛阶为2阶。

matlab程序

function I = ftrapz(fun,a,b,n)%fun,a,b,n分别为被积分函数、积分下限、积分上限、积分区间数目h = (b-a)/n;%求积分区间长度x = linspace(a,b,n+1);%bailinspace(x1,x2,N)用于产生x1,x2之间的N点行线性的矢量；%其中x1、x2、N分别为起始值、终止值、元素个数。%其中x1、x2、N分别为起始值、终止值、元素个数。y = feval(fun,x);I = h * (0.5*y(1) + sum(y(2:n)) + 0.5*y(n+1) );

被积函数f(x) Matlab表示

f=@(x)(sin(1./x)); %@表示参数；%(x)表示一个参数；%(sin(1./x))表示被积分函数；

求解

f=@(x)(sin(1./x)); m=ftrapz(f,1,1.5,5);format