tecsun-hlj 发表于 2007-12-16 19:51

fortran编程计算分叉图遇到困难，请大侠指点

我用用fortran编程计算分叉图，有非线性油膜力
运行总是出现：run-time error M6201:math
                      -sqrt: DOMAIN error

我定义时使用：DOUBLE PRESICION,用的函数是：DSQRT
请指点一下哪里出现了问题，使用开方函数应该注意什么
谢谢！急呀

[ 本帖最后由 无水1324 于 2007-12-19 10:22 编辑 ]

octopussheng 发表于 2007-12-17 13:42

请看fortran的帮助

Examples

SQRT (16.0) has the value 4.0.

SQRT (3.0) has the value 1.732051.

The following shows another example:

! Calculate the hypotenuse of a right triangle
! from the lengths of the other two sides.
REAL sidea, sideb, hyp
sidea = 3.0
sideb = 4.0
hyp = SQRT (sidea**2 + sideb**2)
WRITE (*, 100) hyp
100 FORMAT (/ ' The hypotenuse is ', F10.3)
END

tecsun-hlj 发表于 2007-12-17 19:56

首先谢谢回帖
fortran帮助我看了，好像用法没什么问题,我把sqrt函数内全变为 正的依然出现这个错误提示
搜了一下，看到这是个老问题但没看到有效的解决办法，还请达人指点

octopussheng 发表于 2007-12-18 08:06

变量定义用real*8试试呢！

tecsun-hlj 发表于 2007-12-18 11:52

一样的，我试过了

octopussheng 发表于 2007-12-18 14:54

程序定义大致贴一下吧！这样有点盲人摸象的感觉，呵呵！
页: [1]
查看完整版本: fortran编程计算分叉图遇到困难，请大侠指点