应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【交流】Fortran语言答疑专帖
北京石油化工学院2026年研究生招生接收调剂公告
15811/16返回列表首页上一页91011121314下一页末页
查看: 11518  |  回复: 157
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 程序强帖 ，点击这里进行查看

jihe01

木虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
微尘、梦想(金币+2): 谢谢回贴交流…… 2011-04-19 15:13:08
看了maomao1210开的这个帖子，很受启发，对Fortran学习帮助很大。有个程序问题请教一下，这是网上一个产生随机数的程序，用95编写的，当我在调用FUNCTION时（用6.6版本），出现Access volidation错误，问题出在11 continue后面一行  IDUM=1.0上，但我不明白什么原因，请帮忙看一下，谢谢！


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
赞一下(2人) 回复此楼
101楼2011-04-11 15:24:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)
★ ★
微尘、梦想(金币+2): 谢谢回复！ 2011-04-19 15:13:35
引用回帖:
Originally posted by jihe01 at 2011-04-11 15:24:00:
看了maomao1210开的这个帖子，很受启发，对Fortran学习帮助很大。有个程序问题请教一下，这是网上一个产生随机数的程序，用95编写的，当我在调用FUNCTION时（用6.6版本），出现Access volidation错误，问题出在11 ...

我刚测试了一下（用的是你说的vf6.6），结果一切都不错，估计程序不是出现这个函数中吧，是你操作有问题，或者其他模块的程序有问题。打开资源管理器，把相应的exe关掉，再运行一下，看是否还会出现类似的错误。
赞一下(1人) 回复此楼
102楼2011-04-11 15:52:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jihe01

木虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
dubo(金币+1): 欢迎讨论 2011-08-27 16:39:20
谢谢maomao的测试，可能确实不是这个function的问题，可能是其他程序调用的问题，这个程序有点复杂，我是从网上下的，是关于MCMC方法计算的，当我调试时，总是到这个地方出现措施（ＩＤＵＭ＝１）。附件是原程序和错误提示，请看一下，谢谢。
赞一下(2人) 回复此楼
103楼2011-04-11 19:57:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qwer_2009

铜虫 (小有名气)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
微尘、梦想(金币+2): 欢迎常来程序语言版，鼓励新虫多发话题！ 2011-04-19 15:15:09
program main
implicit none
type node
real x   !设置随机数
type(node),pointer::next
end type
type(node),pointer::head,q,p
real rand
integer i,j
integer::m=0,n=0


call random_seed( )
nullify(head)


do i=1,10
call random_number(rand)
allocate(q)
q%x=rand
q%next=>head
head=>q
end do

q=>head
do while(associated(q))
print*,q%x
q=>q%next
end do

!判断是否大于0.50
q=>head

do while(associated(q))


if(q%x<0.50.and.associated(q,head))then
n=n+1  !判断小于0.5的个数
head=>q%next
deallocate(q)
q=>head
else if(q%x<0.50.and..not.associated(q,head))then
n=n+1
p=>q%next
deallocate(q)
q=>p
else
m=m+1
q=>q%next
end if

end do


write(*,*) "大于0.5的个数m=", m
write(*,*) "小于0.5的个数n=", n

!输出大于0.5的数
q=>head
do while(associated(q))

print*,q%x
q=>q%next
end do

end

请大家帮忙看看这个程序问题出在哪里了，在输出m和n的值时还是正确的，为什么加上最后几句 显示大于0.5的数时就会无限循环下去，而且得到的结果还是错误的呢，谢谢
赞一下(2人) 回复此楼
104楼2011-04-19 11:06:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abc__1203

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
可否提供一段Fortran调用MATLAB函数的实例程序
例如想调用MATLAB里边的besselj函数
谢谢！
赞一下(1人) 回复此楼
105楼2011-05-16 21:29:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suxingtao

铁杆木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
dubo(金币+1): 欢迎讨论 2011-08-27 16:40:22
我有个Fortran插值问题，http://muchong.com/bbs/viewthread.php?tid=3305448&fpage=1，向您请教！
赞一下(2人) 回复此楼
106楼2011-06-16 14:43:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿水哥王

铜虫 (初入文坛)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
dubo(金币+1): 欢迎讨论 2011-08-27 16:40:11
楼主你好  我目前在算一个球星量子点矩阵的本征值和本征矢  我想问下这个矩阵元的积分里怎样编译这个沿径向函数的积分元
赞一下(2人) 回复此楼
沉默越久 爆发越大
107楼2011-06-20 15:42:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shaofain

铁虫 (初入文坛)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
jjdg(金币+1): 感谢参与 2011-08-01 17:02:39
楼主你好：能不能推荐一个关于fortran工程算法的网站或是书？谢谢
赞一下(2人) 回复此楼
108楼2011-08-01 15:24:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anshenyu

铜虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
dubo(金币+1): 欢迎讨论 2011-08-27 16:39:57
请教下fortran能不能处理这样的矩阵代数方程?
P=(I-PH^TH/r)(APA^T+BB^Tq)
其中大写字母均代表矩阵
赞一下(2人) 回复此楼
109楼2011-08-25 16:19:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
微尘、梦想(金币+1): 鼓励发帖交流~ 2011-09-07 11:54:33
本帖仅楼主可见
赞一下
110楼2011-09-05 09:23:29
已阅   申请程序强帖   回复此楼   编辑   查看我的主页
相关版块跳转 我要订阅楼主 maomao1210 的主题更新
15811/16返回列表首页上一页91011121314下一页末页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 327求调剂 +3 Xxjc1107. 2026-04-06 3/150 2026-04-06 20:15 by lijunpoly
[考研] 复试调剂 +13 呼呼？~+123456 2026-04-05 13/650 2026-04-05 22:07 by 醉翁wl
[考研] 生物与医药086000调剂一志愿西北农林320分 +3 美美女士 2026-04-03 3/150 2026-04-05 21:55 by 学员8dgXkO
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4 cs0106 2026-04-04 4/200 2026-04-05 18:46 by imissbao
[考研] 材料调剂 +9 革微桂 2026-04-04 9/450 2026-04-05 08:27 by 544594351
[考研] 考研调剂 +6 15615482637 2026-04-04 6/300 2026-04-04 22:43 by yu221
[考研] 278求调剂 +3 依旧！ 2026-04-02 4/200 2026-04-04 20:27 by 蓝云思雨
[考研] 291求调剂 +4 迷蒙木木 2026-04-01 5/250 2026-04-04 15:59 by sihailian3
[考研] 土木304求调剂 +4 兔突突突， 2026-03-31 4/200 2026-04-04 13:34 by 1753564080
[考研] 400分求调剂 +3 尴尬且挠头 2026-04-04 3/150 2026-04-04 08:41 by jp9609
[考研] 268求调剂 +8 你好tg 2026-04-03 9/450 2026-04-04 05:08 by gswylq
[考研] 320调剂 +4 农业工程与信息� 2026-04-03 4/200 2026-04-03 21:40 by lbsjt
[考研] 机械专硕297 +3 Afksy 2026-04-03 3/150 2026-04-03 14:24 by 1753564080
[考研] 283求调剂 +3 jiouuu 2026-04-03 4/200 2026-04-03 13:28 by jiouuu
[考研] 085600，320分求调剂 +6 大馋小子 2026-04-02 6/300 2026-04-02 21:54 by dongzh2009
[考研] 农学考研求调剂 +3 dkdkxm 2026-04-01 3/150 2026-04-02 16:04 by wangjagri
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +8 双马尾痞老板2 2026-04-02 9/450 2026-04-02 14:45 by 5896
[考研] 314求调剂 +11 1xiaojun23 2026-03-31 12/600 2026-04-02 12:31 by 1xiaojun23
[考研] 340求调剂 +4 希望如此i 2026-03-31 4/200 2026-03-31 16:40 by 690616278
[考研] 福建理工大学材料学院先进合金团队招收考研调剂学生 +3 大华金商都 2026-03-30 4/200 2026-03-31 01:04 by 方英俊602
信息提示
关闭
请填处理意见
关闭