drawimage.m

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

%作出beta的图形
function ave = drawimage(P)
%调用格式: ave = drawimage(P)
%输入:     P - 50个算出的beta值
%输出:     ave - beta的均值
%图形输出: beta的振荡图和均值图


n=length(P)-1;
t([1:n+1])=[0:n];
ave=0;
for  i=1:n+1
    ave=ave+P(i);
end
ave=ave/(n+1);

plot(t,P,'b',[0,n],[ave ave],'r','linewidth',2);
axis([0,n,0,0.5]);
xlabel('k');
ylabel('\beta');
legend('Experimental value','Average value');