function I = prostokaty(ff, a, b, N)
% wyznacza calke z ff(x) po [a,b] zlozona kwadratura
% prostokatow z podzialem na N przedzialow rownej dlugosci
% UWAGA: ff musi byc zwektoryzowana!

% wyznaczam wezly podzialu
x = linspace(a,b,N+1);

% wyznaczam srodki, czyli wezly kwadratury
% for i=1:N
%	v(i) = 0.5*(x(i)+x(i+1));
% end
v = 0.5*(x(1:N)+x(2:N+1));

h = v(2)-v(1); %krok podzialu

y = feval(ff, v); % wartosci w wezlach kwadratury
I = sum(y)*h; % kwadratura!

end
