run-time error M6021:MATH - log:DOMAIN error

楚香饭

楼主，你好。

建议长代码直接上传 *.f90 或 *.for 或 *.f 源代码文件。这样大家好告诉你是第几行的问题。（复制粘贴后，有时候行数不好对应）

根据我的调试，你的问题是，主程序中，Addterm 没有初值。

在第191行，IF (STEP .NE. 1.D-4) THEN 中，STEP ==1.D-4，所以 NE 条件不满足，直接跳到 ENDIF 执行，即217行。虽然在195行对Addterm有给初值，但没有执行。

217行的循环执行完毕后，执行223行的
XN(K)=XN(K)+STEP*(0.5*SUM1-XN(K)*SUM2+XJK*ZETA(K)+ADDTERM-SUBTERM)
时，Addterm 没有初值，于是出现了极其不正常的数据。然后到后面就错了。