(981) LJ to Morse MCS LJ近似曲線（Metropolis,3変数） [
zip,
source ]
メトロポリス法を用いて、Lennard-Jones型ポテンシャルに近いMorse型のポテンシャルを
求めてみます。近さの指標として err = ∫weight(r)*(Morse(r)-LJ(r))^2 dr とします。
メトロポリス法を適用し、温度kTをだんだん下げていくことでerrが最小になるD,r0,Aを
求めます。kTを下げていくとき、sを小さくして、適当な採択率を維持してください。
(created 2006.07.22, last updated 2007.11.24)