Integral

h=0.5; %Lépésköz
x=(-1:h:2);
n=length(x);
fx=x.^3-2*x.^2+5; %Függvény
y=zeros(n-1,2); %Az y értékek halmaza és a növekedést

x2=(n-1);
for i=1:n-1
    x2(i)=(x(i)+x(i+1))/2;
end

for i=1:n-1
    y(i,1) = fx(i);
    y(i,2) = fx(i+1)-fx(i);
end

%%Ábrázolás
figure(1); %Az ablak
b1 = bar(x2,y(1:end,1),'BarWidth',1); %Az x2 (két x érték közti pontok) és az y értékek szerint ábrázolja)
hold('on');
b2 = bar(x2,y(1:end,1)+y(1:end,2),'BarWidth',1);
alpha(0.5);
plot(x,fx,'LineWidth',2); %Az y értékeit ábrázolja x függvényében
hold('off');
xlabel('x'); %X tengely elnevezés
ylabel('y'); %Y tengely elnevezés

%Területszámítás
disp('Az elso képlet szerinti terület');
area1=h*sum(fx(1:end-1));
area2=h*sum(fx(2:end));
fx2=x2.^3-2*x2.^2+5;
area3=h*sum(fx2(1:end));