我的程序 输入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
感谢回复！

 

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

复兴汉服，还我中华衣冠，复我中华礼仪
http://www.tianhan.com.cn/bbs/index.asp
<img>http://memory.116.com.cn/uploadfiles/ttmovie/2005222154103.jpg</img>