乱数発生プログラム fortran

研究

めちゃくちゃ簡単だけど乱数発生プログラムをfortranで書いた。いわゆる性質の悪いやつ。
シミュレーションの初期条件くらいだったらこのくらいでいいと思われる。いままでrand()使ってたけど、なんとなく
使った式は、
A=1103515245
B=12345
M=2147483647
として、以下の漸化式を解くということ、

X_(n+1)=(A*X_n+B)mod M

program test
implicit real(8) (a-h, o-z)
integer(4), parameter :: ix = 100
integer(4), dimension(ix) :: ii,jj

ic=1103515245
jc=12345
mc=2147483647

ii(1) = 1
jj(1) = 2

do i = 1,ix-1
ii(i+1) = mod*1/dble(mc),dble(jj(i+1))/dble(mc)
enddo

end program test

ii,jjは違う種で始めたということです。

*1:ic*ii(i)+jc),mc) jj(i+1) = mod((ic*jj(i)+jc),mc) write(*,*) dble(ii(i+1