应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 关于一个积分算法的问题（有代码）
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）
51/1返回列表
查看: 515  |  回复: 3
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
消愁无计对金樽
1楼2012-03-28 20:52:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luohubin

木虫 (正式写手)
用自己帖子测试一下新学的公式编辑方法
[tex]f(x)=\int_{-\infinity}^xe^x[/tex]
赞一下 回复此楼
消愁无计对金樽
3楼2012-04-13 22:04:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

luohubin

木虫 (正式写手)
这版块这么冷清，就这么几行代码也没人帮个忙吗？
赞一下 回复此楼
消愁无计对金樽
2楼2012-03-31 10:43:29
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luohubin

木虫 (正式写手)
再来一下
[tex]\Large f(x)=\int_{-\infty}^x e^{-t^2}[/tex]
赞一下 回复此楼
消愁无计对金樽
4楼2012-04-13 22:09:58
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭