ワン・ランダウ法
次に示すのはPythonで実装されたワン・ランダウ法のサンプルコードである。遷移確率密度は対称であることを仮定している。
g ( x ′ → x ) = g ( x → x ′ ) {\displaystyle g({\boldsymbol {x}}'\rightarrow {\boldsymbol {x}})=g({\boldsymbol {x}}\rightarrow {\boldsymbol {x}}')}
このコードにおいて、対象となる系は "system" により表されるものとしている。
新たな配置を次の確率で採択/棄却する。 A ( r → r ′ ) = min ( 1 , e S − S ′ g ( r ′ → r ) g ( r → r ′ ) ) {\displaystyle A({\boldsymbol {r}}\rightarrow {\boldsymbol {r}}')=\min \left(1,e^{S-S'}{\frac {g({\boldsymbol {r}}'\rightarrow {\boldsymbol {r}})}{g({\boldsymbol {r}}\rightarrow {\boldsymbol {r}}')}}\right)} ここで、 S = S ( E ( r ) ) {\displaystyle S=S(E({\boldsymbol {r}}))} および S ′ = S ( E ( r ′ ) ) {\displaystyle S'=S(E({\boldsymbol {r}}'))} とする。
このステップを終えたとき、系の遷移先が E i {\displaystyle E_{i}} だったとすると H ( E i ) {\displaystyle H(E_{i})} を 1 だけ増やし、エントロピーを以下のように更新する。