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

来自「拉格朗日插值多项式拟合,牛顿插值多项式,欧拉方程解偏微分方程,使用极限微分求解导」· M 代码 · 共 15 行

%_______________________________________________________________________
%组合梯形公式求解积分
%说明：积分区间(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;