求二元非线性函数的极值或者根，传统数值优化就有很多解法吧。如果是等价于求合力为0的点，可以Newton迭代，Broydent迭代），或者求总势能最低的点，这个可以用BFGS迭代算法，等价于求Hessian矩阵的近似值来得到爬山方向。