Fortran 90 产生不正确的答案，当答案应该是 1 时产生 0。我忽略了什么吗？ [复制]

Question

这是我编写的 Fortran 90 代码。这是我在较大代码中使用的函数，但得到的结果不正确。

我从数学上知道答案应该是 1（或由于计算而接近 1），但我的答案结果是 0。

我做错了吗？有什么我不知道的吗？

这是我保存为 test.f90 的代码，我正在 Windows gfortran 上编译和运行：

program main
    implicit none

    real*8 :: y, t

    t = 0.0
    y = (1/3)*exp(5*t) + (2/3)*exp(-t) + t*t*exp(2*t)

    print*, 'y= ', y

end program main

Answer 1

fortran 认为 1/3 和 2/3 等于 0。你应该在你的数字上加一个点：

program main
implicit none

real*8 :: y, t

t = 0.0
y = (1./3.)*exp(5*t) + (2./3.)*exp(-t) + t*t*exp(2*t)

print*, 'y= ', y