Astronomical Applications Department, U.S. Naval Observatory Sensitivity Integral Page 30
doubleprecision a4
doubleprecision t229
doubleprecision t231
doubleprecision t232
doubleprecision t235
doubleprecision t236
doubleprecision t237
doubleprecision t238
doubleprecision t239
doubleprecision t242
doubleprecision t243
t235 = delta**2
t243 = t235**2
t242 = delta*t235
t236 = nu1**2
t237 = t236**2
t238 = nu1*t237
t239 = nu1*t238
t232 = t242**2
t231 = delta*t243
t229 = t235+3.D0*delta+3.D0
G = 1.D0/(11880.D0*a2*(t232+7.D0*t231+21.D0*t243+35.D0*t242+35.D
#0*t235+21.D0*delta+7.D0)*t239+(10395.D0*a3*(t235+2.D0*delta+2.D0)*
#(t243+4.D0*t242+6.D0*t235+4.D0*delta+2.D0)*nu1*t239+13860.D0*a1*t2
#29*(t235+delta+1.D0)*t238)*(delta+2.D0)+(16632.D0*a0*(5.D0+10.D0*d
#elta+10.D0*t235+5.D0*t242+t243)+9240.D0*a4*t229*(t232+6.D0*t231+15
#.D0*t243+21.D0*t242+18.D0*t235+9.D0*delta+3.D0)*t237)*t237)
return
end
(
)
fortran
,
,
,
H optimized
=
mode
double
=
precision
double
doubleprecision function H(delta,nu1,a0,a1,a2,a3,a4,s,r,g)
doubleprecision delta
doubleprecision nu1
doubleprecision a0
doubleprecision a1
doubleprecision a2
doubleprecision a3
doubleprecision a4
doubleprecision s
doubleprecision r
doubleprecision g
doubleprecision t369
doubleprecision t373
doubleprecision t374
doubleprecision t376
doubleprecision t377
doubleprecision t378
doubleprecision t379
doubleprecision t381
doubleprecision t382
doubleprecision t383
doubleprecision t384
doubleprecision t386
doubleprecision t63
doubleprecision t64
Page 30
