[請益] 與積分有關係的非線性方程

看板Fortran作者trashken (上個暱稱)時間13年前 (2011/09/08 16:52)推噓5(5推 0噓 11→)

留言16則, 3人參與討論串1/1

程式原始碼: PROGRAM guess use IMSL IMPLICIT none external FCN real, parameter :: ERRREL=0.0001 integer, parameter :: N=1 integer, parameter :: ITMAX=100 real :: XGUESS(N)=(/0.4/) real X(N), FNORM CALL NEQNF(FCN, ERRREL, N, ITMAX, XGUESS, X, FNORM) write(*,*) X stop end subroutine FCN (XA, F, N) implicit none integer N real, target ::XA(N) real F(N) real*8 ::El,kT,ul,ur,Tl,Tr,pi,fl,fr,Ea,T,VL,Y REAL,POINTER ::N0 N0=>XA(1) El=450 kT=0 ul=60 ur=60 VL=60 Tl=1 Tr=1 T=Tl+Tr pi=3.14159 fl=1/(exp((Ea-ul)/kT)+1) fr=1/(exp((Ea-ur)/kT)+1) do Ea=0.d0,300.d0,0.1d0 Y=Y+(-(fl*Tl+fr*Tr))/(pi*(Tl+Tr))*((-(1-N0)*(T/2))/((Ea-El)**2+(T/2)**2)-(N0*T /2)/((Ea-El-VL)**2+(T/2)**2)) end do F(1)=Y-N0 RETURN END Subroutine ============================================================================= 可不可以請問一下版上各位大大我目前在解一個與積分有關的方程式也可以說是一個耦合方程式要積分的式子是 (-(fl*Tl+fr*Tr))/(pi*(Tl+Tr))*((-(1-N0)*(T/2))/((Ea-El)**2+(T/2)**2)-(N0*T /2)/((Ea-El-VL)**2+(T/2)**2)) 而積分的變數是Ea 而要解的未知數是N0 而N0也等於這一長串要積分的式子也就是說N0=(-(fl*Tl+fr*Tr))/(pi*(Tl+Tr))*((-(1-N0)*(T/2))/((Ea-El)**2+(T/2)**2 )-(N0*T/2)/((Ea-El-VL)**2+(T/2)**2)) 而我使用函式庫的NEQNF功能而使用一個簡單的DO迴圈來積那一長串的式子再令這一長串的式子為Y 而後在DO迴圈外使用"F(1)=Y-N0"這個式子來求解這個未知數N0 程式可以跑但問題是跑出來的值跟我要模擬的圖相比較之下都怪怪的在這裡可否請教各位如果欲求解的未知數出現在積分式中那要怎麼樣來寫方程式才會比較恰當我這樣用DO迴圈直接求解的方式好像錯了這個問題已經卡住小弟兩三個月了希望有大大可以替小弟解惑 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.115.72.52 ※ 編輯: trashken 來自: 140.115.72.52 (09/08 16:54)