| 查看: 575 | 回复: 0 | ||
jackwqh2010铜虫 (初入文坛)
|
[求助]
compaq Fortran 6.6 调用ImsL 子程序QDAGP
|
|
发现Fortran 调用QDAGP时，不能连续调用两次，否则就报错，代码如下： INTEGER NOUT, NPTS REAL A, ABS, ALOG, B, ERRABS, ERREST, ERROR, ERRREL, & EXACT, F, POINTS(2), RESULT, SQRT INTRINSIC ABS, ALOG, SQRT EXTERNAL F, QDAGP, UMACH C Get output unit number CALL UMACH (2, NOUT) C Set limits of integration A = 0.0 B = 3.0 C Set error tolerances ERRABS = 0.0 ERRREL = 0.01 C Set singularity parameters NPTS = 2 POINTS(1) = 1.0 POINTS(2) = SQRT(2.0) CALL QDAGP (F, A, B, NPTS, POINTS, ERRABS, ERRREL, RESULT, & ERREST) CALL QDAGP (G, A, B, NPTS, POINTS, ERRABS, ERRREL, RESULT, & ERREST) C Print results EXACT = 61.0*ALOG(2.0) + 77.0/4.0*ALOG(7.0) - 27.0 ERROR = ABS(RESULT-EXACT) c ERROR1 = ABS(RESULT1-EXACT) WRITE (NOUT,99999) RESULT,EXACT, ERREST, ERROR 99999 FORMAT (' Computed =', F8.3, 13X, &' Exact =', F8.3, /, /, & ' Error estimate =', 1PE10.3, 6X, 'Error =', 1PE10.3) C END C REAL FUNCTION F (X) REAL X REAL ABS, ALOG INTRINSIC ABS, ALOG F = X**3*ALOG(ABS((X*X-1.0)*(X*X-2.0))) RETURN END REAL FUNCTION G (X) REAL X REAL ABS, ALOG INTRINSIC ABS, ALOG G = X**3*ALOG(ABS((X*X-1.0)*(X*X-2.0))) RETURN END 求大神解答原因,谢谢！ |
回复此楼» 猜你喜欢
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有289人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
-
英国全奖博士招聘-深度学习与量子物理
已经有0人回复
-
间接带隙半导体有效质量求助
已经有0人回复
找到一些相关的精华帖子，希望有用哦~
fortran中最大整数问题
已经有7人回复
- visual studio 20a0 里面怎么运行fortran应该选择什么环境
已经有4人回复
- FORTRAN90标准函数及IMSL标准函数库
已经有166人回复
- FORTRAN中erf函数怎么编译？
已经有9人回复
- intel fortran XE版 调用imsl问题
已经有4人回复
- Intel Visual Fortran 里的 MKL 和 IMSL
已经有7人回复
- 【求助】Fortran如何调用模块中子程序
已经有6人回复
- 【求助】Intel Visual Fortran与Compaq Visual Fortran的区别？
已经有18人回复
- 【分享】带imsl数据库的fortran6.6下载
已经有9人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人