argy.f

来自「网络带宽测试工具」· F 代码 · 共 25 行

      Subroutine argy( a, b, n, y )
! ----------------------------------------------------------------------
! ---  n random numbers in interval (a,b) are stored in array y
!      not sorted.
! ----------------------------------------------------------------------
      Use                numerics
      Implicit           None
 
      Integer         :: n
      Real(l_)        :: a, b, y(n)

      Integer         :: i
      Real(l_)        :: dely, dran1

      Integer         :: iy
      Common /ran1iy/    iy
      External           dran1
! ----------------------------------------------------------------------
      dely = b - a
      Do i = 1, n
         y(i) = a + dely*dran1( iy )
      End Do
! ----------------------------------------------------------------------
      End Subroutine argy