| 查看: 697 | 回复: 5 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 anxt2006 的 5 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
anxt2006木虫 (正式写手)
|
[交流]
[求助]fortran里如何产生从0到1的加权的随机数？
|
||
fortran里如何产生从0到1的加权的随机数？我想用蒙特卡罗积分，但是我的被积函数有尖峰，所以想找一个加权的随机数产生程序，在尖峰的地方多取一些点，在其他地方少取一些点。请高手指教。 [ Last edited by kuhailangyu on 2008-12-9 at 10:32 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
国社科又开始会评了，不知道这次命运如何
已经有13人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有13人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
maomao1210
金虫 (正式写手)
- 程序强帖: 5
- 应助: 2 (幼儿园)
- 金币: 1432.8
- 散金: 242
- 红花: 16
- 沙发: 1
- 帖子: 991
- 在线: 442.1小时
- 虫号: 253215
- 注册: 2006-05-20
- 性别: MM
- 专业: 考古理论
|
SUBROUTINE(JIE,X) IMPLICIT REAL*8 (A-H,O-Z) INTEGER*4 I(2) C READ(*,*) JIE IF (JIE.EQ.0) GO TO 1 I(1)=0 I(2)=JIE 1 IX=1490*I(2)+3805*I(1) IE=3805*I(2) I(2)=IE-IE/32768*32768 IX=IX-IX/32768*32768+IE/32768 I(1)=IX-IX/32768*32768 X=(I(1)*32768+I(2))*9.3132258D-10 END RETURN END 这样就行了，你调用该程序的时候，给定一个种子JIE，然后这个程序就会给你一个数值X。 [ Last edited by maomao1210 on 2008-10-14 at 22:15 ] |
3楼2008-10-13 22:14:13
juloong
木虫 (职业作家)
X射线衍/散射测试
- 应助: 103 (高中生)
- 金币: 6402.5
- 散金: 1920
- 红花: 100
- 帖子: 3142
- 在线: 2054.1小时
- 虫号: 275864
- 注册: 2006-09-02
- 专业: 半导体材料
2楼2008-10-13 18:21:32
maomao1210
金虫 (正式写手)
- 程序强帖: 5
- 应助: 2 (幼儿园)
- 金币: 1432.8
- 散金: 242
- 红花: 16
- 沙发: 1
- 帖子: 991
- 在线: 442.1小时
- 虫号: 253215
- 注册: 2006-05-20
- 性别: MM
- 专业: 考古理论
4楼2008-10-14 08:33:04
anxt2006
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 6103.3
- 散金: 38
- 红花: 1
- 帖子: 381
- 在线: 208.2小时
- 虫号: 578555
- 注册: 2008-07-04
- 专业: 凝聚态物性 II ：电子结构
5楼2008-10-14 15:20:44
30