fortran 计算非整数幂律

likm1110

fortran是不是不能计算非整数幂律？
如下：

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

Program main
  Implicit None
  Real :: a
  Real :: b
  Real :: x
  x = 3.0
  a = -1.0/74.0
  b = 1.1
  Write(*,*) a**b
End Program main

这个时候输出是NaN.为什么？
fortran怎么计算非整数幂律？就是我想用fortran把上面那个数算出来该怎么办？

likm1110

fcode 发表于 2015-11-10 20:33
由于浮点数是有误差的。所以指数很难确保是奇数。负数开方，在实数域内无解。你需要用复数来解决。

参考帖 ...

多谢！昨天也发现了这个问题的原因，我看看方法

fcode

由于浮点数是有误差的。所以指数很难确保是奇数。负数开方，在实数域内无解。你需要用复数来解决。

参考帖子：http://bbs.fcode.cn/thread-770-1-1.html