| 查看: 542 | 回复: 3 | ||
[求助]
求高手帮忙解答这个程序问题 已有1人参与
|
没有提示warning,也没有提示error,但是运行的时候，却说触发了一个断点，让选择继续或中断，选择继续后，根据得不到想要的数据。这到底怎么回事呢？求高手帮忙看看，不胜感激 subroutine bohanshu(n,bhs,energy) IMPLICIT NONE !+++++++++++++++++++++++++++++++++++++++++++++++声明 integer,parameter::num=10 !num1代表n，n表示第n个态 integer,parameter::num1=100 !num2代表x,即坐标 real, parameter::pi=3.14 real::bhs(num,num1) !放波函数的数组 real::energy(num) !存放能量本征值得数组 integer::n real::x,a real ::xx real ::dx integer:: j !++++++++++++++++++++++++++++++++++++++++++++++++赋值 xx=2.0 !空间宽度 dx=xx/num1 a=1.50 !势阱宽度 !++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++运算 do n=1,num !++++++++++++++++++++++++++++n循环 energy(n)=pi**2*(n**2)/(2*a**2) do j=1,num1 !+++++++++++++++++++++++++++x循环 x=dx*j bhs(n,j)=sqrt(2/a)*sin(n*pi*x/a) end do end do return end !******************************************************求微扰矩阵Hnk以及exp(energy(n)-energy(k))*t的子程序 program bobo IMPLICIT NONE !+++++++++++++++++++++++++++++++++++++++++++++++声明 integer,parameter::num=10 !在该程序中选取的最高能量本征态对应的量子数 integer,parameter::num1=100 !位置的份数 integer,parameter::num2=200 !时间的份数 real ,parameter::pi=3.14 real::H !哈密顿量中的微扰部分 real::bhs(num,num1) !放波函数的数组 real::WRJZnk(num) !存放微扰矩阵的数组 real::energy(num) !存放能量本征值的数组 complex::Wnk(num,num2) !用来存放exp(energy(n)-energy(k))*t的数组 complex::b !虚数单位i real::x !位置变量 real::xx !能到达的最大位置 real ::dx !表示位置的变化 integer::k !第k个态的量子数 integer n !是第n个态的量子数，同时也是能量的循环指标 integer ::j !位置循环指标 integer::i !时间循环指标 real::t !时间 real::tx !微扰作用在系统上的截止时间 real::dt !数值模拟的时间间隔 !++++++++++++++++++++++++++++++++++++++++++++++++赋值 tx=10.0 xx=2.0 H=1.105 !微扰 !k=4 !代表初始状态为第4个能量本征态 dx=xx/num1 !位置的数值间隔 dt=tx/num2 !时间的数值间隔 b=(0,1) !++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++运算 open(unit=60,file='WRJZnk.txt') open(unit=70,file='Wnk.txt') call bohanshu(n,bhs,energy) do n=1,num !#################################################################################能量循环 do j=1,num1 !#################################################################################位置循环 WRJZnk(n)=WRJZnk(n)+bhs(n,j)*H*bhs(k,j)*dx !###########################################################计算出了Hnk end do write(60,*)WRJZnk(n) end do do n=1,num !################################################################################能量循环 do i=1,num2 !################################################################################时间循环 t=i*dt Wnk(n,i)=exp(b*(energy(n)-energy(k))*t) write(70,*)Wnk(n,i) end do end do end program bobo |
回复此楼» 猜你喜欢
材料考研调剂
已经有3人回复
-
材料调剂
已经有12人回复
-
英一数一408，总分284，二战真诚求调剂
已经有14人回复
-
085410 一志愿211 22408分数359求调剂
已经有4人回复
-
271求调剂
已经有19人回复
-
385分 生物学（071000）求调剂
已经有3人回复
-
一志愿安徽大学计算机科学与技术学硕，331分求调剂
已经有3人回复
-
318求调剂，计算材料方向
已经有8人回复
-
291求调剂
已经有25人回复
-
一志愿北京科技大学085601材料工程英一数二初试总分335求调剂
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
wxws.2008
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 5015.5
- 散金: 256
- 红花: 9
- 帖子: 354
- 在线: 147.9小时
- 虫号: 639701
- 注册: 2008-10-29
- 性别: GG
- 专业: 荒漠化与水土保持
2楼2015-12-29 20:41:09
3楼2016-01-01 20:47:19
4楼2016-01-09 09:26:44
