trapezc.m

来自「various matlab programs to slove various」· M 代码 · 共 5 行

function int = trapezc(a,b,m,fun)
h=(b-a)/m; x=[a:h:b]; dim = max(size(x)); y=eval(fun); 
if size(y)==1, y=diag(ones(dim))*y; end; 
int=h*(0.5*y(1)+sum(y(2:m))+0.5*y(m+1));