高斯牛顿法，迭代无法终止（死循环）

楚香饭

在我这里，输入 2，3，4，5，6 也一样无法跳出循环。ei 越来越大，最后变成无穷大

楚香饭

你用的什么编译器？我试过 ivf，cvf，gfortran，都无法跳出循环。

楚香饭

OK ，改了以下代码 就可以了

[Fortran] 纯文本查看 复制代码

REAL*8::G(row,col),Gt(col,row),GtG(col,col)=0.0D0,GtGl(col,col)=0.0D0,ds(row)=0.0D0,dl(row)=0.0D0,d(row)=0.0D0,&
Gtd(col)=0.0D0,Gdm(row)=0.0D0,u(col,col)=0.0D0,v(col)=0.0D0,s(ka)=0.0D0,e(ka)=0.0D0,work(ka)=0.0D0,dd(row)=0.0D0

变量要初始化。