Ising Model产品系统

伊辛模型由晶格上的N个自旋Si组成（其中i=1,...,N），
每个都可以指向上，Si = 1，或指向下，Si = -1。每对相邻的
自旋对齐会使系统能量降低J>0的量。

给定一个自旋配置{Si}，总能量为E({Si}) = −JSiSjij其中，该和是对所有不同的最近邻对ij求和。根据玻尔兹曼分布，即观察到给定构型的概率{Si}在平衡状态下的概率P({Si})等于exp[−E({Si})/(kBT)]，系统的物理性质由J与kBT的比值决定（该因子会进入玻尔兹曼权重），我们可以取J、kB和T中的任意两个等于1，而不失一般性——使用J = kB = 1进行模拟

使用Python对伊辛模型的平衡态进行采样能量，磁化，以及相关的统计量进行分析。