版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

jackwqh2010

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 598.5
帖子: 5
在線: 16.9小時
蟲號: 1252958
注冊: 2011-04-02
專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算

[求助] compaq Fortran 6.6 調(diào)用ImsL 子程序QDAGP

發(fā)現(xiàn)Fortran 調(diào)用QDAGP時，不能連續(xù)調(diào)用兩次，否則就報錯，代碼如下：
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

求大神解答原因,謝謝！

回復(fù)此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有107人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有76人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)

1樓 2015-06-07 22:25:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

fortran中最大整數(shù)問題已經(jīng)有7人回復(fù)
visual studio 20a0 里面怎么運行fortran應(yīng)該選擇什么環(huán)境已經(jīng)有4人回復(fù)
FORTRAN90標(biāo)準(zhǔn)函數(shù)及IMSL標(biāo)準(zhǔn)函數(shù)庫已經(jīng)有166人回復(fù)
FORTRAN中erf函數(shù)怎么編譯？已經(jīng)有9人回復(fù)
intel fortran XE版調(diào)用imsl問題已經(jīng)有4人回復(fù)
Intel Visual Fortran 里的 MKL 和 IMSL 已經(jīng)有7人回復(fù)
【求助】Fortran如何調(diào)用模塊中子程序已經(jīng)有6人回復(fù)
【求助】Intel Visual Fortran與Compaq Visual Fortran的區(qū)別？已經(jīng)有18人回復(fù)
【分享】帶imsl數(shù)據(jù)庫的fortran6.6下載已經(jīng)有9人回復(fù)

點擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 jackwqh2010 的主題更新