initiate.m

用Monte-Carlo方法解决二维Ising模型

%随机选取初始状态
function sigma=initiate(n)
%调用格式: sigma = initiate(n)
%输入:     n - 二维Ising模型的边长-1
%输出:     初始状态sigma

sigma=zeros(n,n);
for i=1:n
    for j=1:n
     %   sigma(i,j)=unidrnd(2)*2-3;%sigma的每一格点的值随机为1或-1.
       sigma(i,j)=1;
    end
end

end