data = dlmread('Resultado.txt');
valores = dlmread('Valores Iniciais.txt');

erro1 = zeros(10320, 84);
c = 0;
for i = 1:4
    for j = 1:21
        erro1(:,j+c) = abs(data(:,j+c) - valores(:, i));
    end
    c = c + 21;
end

dlmwrite('erro 1.txt', erro1, 'delimiter', ' ');

erros = zeros(10320, 84);
for i = 1:10320
    c = 0;
    for j = 1:4
        for k = 1:21
            if(i == 162)
                t = 1;
            end
            inic = valores(i,j);
            if(inic == 0)
                inic = 1;
            end
            erros(i,k+c) = (erro1(i,k+c)/inic)*100;
        end
        c = c + 21;
    end
end

dlmwrite('erro 2.txt', erros, 'delimiter', ' ');

statisticAnalysis = zeros(3, 84);
for i = 1:84
    e = erros(:, i);
    statisticAnalysis(1, i) = mean(e);
    statisticAnalysis(2, i) = std(e);
    statisticAnalysis(3, i) = max(e);
end

dlmwrite('Analise Estatística.txt', statisticAnalysis, 'delimiter', ' ');