我用INTEL FORTRAN 13.0 沒有問題
不過發現你的型態有錯誤
應該改為這樣子
a(i)=AA*3.0d0**i-0.4d0*(1.0d0/2.0d0)**i-0.1d0*i-0.05d0
的到的答案是
-0.350000000000000 1
-0.350000000000000 2
-0.400000000000000 3
-0.475000000000000 4
-0.562500000000000 5
-0.656250000000000 6
-0.753125000000000 7
-0.851562500000000 8
-0.950781250000000 9
-1.05039062500000 10
-1.15019531250000 11
-1.25009765625000 12
-1.35004882812500 13
-1.45002441406250 14
-1.55001220703125 15
-1.65000610351563 16
-1.75000305175781 17
-1.85000152587891 18
-1.95000076293945 19
-2.05000038146973 20
-2.15000019073486 21
-2.25000009536743 22
-2.35000004768372 23
-2.45000002384186 24
-2.55000001192093 25
-2.65000000596046 26
-2.75000000298023 27
-2.85000000149012 28
-2.95000000074506 29
-3.05000000037253 30