各位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)