其实大家都没有什么好办法

波塞冬的信徒

背景是[Ti]=0.3%，[O]=0.005%，[N]=0.008%三元素，初始反应到1873K平衡态。计算过程的核心就是求解一个非线性方程组。

牛顿法迭代，就得自己写代码，求偏导，然后计算解线性迭代式，要么用克拉默法则，要么矩阵，还得循环，其实不容易的。

规划求解举例见网页，http://www.360doc.com/content/17/1008/19/38054235_693268454.shtml。
为保证物理意义，就是相关参数大于等于零，添加约束条件，精度过高，也无所谓了，最多就是死机，这代价可以忽略不计。
按照示例的步骤操作，一路顺风顺水。

解出来，跟作者的结果对比，误差还挺大的。

Architect

多试几种解决方案，追踪计算步数及中间步结果。
如果原有算例本身就有问题，那就换一个