沒碰過這問題 不過稍微查了下跟INCLUDE有關係 可能你的環境沒裝好或是你的編譯器和你的檔案版本不一你的版本是不是f77 上網找 f90是D_QDAGdqdag ERRABS和ERRREL這邊是雙精度轉實數誤差?可是精度誤差 第一原來有編譯器會抓第二這不該讓程式跑不動啊由於沒有你的程式碼 我建議把include那行刪掉 但是你的程式運作的目的應該就是使用include的東西另一個可能的改法是 根據maxsub相關的敘述 你必須去更動maxsub的陣列大小 最好是改成動態配置 如果位於副程式（subroutine ）內，就是在宣告時以（:）進行宣告，問題是include的東西能不能改也是個問題
https://docs.roguewave.com/imsl/fortran/6.0/math/default.htm?turl=qdag.htm根據這個網站內提供的程式碼 error tolerancesERRABS = 0.0這個可能也要改