24小时热门版块排行榜

返回列表

jackwqh2010

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 598.5
帖子: 5
在线: 16.9小时
虫号: 1252958
注册: 2011-04-02
专业: 计算数学与科学工程计算

[求助] 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

求大神解答原因,谢谢！

回复此楼