[問題] 對不起,要麻煩Fortran高手幫我解決這個頭痛的問題><
各位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)
推
12/18 00:57, , 1F
12/18 00:57, 1F
→
12/18 00:58, , 2F
12/18 00:58, 2F
推
12/18 01:00, , 3F
12/18 01:00, 3F
→
12/18 01:02, , 4F
12/18 01:02, 4F