function nintegral(m)
% m = pocet kroku, f = funkce (1 ... exp(x), 2 ... sin^3(x))

%a=0; b=2*pi;
a=pi/6; b= pi/2;
function f=f(x)
	f = exp(x);
    %f=sin(x).^3;
end

r=(b-a)/m; % velikost kroku

% obdelnikova metoa
S=0;
for i=a:r:b-r
	S=S+r*(f(i+r/2));
end
sprintf('Obdelnikova metoda:  %16.12g', S)
%disp(['obdelnikova metoda :' S]);
%S

% lichobeznikova metoda
S=0;
for i=a:r:b-r
	S=S+r*(f(i)+f(i+r))/2;
end
sprintf('Lichobeznikova metoda:  %16.12g', S)

% simpsonovo pravidlo
S=0;
for i=a:2*r:b-2*r
	S=S+r/3*(f(i)+4*f(i+r)+f(i+2*r));
end
sprintf('Simpsonovo pravidlo:  %16.12g', S)


end