版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

totoymy

铜虫 (小有名气)

应助: 1 (幼儿园)
金币: 27.5
红花: 2
帖子: 63
在线: 16.3小时
虫号: 1326694
注册: 2011-06-19
性别: MM
专业: 理论和计算化学

[求助] 内循环问题，求大虾帮忙

MODULE PARA
IMPLICIT NONE

INTEGER*4, PARAMETER:: Ndim=12

REAL,DIMENSION(Ndim,5):: pp

INTEGER*4::NOUT

END MODULE PARA

PROGRAM CM
USE MSIMSL
USE PARA
IMPLICIT NONE

INTEGER:: i,j,k,m
INTEGER::t
real::x,y,z
REAL*8::r
REAL::s

OPEN (11, file='input.txt')

READ(11,*) ((pp(i,j),j=1,4),i=1,Ndim)
close(11)

t=1
DO while (t<=100)
call random_seed()
call random_number(r)

k=1
do while(k<=Ndim)
s=s+pp(k,4)
m=s-pp(k,4)

if (m x=x+pp(k,1)
y=y+pp(k,2)
z=z+pp(k,3)
open(12,file="result.txt",access='append')
write(12,*) x,y,z,t,k,r
close(12)
end if
k=k+1
enddo

t=t+1
if (t>100) exit
enddo

ENDPROGRAM
大意就是生成100次随机数，并归属每一个随机数的范围，若属于某特定范围，则对其坐标进行累计加和。。求高手啊高手~~数值见附件，

回复此楼

» 猜你喜欢

Suzuki偶联中硼酸酯原料掉下的硼酸酯是否可以自身偶联，生成B2Pin2杂质已经有3人回复
2026年面上项目中了，2A+B, 会评顺利通过已经有11人回复
大龄残疾硕士的一点执念已经有19人回复
今年E04面上已经有19人回复
刑法学论文投稿求助已经有5人回复
Journal of Environmental Chemical Engineering 已经有3人回复
最后一年，祈求好运已经有3人回复
风电环氧领域已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

关于数学符号的问题？请教大虾已经有16人回复
杂志选择问题，跪求大虾帮助已经有7人回复
一个脱苄问题，向各位大虾求助已经有33人回复
急求一个有机合成问题，取代反应，请各位大虾帮忙！虫虫没有金币了，只散30个金币…… 已经有11人回复
请问下BET大虾有关滞后环的问题已经有5人回复
辰华CHI660C测试扣式锂电池循环伏安电流总是显示很大，安级，请教大虾已经有6人回复
【求助】关于未知树状结构存储的问题？请大虾指教！已经有11人回复
【求助】各位大虾们，请教酰化试剂溶解性质的问题，谢谢啦已经有12人回复
【求助】请教各位大虾，有关环氧胶黏剂脱粘的问题已经有9人回复
【求助】看文献中遇到一个问题，请教各位大虾! 已经有10人回复
急求环境专业的翻译大虾帮忙翻译下中文翻译成英文已经有2人回复
【求助】小妹求助各位大虾红外分析相对强度问题已经有15人回复
【求助】求教橡胶表面处理问题，大虾们帮帮忙啊！已经有6人回复
【求助】过柱子的问题，请各位大虾帮帮忙，谢谢了！已经有22人回复
【求助】培养单晶的大虾来帮帮忙啊已经有21人回复

1楼 2011-10-28 10:57:29

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

totoymy

铜虫 (小有名气)

应助: 1 (幼儿园)
金币: 27.5
红花: 2
帖子: 63
在线: 16.3小时
虫号: 1326694
注册: 2011-06-19
性别: MM
专业: 理论和计算化学

引用回帖:

1楼: Originally posted by totoymy at 2011-10-28 10:57:29:
MODULE PARA
IMPLICIT NONE

INTEGER*4, PARAMETER:: Ndim=12

REAL,DIMENSION(Ndim,5):: pp

INTEGER*4::NOUT

END MODULE PARA

PROGRAM CM
USE MSIMSL
USE PARA
IMPLICIT NONE

INTEGER: ...

帮我说说问题也好。。。应该一共出100组数。。。结果总是不对啊~

赞一下

回复此楼

3楼2011-10-30 12:53:52

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 4 个回答

snoopyzhao

至尊木虫 (职业作家)

程序强帖: 16
应助: 157 (高中生)
贵宾: 0.02
金币: 18844.7
红花: 29
帖子: 3803
在线: 1422.4小时
虫号: 183750
注册: 2006-02-13
专业: 污染生态化学

ben_ladeng: 专家考核存档 2011-10-30 07:11:40

没有相应的库，不应助了……

赞一下(1人)

回复此楼

2楼2011-10-28 16:41:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

程序强帖: 16
应助: 157 (高中生)
贵宾: 0.02
金币: 18844.7
红花: 29
帖子: 3803
在线: 1422.4小时
虫号: 183750
注册: 2006-02-13
专业: 污染生态化学

【答案】应助回帖

★ ★
jjdg(金币+2): 辛苦了 2011-10-30 13:30:00
totoymy(金币+20): 貌似解决了~谢谢大虾，我是在别的文件基础上编的，所以那个库是原来用的。。没有用 2011-10-31 15:11:50

引用回帖:

3楼: Originally posted by totoymy at 2011-10-30 12:53:52:
帮我说说问题也好。。。应该一共出100组数。。。结果总是不对啊~

看了一下，你根本没有必要引用 MSIMSL 这个库。不知道你说的分组是啥个意思……

我改了一下，可以出来结果，但不见得是你想要的：

CODE:

MODULE PARA

IMPLICIT NONE

INTEGER*4, PARAMETER:: Ndim=12   

REAL,DIMENSION(Ndim,5):: pp

INTEGER*4::NOUT

END MODULE PARA

PROGRAM CM

!USE MSIMSL

USE PARA

IMPLICIT NONE

INTEGER:: i,j,k

INTEGER::t

real::x,y,z

REAL*8::r

REAL::s,m

OPEN (11, file='input.txt')

READ(11,*) ((pp(i,j),j=1,4),i=1,Ndim)

close(11)

   call random_seed()

t=1

DO 

   call random_number(r)

!  write (*,*) r

   k=1

   s = 0.0

   do while(k<=Ndim)

      s=s+pp(k,4)

      m=s-pp(k,4)

!     write(*,*) s, m

      if (m
          x=x+pp(k,1)

          y=y+pp(k,2)

          z=z+pp(k,3)

          open(12,file="result.txt",access='append')

          write(12,*) x,y,z,t,k,r

          close(12)  

          t = t + 1

      end if

      k=k+1

   enddo

   if (t>100) exit

enddo

end program cm

赞一下(1人)

回复此楼

4楼2011-10-30 13:18:46

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 4 个回答