Ising模型的matlab数值模拟

clc
clear 
format long
H=300;
L=300;                          %模拟的晶格大小
k=1.3806488*10^（-23）;
Tc=3;        %临界温度
T=1;       %温度
s=round（rand（HL））*2-1;     %-1和1的矩阵       
colormap（[1 0 0;0 1 0]）;        %控制两种磁矩的颜色
h1=imagesc（s）;      % 绘制磁矩矩阵对于的像素图
%axis equal;             %横纵坐标单位长度相同
axis off;               %去掉坐标轴标度
set（gcf‘color‘[0.5 0.5 0.5]）;    %gcf为当前figure对象的句柄，
                                                     %设置颜色为灰色，及背景颜色.
                                                     %仅仅是为了好看

% up=round（rand（1L））*2-1;               %随机边界条件
% down=round（rand（1L））*2-1;
% zuo=round（rand（H1））*2-1;
% you=round（rand（H1））*2-1;

%up=ones（1L）*（-1）;                         %磁矩全部向下的边界条件
%down=ones（1L）*（-1）;
%zuo=ones（H1）*（-1）;
%you=ones（H1）*（-1）;

up=zeros（1