function int_kvad_(n)
%int_kvad_(n)
%
% Funkce vypocte obdelnikovym pravidlem aproximaci 
% integralu funkce x^2 na <0,1>.
%
% Vstupni parametr: n ... pocet dilku deleni intervalu <0,1>
% 29 dilku zaruci presnost 1e-4

h = 1 / n;     % delka dilku deleni
x = 0 : h : 1; % krajni body dilku

c = 0.5 * (x(1:n) + x(2:n+1));  % stredy jednotlivych dilku

f = c.^2;              % funkcni hodnoty ve stredech dilku
S = h * sum(f);        % priblizna hodnota integralu 
err = abs(1/3 - S);    % absolutni chyba vypoctu

disp(' ')
disp(['Obdelnikove pravidlo, n = ', num2str(n)])
disp('--------------------------------------------------')
disp(['Spoctena hodnota integralu: ', num2str(S,'%1.12g')])
disp(['Absolutni chyba: ', num2str(err,'%1.6e')])

% konec


