| 查看: 1090 | 回复: 7 | |||
| 【奖励】 本帖被评价4次，作者学员7eU8hQ增加金币 3.5 个 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[资源]
【分享】产生高斯分布随机数源代码(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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 文献查阅 |
» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有12人回复
-
2025年遐想
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
-
求个博导看看
已经有18人回复