Chinaunix首页 | 论坛 | 博客
  • 博客访问: 1105271
  • 博文数量: 264
  • 博客积分: 7225
  • 博客等级: 少将
  • 技术积分: 5096
  • 用 户 组: 普通用户
  • 注册时间: 2008-11-17 08:53
文章分类

全部博文(264)

文章存档

2011年(33)

2010年(52)

2009年(152)

2008年(27)

我的朋友

分类: LINUX

2009-08-11 15:18:03

【分享】产生高斯分布随机数源代码(新版)
★ ★
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) |
给主人留下些什么吧!~~