trapezrule.m

来自「用matlab写的一些数值算法」· M 代码 · 共 11 行

function  T = trapezrule(f,a,b,m)
% Approximate integral by trapezoidal rule

% Version 4.06.2004.  INCBOX

x = linspace(a,b,m+1);                   % grid points
T = (feval(f,a) + feval(f,b))/2;    % endpoint contrib.
for  i = 1 : m-1
  T = T + feval(f,x(i+1));    % interior point contrib.
end
T = (b-a)/m * T;                      % multiply by  h