midpntc.m

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

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