【分享】产生高斯分布随机数源代码(新版)
★ ★
kuhailangyu(金币+2,VIP+0):非常好,谢谢提供,再接再厉!
kuhailangyu(金币+2,VIP+0):非常好,谢谢提供,再接再厉!
kuhailangyu(金币-2,VIP+0):多按了一次,要回两个,呵呵
Function Gauss(Dummy)
Implicit Double Precision (a-z)
integer i
Parameter ( A1 = 3.949846138, A3 = 0.252408784 )
Parameter ( A5 = 0.076542912, A7 = 0.008355968 )
Parameter ( A9 = 0.029899776 )
! *******************************************************************
Sum = 0.0
Do i = 1, 12
! Call Random_number (ran)
Sum = Sum + Ranf(Dummy)
Enddo
R = (SUM - 6.0 )/4.0
R2 = R * R
Gauss = ((((A9 * R2 + A7) * R2 + A5) * R2 + A3) * R2 +A1)*R
END
阅读(1019) | 评论(0) | 转发(1) |