function Riemannsumme(n)

x=0:0.002:2;
y=(x-0.5).^2+1;
plot([0, x],[0, y],'k');
pause(2);
for k=0:n,
    hold on;
    i=2^k;
    for j=0:i-1;
      patch([2*j/i, 2*j/i, 2*(j+1)/i, 2*(j+1)/i],[0,y(floor((j+1/2)*1000/i)+1),y(floor((j+1/2)*1000/i)+1),0],'b');
    end;%%%%%% Riemannsche Summe mit dem Mittelpunkt jedes Intervalls
    plot([0, x],[0, y],'k');
    hold off;
     pause(1);
    if k <n
    plot([0, x],[0, y],'k');
    end;
end;