二维正方晶格Ising模型的计算机模拟

在无外场的情况下，自由边界的方格上，存在着某种粒子，它具备三种自旋状态-1,0,1。我们认为t_initial=200后系统达到稳定态。系统演化至t_max=1200步。（1）分别绘制单粒子平均磁矩、单粒子平均能量与温度之间的关系图。观察并讨论图中现象。（2）根据（1）中的现象，合理选取不同的T，分别绘制系统达到稳定态后的演化过程中，系统总磁矩的分布曲线