How do you plot a Midpoint Riemann Sum

function[] = graph_integration_method(fhandle,lo,hi,N)

hold on;

fplot(fhandle,[lo, hi]);

x = linspace(lo,hi,N+1);

dx = diff(x);

for index = 1:N

rectangle('position',[x(index) 0 dx(index) fhandle(index)/2]);

end

end