DCUHRE: 多维数值积分 1.0 文件大小：41 授权方式：免费版 软件类别：Fortran 软件语言：简体中文 运行环境：Win9x/Me/NT/2000 开 发 商：Jarle Berntsen  E-mail 相关链接：开发商主页  程序演示 更新时间：2009/11/21 21:38:17 解压密码： 推荐等级： 下载次数：

下载地址

下载地址1

相关软件

当前没有记录！

软件简介

DCUHRE

该程序是一个多维数值积分的Fortran程序，我觉得挺好使的，这块儿共享给大家。里面我加了两个测试函数，一个是计算半径为R的圆球体积；另一个是计算一个归一化的高斯分布。

This is a Fortran code to solve multi-dimensional numerical integration

简单解释一下用法：

PARAMETER (NDIM = 3, NW = 5000, NF = 2)

NDIM: 积分的维数  
NW: 是算法里面要用到的一个数组的长度，要足够大
NF :  要计算的积分个数 

 A(1)=0
 B(1)=12
 A(2)=0
 B(2)=PI
 A(3)=0
 B(3)=2*PI
数组A()和B()分别给出被积空间的下限和上限
    

下面就是用户自己定义自己的被积函数了，这块儿给了球坐标系下两个测试函数F(1)和F(2)
第一个积分结果应该是4*Pi/3*12^3=7238.229
第二个积分结果应该是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
!===========================================================

 

网友评论 发表评论