最优控制 + 化学势 双双崩溃中。。。。

最优控制终于算好了一维和一维线性变化的例子。。。许久没有做这样的数值了，感觉有些生疏。。。。另外 MD 真是慢。。。butane 还是有差距，在想到底是什么原因。。


化学势崩溃中。。。设计了几个“重水”，之前的理论发现严重不靠谱。。。木有什么新思路。做几个疯狂的算例看看，比如：
1，1，1
16，16，16
16，16，1
16，16，0.1