| 查看: 692 | 回复: 3 | ||
luohubin木虫 (正式写手)
|
[求助]
关于一个积分算法的问题（有代码）
|
|
请教大家一个散点积分的问题。我将fortran代码贴在下面。 说明一下参数：f是一维矩阵，存放n个散点。散点间等距，为d。fi是积分值。 问题是，我看不出这个积分是基于什么方法的，在数值分析书上没找到对应的。希望大家给我点建议。谢谢！ SUBROUTINE gensim(f,d,n,fi) IMPLICIT NONE REAL(KIND=8), DIMENSION(*) :: f REAL(KIND=8) :: d, fi, ff, corr INTEGER :: n, n1, j, n2 C IF(n.GE.3) THEN IF(MOD(n,2).EQ.0) THEN n1=n-1 ff=-f(n-2)+8.d0*f(n-1)+5.d0*f(n) corr=ff*d/12.d0 ELSE n1=n corr=0.d0 ENDIF n2=n1-1 ff=0.d0 DO 20 j=2,n2,2 20 ff=ff+(f(j-1)+4.d0*f(j)+f(j+1)) fi=d*ff/3.d0+corr ELSEIF(n.EQ.2) THEN ff=5.d0*f(1)+8.d0*f(2)-f(3) fi=d*ff/12.d0 ELSEIF(n.LE.1) THEN fi=0.d0 ENDIF C RETURN END |
回复此楼» 猜你喜欢
估计是周四
已经有3人回复
-
基金啊基金
已经有3人回复
-
科研孤儿太难了
已经有15人回复
-
让我中一个面上吧！
已经有10人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 一个算法多做个领域的问题，算不算一稿多投
已经有17人回复
- vasp算法的一些基础问题
已经有10人回复
- 我编的Simpson积分法fortran程序给不出结果，大侠们看看哪里出了问题？
已经有4人回复
- 同一个问题，不同算法的两篇论文，导师说这两篇投稿时间不能太近？
已经有26人回复
- 关于MATLAB和OpenCV对相同算法的运算速度对比问题
已经有20人回复
- 请教关于遗传算法编程的问题（俺是小白）
已经有4人回复
- 关于免疫算法的几个问题。（搜索无重复）
已经有8人回复
- 【求助】DLS中数据算法问题
已经有16人回复
- 关于算法比较的问题，请大家帮忙
已经有7人回复
- 【求助】k-medoids聚类算法关于中心替换的问题
已经有5人回复
- 【讨论】预调件共轭梯度法（PCG）
已经有50人回复
- 【求助】求积分程序中被积函数问题
已经有15人回复
luohubin
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1562.8
- 散金: 50
- 红花: 24
- 帖子: 604
- 在线: 132.9小时
- 虫号: 429401
- 注册: 2007-08-04
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2012-03-31 10:43:29
luohubin
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1562.8
- 散金: 50
- 红花: 24
- 帖子: 604
- 在线: 132.9小时
- 虫号: 429401
- 注册: 2007-08-04
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2012-04-13 22:04:43
luohubin
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1562.8
- 散金: 50
- 红花: 24
- 帖子: 604
- 在线: 132.9小时
- 虫号: 429401
- 注册: 2007-08-04
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2012-04-13 22:09:58
40