| 查看: 924 | 回复: 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 |
回复此楼» 猜你喜欢
情人节自我反思：在爱情中有过遗憾吗？
已经有4人回复
-
球磨粉体时遇到了大的问题，请指教！
已经有12人回复
-
江汉大学解明教授课题组招博士研究生/博士后
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请高手指点一下怎么确定一个符号矩阵中的等值元素
已经有6人回复
- 陌生的fortran语言，哪位可以指点一二？悬赏100金币
已经有5人回复
- 求助计算机高人，编写个Fortran程序
已经有13人回复
- 请教一个fortran小程序编译出错的问题，谢谢
已经有9人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】求高手关于Fortran数据读写的问题。【已完结】
已经有12人回复
- 【求助】求指点，在Fortran里面，怎样实现“数组维度可调”的数组?
已经有5人回复
- 【求助】求助个Fortran计算，应该不难，半天搞不定，请高手帮忙
已经有11人回复
» 抢金币啦！回帖就可以得到:
-
985教授征女友
+1/246
-
西湖大学拓扑光学、非厄米光学、太赫兹方向博士后招聘
+2/236
-
西湖大学2026年秋季入学物理学、光学、电子信息方向博士生有名额速来！！！
+2/224
-
陆军军医大学第二附属医院（新桥医院）冉茜课题组招聘科研人员
+1/79
-
龙凤Tai——写给恋人的第100封情书
+1/77
-
北京-89175-事业单位-诚征女友
+1/63
-
有没有人做过这种结构的顺式体向反式体的转化？
+1/33
-
澳大利亚麦考瑞大学(Macquarie University)国际博士硕士全额奖学金-计算机-26年中开学
+1/32
-
同济大学脑机智能团队脑机接口方向招生招聘
+1/27
-
上海交通大学-宁波东方理工大学联合培养博士生
+1/14
-
中科院深圳先进技术研究院招聘免疫及计算生物学方向博士后、科研助理或联培学生
+1/9
-
中科院深圳先进技术研究院招联培学生 -- 多中心大队列数据已就绪，助你快车道产出成果
+1/9
-
德国图宾根大学诚招全奖岗位制博士（地下流固化学反应耦合数值模拟方向）
+1/8
-
中北大学冯瑞教授*开山大弟子*招募
+1/6
-
墨尔本大学(QS13)急招CSC博士(补齐全奖)/访问学者/博士后 (材料/生物医学/器官芯片等)
+1/5
-
德国图宾根大学诚招全奖岗位制博士（地下流固化学反应耦合数值模拟方向）
+1/5
-
墨尔本大学(QS13)急招CSC博士(补齐全奖)/访问学者/博士后(生物医学材料/器官芯片)
+1/4
-
26储能博士申请自荐
+1/1
-
澳科大招收2026年秋季药物递送/生物材料方向全奖博士研究生（3月5日18:00截止）
+1/1
-
苏州大学国家级青年人才团队2026年博士招生（有机光电功能材料方向）
+1/1
5楼2011-04-04 23:08:32
★ ★ ★ ★
田山东(金币+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