| 查看: 581 | 回复: 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 求大神解答原因,谢谢！ |
回复此楼» 猜你喜欢
光学工程学硕调剂信息
已经有32人回复
-
欢迎加入课题组
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有88人回复
-
散金币，求好运，祝面上顺利！
已经有79人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
