function mp = midpt(f,rng,N)

a = rng(1); b = rng(2);
h = (b-a)/N;
x = a + h*((1:N)-0.5);
mp = h*sum(f(x));