来自朴柱宏的问题
fortran77中为什么不能像正数那样求负数的三次方根,该怎么求?
fortran77中为什么不能像正数那样求负数的三次方根,该怎么求?
1回答
2020-03-17 16:29
fortran77中为什么不能像正数那样求负数的三次方根,该怎么求?
fortran77中为什么不能像正数那样求负数的三次方根,该怎么求?
主要因为幂函数x**a是定义在[0,+无穷)上的,在Fortran77中除了特殊的指数a之外x**a一般都按照exp(a*log(x))来实现.如果x.GE.0.0D0,那么x**(1.0D0/3.0D0)当然没什么问题.如果x.LT.0.0D0,Fortran是不会知道x**(1.0D0/3....