%估計M*M格點的臨界beta
function beta = solvebeta(M)
%調用格式: beta=solvebata(M)
%自變量:   M - 二維Ising模型的邊長
%輸出:    beta - 臨界溫度值


%請選擇合適的策略,否則程序運行非常緩慢!
t=20; %重復次數
p=50; %beta個數
N=500;N1=500;N2=3;number=1; %組策略



P=zeros(t,1);
for k=1:t
    [U,C]=Ising(M,p,N,N1,N2,number,2,0);
    P(k)=(maxN(C)-1)/p;
end
beta=drawimage(P);
end

function P=maxN(C)
n=length(C);
P=1;
for j=2:n
    if C(j)>C(P);
        P=j;
    end
end
end