[問題] 常態分配
老闆要求開始跑模擬了,可是fortran是我五年前學過的
而且只學過初級的 > <
這次我要模擬出mean =10 20 30 40 50 標準差 = 1 的常態分配
X~N(10,1), N(20,1), N(30,1),N(40,1),N(50,1)各要模擬5個值
每組要10次
可是我現在只能做出mean=10
可是不知道該怎麼在寫出個迴圈跑mean
這個程式是之前有位版大寫的
Program normal
use msimsl
IMPLICIT NONE
integer I,j,M,NR,ISEED,NOUT
real R(10),RAND(5),XM,XSTD
DO 100 I=1,10 !→產生樣本的組數為M 組
DO 50 J=1,5 !→產生每組有5 個隨機常態的樣本點
NR =5
XM=10
XSTD=1
CALL RNSET(ISEED)
CALL UMACH(2,NOUT)
R(I) =RNNOF()
CALL RNGET(ISEED)
CALL RNNOA (NR, R)
CALL SSCAL (NR,XSTD,R,1)
CALL SADD (NR,XM,R,1)
RAND(J)=R(5)
50 CONTINUE
WRITE (NOUT,99999) RAND
99999 FORMAT ('Normal Random Deviates: ', 5F8.4)
100 CONTINUE
end Program normal
請問我該怎麼再加一個迴圈跑mean呢?
麻煩大家幫我一下,謝謝~~
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.75.100
Fortran 近期熱門文章
PTT數位生活區 即時熱門文章
-1
12