各位fortran高手, 刚刚我花了好多时间?在简查我的程式 可是始终不知道问题所在 某一部份的函数计算 执行後它一直说我INF 可是我check过整个函数内容 计算部份并不会造成INF 所以难道有其它原因造成INF吗 可是我简察不出来 这麽长的程式要发问 真的觉得很麻烦和抱歉 我下面有注明我的问题 能不能有高手能帮我检察出问题呢 感激不尽>< program paper1 double precision L,ione,bone,Cone,r,ra,u,e,ea,eb,ec,Tzero double precision c,k,Abig,v,mrate,q,Hone,paa,pbb,Tone,Ttwo double precision m,G,pi,b,h,lightv,boltzman,Tr,z,calculusansA double precision pa,pb,Asmall,Done,Dtwo,gamme,zete,number double precision gnumber,znumber,p,eav,rbv,alpha,beta double precision Tonev,Ttwov,Trv,n,Msun,SBConstant,Tcombine integer o ******************************************************************** ******************************Constant Number Part************* pi = 3.1415926D0 G = 6.6 h = 6.6260693D-27 k = 1.3806505D-26 pa = 0D0 pb = 0D0 paa = 0D0 pbb = 0D0 c = 299792458D2 Tzero = 500D0 r = 1D0 SBConstant = 5.6704004D-7 Msun = 1.988435D33 Rstar = 3*(10**11)*((Mrate/1D-5/Msun)**0.33333333) L = G*M*Mrate/365/24/60/60/Rstar Asmall = 7.5657D-15 Tcombine = ((alpha*Tonev)**4+(beta*Ttwov)**4)**0.25 Cone = Mrate/365/24/60/60/4/pi/sqrt(2*G*M) write(*,*) 'number?,gnumber?,znumber?,Mrate,M,n,A,alpha,beta,r' read(*,*) number , gnumber , znumber ,Mrate , M ,n,Abig,alpha,beta -,r write(*,*) Done(Abig,Asmall,Cone,L,c,k,h,n,gnumber,znumber) <-----执行後老是说INF double precision function Done(Abig,Asmall,Cone,L,c,k,h,n, <-----就是这个函数 -gnumber,znumber) double precision Abig,Asmall,Cone,L,c,k,h,n,gnumber,znumber,pi Done = (3*6*1.08232323*Abig* <----但我check过计算式 -Cone*L/10/pi/2/1.2020564/Asmall/c* 没问题啊>< -(k/h)**1)**0.3333333 end double precision function Dtwo(r,n,L,Asmall,c,Tzero,gnumber, -znumber) double precision r,n,L,Asmall,c,Tzero,gnumber,znumber Dtwo = (6*1.08232323*Tzero**1*L/4/pi/ -24/1.03692776/Asmall/c)**0.2 end --

※ 发信站: 批踢踢实业坊(ptt.cc) ◆ From: 140.109.177.125 ※ 编辑: pianoboy 来自: 140.109.177.125 (12/18 00:02)