sqrt 定义域出错的问题

wansq

还是这段老程序，出了好几次错了，出一次改一次，又才出一次。。。运行出来，是开方的问题，但是不知道具体是哪里，怎么改啊，求助求助！

楚香饭

检查一下 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

chuxf 发表于 2014-3-28 13:15
检查一下 144，145 行。
[mw_shl_code=fortran,true]PSI=2*ASIN((SQRT(XXI-XX(J))**2+(YYI-YY(J))**2)/(2*R ...

再次谢谢你，我改过来了

wansq

chuxf 发表于 2014-3-28 13:15
检查一下 144，145 行。
[mw_shl_code=fortran,true]PSI=2*ASIN((SQRT(XXI-XX(J))**2+(YYI-YY(J))**2)/(2*R ...

请问一下，您是经过调试，查找到这一行的是吧

楚香饭

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

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

wansq

chuxf 发表于 2014-3-28 13:33
是的。直接调试，出错后就自动停在这一行

哦，明白了，如果下次出问题，我试着调试一下，就知道问题出在哪里了