容我试验一下代码的使用哈

likm1110 · 发表于 2015-3-13 09:13:22

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

01Module ran001
02   Implicit None
03   Save
04   Integer :: n = 9876
05End Module ran001
06!-------------------------------------------------------------
07Subroutine random0(ran)
08   Use ran001
09   Implicit None
10   Real, Intent(out) :: ran
11   n = Mod(8121 * n + 28411, 134456)    ! A number between 0 and 134455
12   ran = real(n) / 134456.      ! Generate a number [0,1.)
13End Subroutine random0
14!-------------------------------------------------------------
15Subroutine seed(iseed)
16   Use ran001
17   Implicit None
18   Integer, Intent(in) :: iseed
19   n = abs(iseed)
20End Subroutine seed