clear
f=[4.3 4.1 4.0 4.2 4.5 4.9 4.4 3.8 3.6 3.2
    4.5 4.3 4.2 4.4 4.6 4.8 4.6 4.2 4.0 3.9
    4.6 4.3 4.2 4.2 4.4 4.7 4.4 4.2 4.1 4.2
    4.5 4.1 4.1 4.0 3.9 3.8 3.6 3.4 3.7 3.9
    4.2 3.9 3.8 3.6 3.7 3.9 3.7 3.5 3.6 3.7
    4.0 3.8 3.7 3.7 3.6 3.5 3.2 3.4 3.6 3.4];

%Arvutame funktsiooni g väärtused
for i=1:6
    g(i)=0;
    for j=1:9
        g(i)=g(i)+f(i,j)+f(i,j+1);
    end
    g(i)=0.1*g(i); %0.2/2, 0.2 on samm
end

%Arvutame kubatuurvalemi vastuse
integraal=0;
for i=1:5
    integraal=integraal+g(i)+g(i+1);
end
integraal=0.1*integraal %0.2/2, kus 0.2 on samm





%suurus g_i
for i=1:6
    g(i)=0;
    for j=1:9
        g(i)=g(i)+f(i,j)+f(i,j+1);
    end
    g(i)=0.2/2*g(i); %0.2 on sammupikkus
end

%teise muutuja järgi integraal
S=0;
for i=1:5
    S=S+g(i)+g(i+1);
end
S=0.2/2*S; %0,2 on teine sammupikkus
vastus=S