下載地址
相關(guān)軟件
當前沒(méi)有記錄!
軟件簡(jiǎn)介
DCUHRE
該程序是一個(gè)多維數值積分的Fortran程序,我覺(jué)得挺好使的,這塊兒共享給大家。里面我加了兩個(gè)測試函數,一個(gè)是計算半徑為R的圓球體積;另一個(gè)是計算一個(gè)歸一化的高斯分布。
This is a Fortran code to solve multi-dimensional numerical integration
簡(jiǎn)單解釋一下用法:
PARAMETER (NDIM = 3, NW = 5000, NF = 2)
NDIM: 積分的維數
NW: 是算法里面要用到的一個(gè)數組的長(cháng)度,要足夠大
NF : 要計算的積分個(gè)數
A(1)=0
B(1)=12
A(2)=0
B(2)=PI
A(3)=0
B(3)=2*PI
數組A()和B()分別給出被積空間的下限和上限
下面就是用戶(hù)自己定義自己的被積函數了,這塊兒給了球坐標系下兩個(gè)測試函數F(1)和F(2)
第一個(gè)積分結果應該是4*Pi/3*12^3=7238.229
第二個(gè)積分結果應該是1
==== user's function =============================================
SUBROUTINE FTEST(NDIM, X, NFUN, F)
INTEGER N, NDIM, NFUN
DOUBLE PRECISION X(NDIM), F(NFUN), SUM
PI=3.1415926
w2=1
F(1) = X(1)*X(1)*SIN(X(2))
F(2) = X(1)*X(1)*SIN(X(2)) * EXP(-X(1)**2/(2*w2)) / (2*Pi*w2)**1.5
END
!===========================================================