trapezoidal.m

辛普森方法

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;