function T = Romberg(a,b,f,n)
%ROMBERG Romberg-Verfahren

h(1)=b-a;
T(1,1)=h(1)*(f(a)/2+f(b)/2);
for i=2:n
    h(i)=h(i-1)/2;
    T(i,1)=1/2*T(i-1,1);
    for k=1:2^(i-2) % Index in der Folge ungerader Zahlen zwischen 1 und 2^(i-1)
        T(i,1)=T(i,1)+h(i)*f(a+(2*k-1)*h(i));
    end
    for j=2:min(i,n)
        T(i,j)=T(i,j-1)+(T(i,j-1)-T(i-1,j-1))/(4^j-1);
    end
end

end

