GAMMA 的结果低于同类

Question

我想使用下面的程序计算 gamma(-170.1)：

program arithmetic  
! program to do a calculation  
real(8) :: x  
x = GAMMA(-170.1)  
print *, x  
end program  

但出现错误：

测试.f95:4.10:

x = 伽马(-170.1) 1 错误：GAMMA 的结果在 (1) 处下溢

错误：当我使用 gfortran 进行编译时， 。根据 Maple gamma(-170.1) = 5.191963205*10^(-172) ，我认为应该在我定义的变量 x 的指数范围内。

Answer 1

对您的程序进行以下修改应该可以工作。请记住，在 Fortran 中，在分配给 LHS 之前先计算 RHS，并且浮点文字是默认类型，即单精度。因此，使 GAMMA 的参数成为双精度，编译器会选择双精度 GAMMA。

program arithmetic  
! program to do a calculation  
integer, parameter :: dp = kind(1.0d0)
real(dp) :: x  
x = GAMMA(-170.1_dp)  
print *, x  
end program

Answer 2

-170.0 可以被视为浮点数。如果是这样，将其更改为双精度应该可以解决问题。