子例行程序

胆怯滴大侠

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

program hello
    implicit none
    integer n,s

    subroutine fcount(n,s)
    s=0
    do i=1,n
        s=s+1./i
    end do
    end

    call fcount(50,s1)
    call fcount(100,s2)
    call fcount(150,s3)
    call fcount(200,s4)
    s=s1+s2+s3+s4
    print *, "s=", s

end program

求助，请教大佬们，subroutine fcount(n,s)这里报错了，这个问题怎么解决呢？

vvt

程序单元里不能直接放入另一个程序单元，需要符合一定的规则。

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

program hello
  implicit none
  real :: ss,s(4)

  call fcount(50,s(1))
  call fcount(100,s(2))
  call fcount(150,s(3))
  call fcount(200,s(4))
  ss=sum(s)
  print *, "s=", ss

contains

subroutine fcount(n,s)
  integer :: n , i
  real :: s
  s=0
  do i=1,n
    s=s+1./i
  end do
end subroutine fcount

end program hello

胆怯滴大侠

vvt 发表于 2021-4-18 16:19
程序单元里不能直接放入另一个程序单元，需要符合一定的规则。

[mw_shl_code=fortran,true]program hello

非常感谢

vvt

考虑到 Fortran 是一门向量化的语言，给出一个高级的版本，非常简单且高效。

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

program hello
  implicit none
  write(*,*) "s=", sum(fcount([50,100,150,200]))

contains

Elemental Real Function fcount(n) result( s )
  integer , intent(IN) :: n
  integer :: i
  s=sum(1./[(i,i=1,n)])
End Function fcount

End program hello