| 查看: 586 | 回复: 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 求大神解答原因,谢谢！ |
回复此楼» 猜你喜欢
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有199人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有5人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
