for (s = 0; s < Z; s++)
    {
          for(m=0 ; m < Z ; m++)
        {
            valor=0.0;
            for(i=0; i < X; i++)
        for(j=0; j < Y; j++)
            {
R2[m+Z*s]+=HIM[i*X+j+(s*X*Y)]*HIM[i*X+j+(m*X*Y)]/(X*Y*Z);
            }
        }
    }