trapezoidal.m

来自「辛普森方法」· M 代码 · 共 4 行

function INTF = trapezoidal(f,a,b)
%integral of f(x) over [a,b] by trapezoidal rule

INTF = (feval(f,a) + feval(f,b))*(b-a)/2;