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

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

我用用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

程序定义大致贴一下吧！这样有点盲人摸象的感觉，呵呵！

账号		自动登录	找回密码
密码			我要加入

[稳定性与分岔] fortran编程计算分叉图遇到困难，请大侠指点

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。