开始写 non eq msm 的验证程序了。

先最简单的平衡态体系。考虑了下单位问题。用 Euler 格式实现。写了分布统计程序。想了想扰动应该怎么加，时间尺度问题。

发现在格式上有不少改进空间。