二维伊辛模型的Matlab代码

二维伊辛模型的Matlab代码，可以给学习蒙特卡洛算法的同学们练手

function y = Cal_E（H）
global L;
global SPIN;
E = 0;
for i = 2:L+1
    for j = 2:L+1
        E =-（SPIN（ij）*SPIN（i-1j）+SPIN（ij）*SPIN（i+1j）+SPIN（ij）*SPIN（ij+1）+SPIN（ij）*SPIN（ij-1）+SPIN（ij）*SPIN（i+1j+1）+SPIN（ij）*SPIN（i-1j-1））/2-SPIN（ij）*H+E;
    end
end
y=E;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        276  2013-10-22 22:48  Monte Carlo Ising Model\Cal_E.m

     文件        282  2013-10-22 22:54  Monte Carlo Ising Model\Cal_E1.m

     文件        371  2013-10-22 22:48  Monte Carlo Ising Model\Flip_Spin.m

     文件       3292  2013-10-23 04:08  Monte Carlo Ising Model\main.m

     文件        756  2013-10-23 03:39  Monte Carlo Ising Model\Monte_Carlo.m

     文件        365  2013-10-23 02:48  Monte Carlo Ising Model\Monte_Carlo_warmup.m

     目录          0  2014-08-23 17:07  Monte Carlo Ising Model

----------- ---------  ---------- -----  ----

                 5342                    7