clc;
%krok (h):
h=0.01;
%granice całkowania (a,b):
a=2;
b=8;
%dziedzina funkcji:
x=a:h:b;
%funkcja f(x) = x^2-2x+sin(x) :
y=x.^2-2*x+sin(x);
%ilość kroków (n):
n=length(y);
%początkowe pole = 0 (s):
s=0;
%iteracja do wyczerpania kroków:
for i=1:n-1
    d=h*((y(i)+y(i+1))/2);
    s=d+s;
end;
%wynik:
s
plot(x,y)