clear
clc
%przedział:
a=-1;
b=4;

%metoda trapezów dla kroku 0.2 :
h=0.2;
x=a:h:b;
y=x.^2+2;
m=0;
n=(b-a)/h;
for k=1:n
    ysr=(y(k)+y(k+1))/2;
    m=m+ysr;
end;
p=h*m %podany wynik
plot(x,y)
grid on;

%metoda trapezów dla kroku 0.5 :
h=0.5;
x=a:h:b;
y=x.^2+2;
m=0;
n=(b-a)/h;
for k=1:n
    ysr=(y(k)+y(k+1))/2;
    m=m+ysr;
end;
p=h*m %podany wynik

%metoda trapezów dla kroku 1 :
h=1;
x=a:h:b;
y=x.^2+2;
m=0;
n=(b-a)/h;
for k=1:n
    ysr=(y(k)+y(k+1))/2;
    m=m+ysr;
end;
p=h*m %podany wynik