%script general para hacer una corrida  a un set de parámetros,
%beta, tamaño de la red,

%Tamaño en x,
%Lx=10;
%Tamaño en y,
%Ly=10;
% beta = 1/T
%beta=0.1;

%propongo un estado inicial al azar

%Sij es una matriz de 1 y -1 indicando las dos proyecciones de
%espin

Sij=2*(rand(Lx,Ly)>0.5) -1;

%pcolor(Sij);
%npre = 100;
%npasos = 1000;
energia=zeros(npasos+1,1);
magnet=zeros(npasos+1,1);

%pretermalizo
for n=1:npre
[Sij,DE,DM]=ising2Dpaso(Sij,beta);
end

energia(1) = En(Sij);
magnet(1) = sum(sum(Sij));
for n=1:npasos
[Sij,DE,dM] = ising2Dpaso(Sij,beta);
energia(n+1) = energia(n) + DE;
magnet(n+1) = magnet(n) + dM;

%if(mod(n,10) == 0)
%pcolor(Sij);shading flat;
%title(['beta = ' num2str(beta) ' e = ' num2str(energia(n)/(Lx*Ly)) ' m = ' num2str(magnet(n)/(Lx*Ly))]);
%drawnow;
%end


end