以文本方式查看主题 - 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 很容易溢出,注意其取值范围 |