这个Ising模型通过应用Metropolis算法-蒙特卡罗方法来模拟磁系统(正、负或随机自旋)。运行主文件，输入晶格尺寸(最好是100)，然后选择一个输入自旋作为初始配置。设置两种不同的温度(T=2.0和T=2.5)。

例如，在T=2(低温)初始化一个选择的正自旋时，大多数自旋是黑色的，这是因为自旋翻转的几率很小，而且材料具有铁磁性。对于T=2.5的情况，由于温度较高，它产生了自旋翻转的趋势。因此，材料失去了磁化。自旋似乎没有对齐，由此产生的构型似乎是随机无序的。这是因为顺磁行为。

模拟的下一部分是可观测值的计算:平均磁化强度、平均能量、平均磁化率和比热。为了计算平均能量和磁化,我们必须找到一个时间依赖的能源和磁化强度值变化小(能源和磁化强度的时间变化随着时间的增加)因此,我们选择精度p和检查时间间隔(时间步的数量应满足精度)。这些间隔应该取决于初始配置。每次步长比较当前步长变量和以前步长变量(平均能量和平均磁化强度)，并与精度进行比较。如果它满足整个区间的条件，我们将当前的时间步长称为近似的热化时间，并开始计算变量，平均磁化强度，平均能量，平均磁化率和比热。