以文本方式查看主题

-  Fortran中文网  (http://fortran.cn/bbs/index.asp)
--  Fortran语言使用经验交流  (http://fortran.cn/bbs/list.asp?boardid=2)
----  溢出怎么改  (http://fortran.cn/bbs/dispbbs.asp?boardid=2&id=433)

--  作者:zhang_018
--  发布时间:2006/10/23 20:45:51

--  溢出怎么改
我的程序 输入U=0.0,t=1.0e-23显示溢出错误,该怎么改

implicit real(8)(B,d,Z,m)

real,parameter:: kB=1.38e-23,miuB=9.274e-24,  kBT=1.0E-23

read* ,U,t

do B=0.1,2.0,0.1      

delta=sqrt(U**2/4.0+4.0*t**2)

ZBG=1.0+exp(-U/kBT)+exp(-2*B/kBT)+exp(2*B/kBT)+exp(-(U/2.0+delta)/kBT)+exp(-(U/2.0-delta)/kBT)

mBG_miuB=2.0*(exp(2.0*B/kBT)-exp(-2.0*B/kBT))/ZBGM

write(10,*)mBG_miuB

enddo

end

感谢回复!


--  作者:齐东野人
--  发布时间:2006/10/24 3:35:28

--  
exp 很容易溢出,注意其取值范围

京ICP备05056801号