| 查看: 602 | 回复: 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 求大神解答原因,谢谢！ |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有193人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
找到一些相关的精华帖子，希望有用哦~
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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
