The performance of current solutions for nonlinear equations, such as N-R and N+1 strongly depends on precision and initial values of model, which causes unstableness of convergence.

• 目前采用的牛顿-拉夫逊法和N+1残量法对方程组的求解依赖于模型的精度和初猜值，收敛性不稳定。