[交流] 【求助】一个求积分的fortran程序，请高手指点（还是在物理区问好些）

我算了一个很简单的求积分的程序，得到的结果为什么是0呢？用imsl下的数据库函数QDAGS (F, a, b, RES, ERRA,ERRREL, ERREST)求的。高手帮忙指点下吧。不论函数f(x)取什么形式，得到结果总是0.
program jifenQDAGS
   use imsl
   implicit none
   real a,b,erra,errest,f,res,errrel
   external f
   ERRA=1.0e-3
   ERRREL=1.0e-2
   a=0.0
   b=3.0
CALL QDAGS (F, a, b, RES, ERRA,ERRREL, ERREST)
WRITE (*,*) RES, ERREST
END
REAL FUNCTION F (X)
REAL X
F = x*x
RETURN
END

回复此楼

» 猜你喜欢

上海工程技术大学【激光智能制造】课题组招收硕士已经有6人回复
带资进组求博导收留已经有11人回复
自荐读博已经有5人回复
求个博导看看已经有16人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

请高手指点一下怎么确定一个符号矩阵中的等值元素已经有6人回复
陌生的fortran语言，哪位可以指点一二？悬赏100金币已经有5人回复
求助计算机高人，编写个Fortran程序已经有13人回复
请教一个fortran小程序编译出错的问题，谢谢已经有9人回复
写了一个fortran90的小程序，编译通不过，请大侠帮忙已经有59人回复
【求助】求高手关于Fortran数据读写的问题。【已完结】已经有12人回复
【求助】求指点，在Fortran里面，怎样实现“数组维度可调”的数组? 已经有5人回复
【求助】求助个Fortran计算，应该不难，半天搞不定，请高手帮忙已经有11人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-03-13 18:29:34

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhuwei04

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 2536.5
帖子: 372
在线: 91.6小时
虫号: 162510

★ ★ ★ ★
田山东(金币+1):谢谢参与
lpszk(金币+3): 2011-03-13 22:32:33
田山东(金币+15): 多谢您了，我太粗心了。 2011-03-14 08:47:42

引用回帖:

Originally posted by 田山东 at 2011-03-13 18:29:34:
我算了一个很简单的求积分的程序，得到的结果为什么是0呢？用imsl下的数据库函数QDAGS (F, a, b, RES, ERRA,ERRREL, ERREST)求的。高手帮忙指点下吧。不论函数f(x)取什么形式，得到结果总是0.
program jifenQDAG ...

你把RES的位置放错了，修改如下：
program jifenQDAGS
   use imsl
      implicit none
      real a,b,erra,errest,f,res,errrel
      external f
   ERRA=1.0e-3
   ERRREL=1.0e-2
      a=0.0
      b=3.0
CALL QDAGS (F, a, b, ERRREL, ERRA,RES, ERREST)
WRITE (*,*) RES, ERREST
      END
REAL FUNCTION F (X)
      REAL X
      F = x*x
      RETURN
END

赞一下(2人)

回复此楼

3楼2011-03-13 20:19:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 5 个回答

zhuwei04

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 2536.5
帖子: 372
在线: 91.6小时
虫号: 162510

结果

★ ★ ★ ★ ★ ★ ★
田山东(金币+2): 再次感谢 2011-03-14 08:47:54
小木虫:)(金币+7): 3q 2011-04-05 07:20:30

引用回帖:

Originally posted by zhuwei04 at 2011-03-13 20:19:11:
你把RES的位置放错了，修改如下：
program jifenQDAGS
   use imsl
      implicit none
      real a,b,erra,errest,f,res,errrel
      external f
   ERRA=1.0e-3
   ERRREL=1.0e-2
...

附上结果图

赞一下(1人)

回复此楼

4楼2011-03-13 20:21:22

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖