组合梯形公式求解积分.m

拉格朗日插值多项式拟合,牛顿插值多项式,欧拉方程解偏微分方程,使用极限微分求解导数（微分）,微分方程组的N=4龙格库塔解法,雅可比爹迭代法解方程AX=B,最小二乘多项式拟合,组合辛普生公式求解积分,用

%_______________________________________________________________________
%组合梯形公式求解积分
%说明：积分区间(a,b) ,被积函数f(x)，结果为h/2(f(a)+f(b))+h∑f(xk),k=1到M-1,k=0,1,2,...M
%xk=a+kh,注意x0=a,xM=b
%________________________________________________________________________
function [s]=trapr1(f,a,b,M)
%M为采样点个数
h=(b-a)/M;
s=0;
for k=1:(M-1)
    x=a+h*k;
    s=s+feval(f,x);
end
s=h*(feval(f,a)+feval(f,b))/2+h*s;