| 查看: 927 | 回复: 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 O5 51O 54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有10人回复
-
球磨粉体时遇到了大的问题，请指教！
已经有13人回复
-
情人节自我反思：在爱情中有过遗憾吗？
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
江汉大学解明教授课题组招博士研究生/博士后
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请高手指点一下怎么确定一个符号矩阵中的等值元素
已经有6人回复
- 陌生的fortran语言，哪位可以指点一二？悬赏100金币
已经有5人回复
- 求助计算机高人，编写个Fortran程序
已经有13人回复
- 请教一个fortran小程序编译出错的问题，谢谢
已经有9人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】求高手关于Fortran数据读写的问题。【已完结】
已经有12人回复
- 【求助】求指点，在Fortran里面，怎样实现“数组维度可调”的数组?
已经有5人回复
- 【求助】求助个Fortran计算，应该不难，半天搞不定，请高手帮忙
已经有11人回复
» 抢金币啦！回帖就可以得到:
-
医学超声影像负责人招聘-中国科学院赣江创新研究院
+1/975
-
医学超声影像负责人招聘-中国科学院赣江创新研究院
+1/974
-
西湖大学拓扑光学、非厄米光学、太赫兹方向博士后招聘
+2/236
-
西湖大学2026年秋季入学物理学、光学、电子信息方向博士生有名额速来！！！
+2/224
-
人间烟火，实则就是追求最简单的快乐
+1/74
-
一个陌生女人的来信
+1/73
-
海南大学海洋技术与装备学院-科研助理招聘（可读博）--膜分离水处理方向
+1/36
-
罗格斯大学纽瓦克校区(Rutgers-Newark) 招收 PHD，计算材料物理方向
+1/35
-
清华大学深圳国际研究生院招聘-博士后（长期有效）
+1/34
-
澳大利亚麦考瑞大学(Macquarie University)国际博士硕士全额奖学金-计算机-26年中开学
+1/15
-
上海交通大学-宁波东方理工大学联合培养博士生 – 力学
+1/14
-
上海交通大学-宁波东方理工大学联合培养博士生
+1/13
-
澳门大学生物医学影像实验室诚招博士生（2026秋季入学）
+1/7
-
全奖博士 英国利物浦大学 × 台湾清华大学 联合培养
+1/7
-
上海大学生物有机电子材料及器件团队博士研究生招聘
+1/6
-
德国图宾根大学诚招全奖岗位制博士（地下流固化学反应耦合数值模拟方向）
+1/5
-
26博士申请
+1/5
-
【科研助理招聘-北京理工大学-集成电路与电子学院-国家杰青团队】
+1/4
-
代教授（南昌航空大学）招收CO2光&光热催化还原方向的博士生
+1/4
-
澳科大招收2026年秋季药物递送/生物材料方向全奖博士研究生（3月5日18:00截止）
+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