% 2000/3
clear all;
close all;
clc;

dx = 1;

x = -10:dx:10;

y = x.*x;

%Integrace
iy = 0;

for I=1:size(y,2)
    So = y(I)*dx;
    if I>1
        iy(I) = iy(I-1) + So;
    end
end

subplot(1,3,1);
plot(x,y, 'b')
xlabel('x');
ylabel('y');
title('Interovaná funkce');
grid on;
hold on;
bar(x,y, dx);
subplot(1,3,2);
plot(x, iy, 'g')
legend('Vypočítaná interovaná funkce');

iy2 = x.*x.*x/3;
subplot(1,3,3);
scatter(x, iy, 'r')