| 查看: 1074 | 回复: 4 | |||
[交流]
【求助】一个求积分的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 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请高手指点一下怎么确定一个符号矩阵中的等值元素
已经有6人回复
- 陌生的fortran语言，哪位可以指点一二？悬赏100金币
已经有5人回复
- 求助计算机高人，编写个Fortran程序
已经有13人回复
- 请教一个fortran小程序编译出错的问题，谢谢
已经有9人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】求高手关于Fortran数据读写的问题。【已完结】
已经有12人回复
- 【求助】求指点，在Fortran里面，怎样实现“数组维度可调”的数组?
已经有5人回复
- 【求助】求助个Fortran计算，应该不难，半天搞不定，请高手帮忙
已经有11人回复
★ ★ ★ ★
田山东(金币+1):谢谢参与
lpszk(金币+3): 2011-03-13 22:32:33
田山东(金币+15): 多谢您了，我太粗心了。 2011-03-14 08:47:42
田山东(金币+1):谢谢参与
lpszk(金币+3): 2011-03-13 22:32:33
田山东(金币+15): 多谢您了，我太粗心了。 2011-03-14 08:47:42
|
你把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 |
3楼2011-03-13 20:19:11
4楼2011-03-13 20:21:22
5楼2011-04-04 23:08:32
简单回复
2011-03-13 19:03
回复
田山东(金币+1):谢谢参与
田山东(金币+1): 这样可不好啊！帮我看下是怎么回事吧 2011-03-13 19:09:15
小木虫:)(金币-1): 禁止灌水 2011-04-05 07:20:00
田山东(金币+1): 这样可不好啊！帮我看下是怎么回事吧 2011-03-13 19:09:15
小木虫:)(金币-1): 禁止灌水 2011-04-05 07:20:00
