一番目の格子と周りの格子との間の交換相互作用を計算します。
但し、系には周期的境界条件がかけられているので格子点１の左隣は３番の格子点、下は７番の格子点 となります。 ちなみにこの場合
E_before=-4J となります。

次に、１番目の格子点をひっくり返した時の交換相互作用によるエネルギーを 計算します。 ちなみにこの場合
E_after=2
です。

格子点１をひっくり返す後前のエネルギー差を計算します。
ΔE=E_after-E_before
ちなみにこの場合は
ΔE=4
次にこの値をボルツマン因子に代入します。 ちなみにこの場合は P=exp(-4/kT) となります。

０から１の乱数を発生させてボルツマン因子の方が乱数より大きい値 をとるときにスピンをひっくり返します。

ちなみにボルツマン因子とエネルギーの関係は下の図のようになります。