24小时热门版块排行榜

返回列表

【奖励】本帖被评价4次，作者学员7eU8hQ增加金币 3.5 个

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zsuzwr

木虫 (正式写手)

应助: 2 (幼儿园)
贵宾: 0.15
金币: 2811.7
帖子: 554
在线: 275小时
虫号: 209571

[资源] 【分享】产生高斯分布随机数源代码(Fortran)

subroutine gauss(x1,x2,x,w,n)
   integer n
   real x1,x2,x(n),w(n)
   double precision xm,xl,z,p1,p2,p3,pp,z1
   integer m,i,j
   parameter (eps=3.d-14)
   m=(n+1)/2
   xm=0.5d0*(x2+x1)
   xl=0.5d0*(x2-x1)
   do 12 i=1,m
      z=cos(3.141592654d0*(i-.25d0)/(n+.5d0))
1    continue
      p1=1.d0
      p2=0.d0
      do 11 j=1,n
         p3=p2
         p2=p1
         p1=((2.d0*j-1.d0)*z*p2-(j-1.d0)*p3)/j
11       continue
      pp=n*(z*p1-p2)/(z*z-1.d0)
      z1=z
      z=z1-p1/pp
      if(abs(z-z1).gt.eps)go to 1
      if (abs(z).lt.eps) z=0.
      x(i)=xm-xl*z
      x(n+1-i)=xm+xl*z
      w(i)=2.d0*xl/((1.d0-z*z)*pp*pp)
      w(n+1-i)=w(i)
12 continue
   return
   end

[ Last edited by woshilsh on 2009-2-20 at 16:32 ]

回复此楼