| 查看: 1091 | 回复: 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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 文献查阅 |
» 猜你喜欢
不自信的我
已经有5人回复
-
磺酰氟产物，毕不了业了！
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有16人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
★★★★★ 五星级,优秀推荐
★ ★ ★
woshilsh(金币+3,VIP+0):多谢交流，常来，互相帮助 2-20 16:32
woshilsh(金币+3,VIP+0):多谢交流，常来，互相帮助 2-20 16:32
|
楼主，我怎么运行有好多错误啊 --------------------Configuration: shuijishu - Win32 Debug-------------------- Compiling Fortran... D:\exe\shuijishu.for f90: Error: Illegal character in statement label field [s] f90: Error: Illegal character in statement label field f90: Error: Illegal character in statement label field f90: Error: Illegal character in statement label field [r] f90: Error: Illegal character in statement label field [o] f90: Error: First statement in file must not be continued D:\exe\shuijishu.for(1) : Error: Syntax error, found END-OF-STATEMENT when expecting one of: ( % . = => subroutine gauss(x1,x2,x,w,n) -----------------------------^ D:\exe\shuijishu.for(30) : Error: A RETURN statement is invalid in the main program. return ------^ D:\exe\shuijishu.for(3) : Error: An automatic object is invalid in a main program. [X] real x1,x2,x(n),w(n) -----------------^ D:\exe\shuijishu.for(3) : Error: An automatic object is invalid in a main program. [W] real x1,x2,x(n),w(n) ----------------------^ D:\exe\shuijishu.for(2) : Error: A specification expression object must be a dummy argument, a COMMON block object, or an object accessible through host or use association [N] integer n --------------^ Error executing df.exe. shuijishu.obj - 11 error(s), 0 warning(s) 我对fortran不是太懂，你可以给我的意见吗，谢谢 |
3楼2009-02-20 15:37:15
5楼2009-02-21 09:59:30