sqrt 定义域出错的问题

楚香饭

检查一下 144，145 行。

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

PSI=2*ASIN((SQRT(XXI-XX(J))**2+(YYI-YY(J))**2)/(2*RO1(J)))
      CSI=2*ASIN((SQRT(XXXI-XXX(J))**2+(YYYI-YYY(J))**2)/(2*RO2(J)))

你确定是要对 XXI-XX(j) 先开根号，再平方？？
而不是平方和，然后开根号？

楚香饭

wansq 发表于 2014-3-28 13:29
请问一下，您是经过调试，查找到这一行的是吧

是的。直接调试，出错后就自动停在这一行