function En=En(Sij)

[Lx,Ly]=size(Sij);
En = 0;
for i=1:Lx
    for j=1:Ly
        
    %busco primeros vecinos, con condiciones periodicas.
    iright = mod(i,Lx) + 1; 
    ileft = mod(i+Lx-2,Lx) + 1;

    jup = mod(j,Ly) + 1;
    jdown = mod(j+Ly-2,Ly) + 1;

    En = En - Sij(i,j)*(Sij(i,jup) + Sij(i,jdown) + Sij(iright,j) + Sij(ileft,j));
    
    end
end

   %Como conté  el doble, divido.
   En = En/2;

end