Rumpの例題

Rumpの例題
a = 77617,  b = 33096;
に対して
f = (333.75 - a^2) * b^6 + a^2 * (11 * a^2 * b^2 -121 * b^2 -121 * b^4 - 2) + 5.5 * b^8 + a / (2 * b)

を計算する.
 
fの正しい結果は
 f = -0.827396...
 
一方で,これらをOctaveFortranで計算してみると...
 
Octaveのプログラム(prog01.m)
format long;
a = 77617; b = 33096;
f = (333.75 - a^2) * b^6 + a^2 * (11 * a^2 * b^2 -121 * b^2 -121 * b^4 - 2) + 5.5 * b^8 + a / (2 * b)
 
実行結果
> prog01
f =  1.172603940053179
 
答えが違う....符号すら合っていません.
 
Fortranのプログラム(prog02.f90)
program prog02
implicit none
real(kind(1E0)):: ae,be,fe
real(kind(1D0)):: ad,bd,fd
real(kind(1Q0)):: aq,bq,fq
ae = 77617E0; be = 33096E0
ad = 77617D0; bd = 33096D0
aq = 77617Q0; bq = 33096Q0
fe = (333.75E0 - ae**2) * be**+ ae*** (11 * ae*** be**2 - 121 * be**4 - 2) + &
        5.5E0 * be**+ ae / (* be)
fd = (333.75D- ad**2) * bd**+ ad*** (11 * ad*** bd**2 -121 * bd**4 - 2) + &
        5.5D* bd**+ ad / (* bd)
fq = (333.75Q- aq**2) * bq**+ aq*** (11 * aq*** bq**2 -121 * bq**4 - 2) + &
        5.5Q* bq**+ aq / (* bq)
write(*,*) fe
write(*,*) fd
write(*,*) fq
end program prog02
 
実行結果
$ gfortran prog02.f90
$ ./a.out
   1.17260396    
   1.1726039400531787     
   1.17260394005317863185883490452018380      
 
これも答えが違い,符号すら合っていません....