run.m

This is a Matlab v6.5 implementation of the Romberg integration algorithm.

% quadrature

clc;
disp ('Calkowanie numeryczne:');

x = linspace (-5, 5, 100);
grid;
plot (x, f(x), 'r');
grid; pause; close;

a = input ('Podaj a: ');
b = input ('Podaj b: ');
n = input ('Podaj liczbe podprzedzialow: ');

if (a > b)
	disp ('a musi byc mniejsze od b!');
else
	% wewnatrz funckji przedzial <a, b> zostanie podzielony na n podprzedzialow
	s = rom (a, b, n);

	disp (['calka = ', num2str(s)]);
	rys = input ('Narysowac wykres funkcji? (t, n): ', 's');
	if (rys == 't')
		% tutaj rysuje wykres funkcji
		% przez 500 prostokotow zeby zamalowac pole pod funkcja
		x = linspace (a, b, 500);
		bar (x, f(x));
		grid; pause; close;
	end;
end;