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）;??????%?繪制磁矩矩陣對(duì)于的像素圖
%axis?equal;?????????????%橫縱坐標(biāo)單位長(zhǎng)度相同
axis?off;???????????????%去掉坐標(biāo)軸標(biāo)度
set（gcf‘color‘[0.5?0.5?0.5]）;????%gcf為當(dāng)前figure對(duì)象的句柄，
?????????????????????????????????????????????????????%設(shè)置顏色為灰色，及背景顏色.
?????????????????????????????????????????????????????%僅僅是為了好看

%?up=round（rand（1L））*2-1;???????????????%隨機(jī)邊界條件
%?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