midpntc.m

various matlab programs to slove various problems

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);