24小时热门版块排行榜

返回列表

田山东

捐助贵宾 (著名写手)

应助: 5 (幼儿园)
金币: 77
帖子: 1577
在线: 138.1小时
虫号: 800307

[交流] 【求助】一个求积分的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

回复此楼

» 猜你喜欢

拟解决的关键科学问题还要不要写已经有11人回复
救命帖已经有5人回复
限项规定已经有5人回复
为什么nbs上溴没有产物点出现呢已经有9人回复
招博士已经有3人回复
存款400万可以在学校里躺平吗已经有35人回复
最失望的一年已经有18人回复
求推荐博导已经有4人回复
求推荐英文EI期刊已经有5人回复
疑惑？已经有5人回复

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

请高手指点一下怎么确定一个符号矩阵中的等值元素已经有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的回帖

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的回帖