Drand48Gam Drand48Gam_sub GenrandGam GenrandGam_sub RanfArrGam RanfArrGam_sub

GenrandGam_sub() Sub GenrandGam_sub ( Ret As  Double, Optional Alpha As  Double = 1, Optional Beta As  Double = 1  ) 53 ビット実数ガンマ乱数 (Marsaglia and Tsang の圧搾法) (メルセンヌ・ツイスター) (サブルーチン版) 目的 ガンマ分布 Γ(α, β) に従うガンマ乱数 (53 ビット実数) を生成する. メルセンヌ・ツイスターの乱数生成ルーチン GenrandRes53() により一様乱数を生成し, Marsaglia and Tsang の圧搾法によりガンマ乱数に変換する. 引数 [out] Ret ガンマ乱数 (53 ビット実数). Alpha あるいは Beta が正しくないときは0を返す. [in] Alpha (省略可) ガンマ分布 Γ(α, β) のパラメータ α(形状母数). (α > 0) (省略時 = 1) [in] Beta (省略可) ガンマ分布 Γ(α, β) のパラメータ β(尺度母数). (β > 0) (省略時 = 1) 参考文献 Mersenne Twister Home Page 四辻哲章「確率分布乱数生成法」第2版, プレアデス出版, 2019.