| 查看: 381 | 回复: 0 | |||
[交流]
【求助】程序是产生随机数,请看以下问题，谢谢！
|
|
程序是产生随机数的,f90，我用fortran 6.6运行出现 Access Violation 好像是IDUM=1的问题 FUNCTION RAN1(IDUM) DIMENSION R(97) PARAMETER(M1=259200,IA1=7141,IC1=54773,RM1=1./M1) PARAMETER(M2=134456,IA2=8121,IC2=28411,RM2=1./M2) PARAMETER(M3=243000,IA3=4561,IC3=51349) DATA IFF/0/ IF(IDUM.LT.0.OR.IFF.EQ.0) THEN IFF=1 IX1=MOD(IC1-IDUM,M1) IX1=MOD(IA1*IX1+IC1,M1) IX2=MOD(IX1,M2) IX1=MOD(IA1*IX1+IC1,M1) IX3=MOD(IX1,M3) DO 11 J=1,97 IX1=MOD(IA1*IX1+IC1,M1) IX2=MOD(IA2*IX2+IC2,M2) R(J)=(real(IX1)+real(IX2)*RM2)*RM1 11 CONTINUE IDUM=1.0 ENDIF IX1=MOD(IA1*IX1+IC1,M1) IX2=MOD(IA2*IX2+IC2,M2) IX3=MOD(IA3*IX3+IC3,M3) J=1+(97*IX3)/M3 IF(J.GT.97.OR.J.LT.1)PAUSE RAN1=R(J) R(J)=(real(IX1)+real(IX2)*RM2)*RM1 RETURN END function |
回复此楼» 猜你喜欢
278求调剂
已经有9人回复
-
297求调剂
已经有12人回复
-
0703化学297求调剂
已经有3人回复
-
材料学硕301分求调剂
已经有6人回复
-
280求调剂
已经有11人回复
-
311求调剂
已经有3人回复
-
求调剂
已经有3人回复
-
313求调剂
已经有4人回复
-
317求调剂
已经有15人回复
-
307求调剂
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
请教Windows系统下perl脚本的ARGV应用问题，谢谢
已经有5人回复
- fortran 随机数 每次运行都是一样的！
已经有14人回复
- 请教一个fortran小程序编译出错的问题，谢谢
已经有9人回复
- 产生高斯分布随机数的原理
已经有3人回复
- 求助！！请高手看看这段英文有什么问题，帮忙改改，谢谢！！
已经有3人回复
- 有关fortran的一次而问题，希望大家能帮帮忙，谢谢
已经有4人回复
- 中国物理快报审稿问题，懂的请看一下，谢谢！
已经有19人回复
- 【求助】如何用汇编语言产生随机数
已经有5人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
