% rom (a, b, n)		- calkowanie metoda Romberga (trapezow)
%			  <a, b> - przedzial calkowania
%			  n - dokladnosc (ilosc podprzedzialow)
%			  funkcja f musi byc zdefiniowana
function r = rom (a, b, n)

    if (nargin == 0)
	help rom;
    else
	h = (b - a)/n;

	s = 0;
	t = a + h;
	for i = 1:(n-1)
	    s = s + f(t);
	    t = t + h;
	end;
	s = 0.5 * h * (f(a) + f(b) + 2 * s);

	r = s;
    end;