应《网络安全法》要求，自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下一页末页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 土木水利专硕276分求调剂 +3 我想上学！！6 2026-04-05 6/300 2026-04-06 18:51 by 我想上学！！6
[考研] 一志愿安徽某211 0703化学总分339求调剂 +7 晚风不晚 2026-04-04 7/350 2026-04-06 14:06 by houyaoxu
[考研] 材料工程302分求调剂 +9 zyx上岸！ 2026-04-04 9/450 2026-04-05 22:08 by 醉翁wl
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-04-04 5/250 2026-04-05 14:06 by imissbao
[考研] 一志愿西北农林畜牧专硕336分求调剂 +3 5ourr 2026-04-03 3/150 2026-04-05 10:40 by JOKER0401
[考研] 271分求调剂学校 +12 zph158488！ 2026-04-02 13/650 2026-04-05 10:13 by lqwchd
[考研] 男生，一志愿沪9生物学071000，初试308求调剂 +3 刘墨墨 2026-04-04 3/150 2026-04-05 08:26 by barlinike
[考研] 413求调剂 +4 柯某某 2026-03-31 4/200 2026-04-04 22:18 by 学员6BFVa3
[考研] 277求调剂 +4 12A3 2026-04-02 5/250 2026-04-04 20:28 by 蓝云思雨
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 14:46 by 湘农储能材料
[考研] 309求调剂 +6 刘刘刘1231 2026-04-02 7/350 2026-04-04 13:41 by liucky
[考研] 求调剂，一志愿北京中医药大学 +3 小小达不溜 2026-04-02 3/150 2026-04-03 22:55 by 冲矢昴星团
[考研] 285求调剂 +7 AZMK 2026-04-02 9/450 2026-04-03 11:12 by wanwan00
[考研] 274求调剂 +10 薛定谔的虎。 2026-04-01 10/500 2026-04-03 10:13 by tianyyysss
[考研] 302求调剂 +9 zyx上岸！ 2026-04-02 9/450 2026-04-02 23:07 by 马儿快快地跑
[考研] 260求调剂 +6 朱芷琳 2026-04-02 6/300 2026-04-02 20:27 by 6781022
[考研] 314求调剂 +11 1xiaojun23 2026-03-31 12/600 2026-04-02 12:31 by 1xiaojun23
[考研] 一志愿346上海大学生物学 +3 上海大学346调剂 2026-04-01 3/150 2026-04-02 08:36 by w虫虫123
[考研] 307分求调剂 +14 (o~o) 2026-03-31 15/750 2026-04-01 20:43 by longlotian
[考研] 一志愿浙江大学工科动力工程370,数一121,专业课135，现在能去哪里 +3 080700调剂 2026-03-30 4/200 2026-03-31 12:00 by KLMY666
信息提示
关闭
请填处理意见
关闭