| 查看: 1104 | 回复: 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 |
» 猜你喜欢
申博发邮件
已经有8人回复
各位大神,目前国内有哪些比较好用的逆合成软件?
已经有9人回复
国社科系统bug了,是不是要放榜了?
已经有7人回复
上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告
已经有8人回复
上海工程技术大学激光智能制造课题组招收博士研究生
已经有8人回复
» 本主题相关价值贴推荐,对您同样有帮助:
请高手指点一下怎么确定一个符号矩阵中的等值元素
已经有6人回复
陌生的fortran语言,哪位可以指点一二?悬赏100金币
已经有5人回复
求助计算机高人,编写个Fortran程序
已经有13人回复
请教一个fortran小程序编译出错的问题,谢谢
已经有9人回复
写了一个fortran90的小程序,编译通不过,请大侠帮忙
已经有59人回复
【求助】求高手关于Fortran数据读写的问题。【已完结】
已经有12人回复
【求助】求指点,在Fortran里面,怎样实现“数组维度可调”的数组?
已经有5人回复
【求助】求助个Fortran计算,应该不难,半天搞不定,请高手帮忙
已经有11人回复
» 抢金币啦!回帖就可以得到:
测试█TEM/ EPR/ XPS/PY-GCMS/TG-IR/XRF/BET/MIP/核磁/EA/ICP,VX: 761711562。
+1/98
【专业EPR】:自由基(超氧、羟基、单线态氧),空位缺陷;定量计算,VX: 761711562。
+1/97
转让有机合成实验室
+1/48
中国地质大学(北京)地质资源与地质工程专业 招收2027级硕士研究生
+1/26
【招聘信息】设备设计工程师,中科院孵化企业,技术大牛带队,新材料装备方向,2-3名
+1/16
【推免招生】招收材料、化工、环境、冶金类推免生
+1/12
(捡漏)中山大学2027年接收推荐优秀免试研究生报名通知(第二批)
+1/12
凝聚青年力量,共筑学术新程,Biofunctional Materials开启青年编委全球招募
+1/12
激荡学术火花,赋能青年成长,Biofunctional Materials青年编委招募
+1/11
信息工程大学推免地方研究生招生[9月18日截止,预报从速]
+1/9
课题组/公司长期承接多肽定制合成服务,欢迎有需要的虫友联系。
+1/7
寻找博士后的机会
+1/6
集结青年新锐,共创学术沃土,Biofunctional Materials 面向全球招募青年编委
+1/5
韩国全北国立大学禽病研究所博士研究生招生简章
+2/4
美国Ohio University医学院骨性关节炎实验室招收生物医学全奖博士生-2027年秋入学
+1/4
【招募学术报告嘉宾】青年学者论坛,线上讲座嘉宾招募中(颁发演讲证书)
+1/4
【生物材料与生物打印专题征稿】Scopus期刊,免APC
+1/3
柔性器件多物理场仿真:COMSOL热-力-电耦合与代理模型
+1/3
碎碎念
+1/3
口腔特刊客座招募客编Biofunctional Materials牙科生物功能材料专题,Scopus期刊
+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










回复此楼