| 查看: 1051 | 回复: 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 |
回复此楼» 猜你喜欢
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有27人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有4人回复
-
UV压敏胶开发
已经有6人回复
-
有没有H口的？有收到消息的吗？
已经有3人回复
-
十年后又回来了，论文投稿求助
已经有3人回复
-
纯娱乐，不喜欢勿喷
已经有10人回复
-
面上提前没消息，有中的吗
已经有18人回复
-
娱乐
已经有4人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有16人回复
-
一个有机合成实验室都需要哪些设备？
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请高手指点一下怎么确定一个符号矩阵中的等值元素
已经有6人回复
- 陌生的fortran语言，哪位可以指点一二？悬赏100金币
已经有5人回复
- 求助计算机高人，编写个Fortran程序
已经有13人回复
- 请教一个fortran小程序编译出错的问题，谢谢
已经有9人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】求高手关于Fortran数据读写的问题。【已完结】
已经有12人回复
- 【求助】求指点，在Fortran里面，怎样实现“数组维度可调”的数组?
已经有5人回复
- 【求助】求助个Fortran计算，应该不难，半天搞不定，请高手帮忙
已经有11人回复
» 抢金币啦！回帖就可以得到:
-
山东征女友，坐标济南
+1/175
-
2-氯-1,3,2-二氧磷杂环戊烷（CAS：822-39-9）
+1/98
-
威格科技张德龙：二十年深耕，铸就国产惰性气氛设备新高度
+1/82
-
具身智能机器人公司招聘机械结构岗位（工程师、专家、负责人均可）薪酬最高150万
+1/81
-
南京大学医工交叉方向招聘脂质纳米材料与RNA递送方向博士后1-2名
+1/77
-
药化高级研究员
+1/77
-
小木虫是我和国基唯一的人脉，现在都感觉快断了
+1/58
-
有机合成高级研究员
+1/27
-
双一流高校招收博士（现代技术管理/环境技术经济管理）及硕士研究生（长期有效）
+1/22
-
天津工业大学招26师资博后/27博士/科研助理
+1/17
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/16
-
北理工集成电路杰青团队 | 诚招科助理
+1/15
-
苏州大学李刚教授课题组博士后招聘公告
+1/12
-
福建师范大学海峡柔性电子学院-2027年9月入学博士-考核制-电催、能源
+1/10
-
多本Scopus收录国际英文期刊征稿，全免 APC！含医学期刊
+1/4
-
北京师范大学（珠海校区）赵亚然副教授招收2027年电催化推免硕士/博士研究生
+1/4
-
2027年硕士推免生 招生--哈尔滨工业大学王磊
+1/3
-
国家自然科学基金申报者发榜前心理状态调查，麻烦帮忙填一下哇！
+1/3
-
组委会诚挚滴邀请各界学者踊跃投稿参会
+1/2
-
致每一位等待国自然结果的科研人：我们想听听你的心声
+1/1
★ ★ ★ ★
田山东(金币+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
