%Algoritmus prevzat z : Vysoka skola banska-
%Technicka univerzita Ostrava, Fakulta stavebni
% Prednaska z predmetu: Algoritmizace inzenyrskych vypoctu
%Tema c. 6: Numericka integrace urciteho integralu
%http://fast10.vsb.cz/krejsa/studium/algoritmy_06.pdf
%<cit. 31.1.2016>
function [integral]=romberg(f,a,b,n)
h=(b-a)./(2.^(0:n-1));
integral(1,1)=(b-a)*(f(a)+f(b))/2;
for j=2:n
s=0;
for i=1:2^(j-2)
s=s+f(a+(2*i-1)*h(j));
end
integral(j,1)=integral(j-1,1)/2+h(j)*s;
for k=2:j
integral(j,k)=(4^(k-1)*integral(j,k-1)-integral(j-1,k-1))/(4^(k-1)-1);
end
end
end
