作者microcat (微喵)
看板Fortran
標題[問題] fortran產生隨機亂數的方法?
時間Thu Mar 12 17:06:41 2009
請問如何用 fortran 產生 常態分配 及 卡方分配 的亂數
normal chi-square
以前用統計軟體甚是方便,突然間要用 fortran ,
熊熊不知道該如何產生?
參考及修改了別人的程式,依舊無法執行,
請多幫忙了,謝謝。
已修正的程式:(感謝大家幫忙!)
可正常執行!
Program CHI
USE IMSL
IMPLICIT NONE
INTEGER ISEED, NOUT, NR
REAL DF, R(5)
CALL UMACH (2, NOUT)
DF = 5.0
NR = 5
ISEED = 0
CALL RNSET (ISEED)
CALL RNCHI (DR, DF, R)
WRITE (NOUT,99999) R
99999 FORMAT (' Chi-squared random deviates with 5 df: ', 5F7.3)
END Program CHI
http://tkustat.stat.tku.edu.tw/randgen.php
淡江的網站還挺好用的!
不知道能不能寫出一樣好用的程式?!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.109.72.144
※ 編輯: microcat 來自: 140.109.72.144 (03/12 17:17)
1F:→ charlesdc:先看你的編譯器有沒有支援吧 03/12 21:13
2F:推 Cypresslin:有副程式或是函數庫嗎? 03/12 22:04
3F:→ microcat:請問要如何知道編譯器的版本?謝謝 03/13 00:06
4F:推 jsb:這個是 IMSL 的副程式。如果你是用M$的編譯器, 03/13 09:47
5F:→ jsb:在 Program XXX 的下一行加個 USE MSIMSL 就行。 03/13 09:48
6F:→ microcat:謝謝各位的指導 03/13 10:44
※ 編輯: microcat 來自: 140.109.72.144 (03/13 10:51)