Fortran只给出第一个任期的结果？

Question

我正要计算cos(x)+1/4*cos(2x),但结果总是给我结果cos(x).我的代码中的错误在哪里？

  program write
  implicit none
  integer, parameter :: N=8
  integer :: j
  real :: h, L
  real, dimension(0:N-1) ::x, fx
  real(8), parameter :: pi=4.0_8*atan(1.0_8)
  L=2*pi

  h=L/N
  do j=0,N-1
     x(j)=h*j
  end do
  do j=0,N-1
     fx(j)=cos(x(j))+1/4*cos(2*x(j))
  end do

  write(*,*),fx


  end program write

Answer 1

你的问题是1/4.

由于两个1和4是整数,1/4被解释为一个整数除法,任何提醒被丢弃.简而言之:1/4 == 0,而1.0/4 == 1/4.0 == real(1)/4 == 0.25.

请注意real(1/4) == real(0) == 0.0.