clc;
close all;
clear all;


for z=2
    for f=2
        for k=2
            for p=1
                
filename = 'piros_8cm_maxI.txt';
PEO = dlmread(filename);  

for i=1:6
    for j=1:6;
    PEO(i,j)=PEO(i,j)*((19.89-(2.73*(p-1)))/19.2);
    end;
end;

for i=1:6
    for j=1:6;
    P(i,j)=0;
    end;
end;
       
dvektor=[1 3 5];
rvektor=[1 3 5];

for d=dvektor
   for r=rvektor
      for i=1:6
           for j=1:6
                P(i,j)=P(i,j)+PEO(i+(6-r),j+(6-d));
           end;
       end;
   end;
end;

filename = 'kek_8cm_maxI.txt';
KEO = dlmread(filename);  

for i=1:6
    for j=1:6;
    KEO(i,j)=KEO(i,j)*((6.91-(0.8641*(k-1)))/6.9);
    end;
end;

for i=1:6
    for j=1:6;
    K(i,j)=0;
    end;
end;
       
dvektor=[2 4 6];
rvektor=[1 3 5];

for d=dvektor
   for r=rvektor
      for i=1:6
           for j=1:6
                K(i,j)=K(i,j)+KEO(i+(6-r),j+(6-d));
           end;
       end;
   end;
end;

filename = 'feher_8cm_maxI.txt';
FEO = dlmread(filename);  

for i=1:6
    for j=1:6;
    FEO(i,j)=FEO(i,j)*((13.9-(1.7959*(f-1)))/13.5);
    end;
end;

for i=1:6
    for j=1:6;
    F(i,j)=0;
    end;
end;
       
dvektor=[2 4 6];
rvektor=[2 4 6];

for d=dvektor
   for r=rvektor
      for i=1:6
           for j=1:6
                F(i,j)=F(i,j)+FEO(i+(6-r),j+(6-d));
           end;
       end;
   end;
end;

filename = 'zold_8cm_maxI.txt';
ZEO = dlmread(filename);  

for i=1:6
    for j=1:6;
    ZEO(i,j)=ZEO(i,j)*((5.56-(0.65626*(z-1)))/5.1);
    end;
end;

for i=1:6
    for j=1:6;
    Z(i,j)=0;
    end;
end;
       
dvektor=[1 3 5];
rvektor=[2 4 6];

for d=dvektor
   for r=rvektor
      for i=1:6
           for j=1:6
                Z(i,j)=Z(i,j)+ZEO(i+(6-r),j+(6-d));
           end;
       end;
   end;
end;



for i=1:6
    for j=1:6
        ossz(i,j)=Z(i,j)+F(i,j)+K(i,j)+P(i,j);
    end;
end;

Emax=max(ossz(:));
Emin=min(ossz(:));

elteres(1000*z+100*f+10*k+p) = ( Emax-Emin)/(Emax+Emin);



            end;
        end;
    end;
end;