function Unterober(n)

x=0:0.002:2;
y=(x-0.5).^2+1;
plot([0, x],[0, y],'b');
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,min(y(floor(j*1000/i)+1:floor((j+1)*1000/i)+1)),min(y(floor(j*1000/i)+1:floor((j+1)*1000/i)+1)),0],'r');        
      patch([2*j/i, 2*j/i, 2*(j+1)/i, 2*(j+1)/i],[min(y(floor(j*1000/i)+1:floor((j+1)*1000/i)+1)),max(y(floor(j*1000/i)+1:floor((j+1)*1000/i)+1)),max(y(floor(j*1000/i)+1:floor((j+1)*1000/i)+1)),min(y(floor(j*1000/i)+1:floor((j+1)*1000/i)+1))],'g');
    end;
    plot([0, x],[0, y],'b');
    hold off;
    pause(1);
    if k <n
    plot([0, x],[0, y],'b');
    end;
end;