(defn rump [a b]
  (+ (* (rationalize 333.75) (expt b 6))
      (* (expt a 2)
          (- (* 11 (expt a 2) (expt b 2)) (expt b 6) (* 121 (expt b 4)) 2))
      (* (rationalize 5.5) (expt b 8))
      (/ a (* 2 b))))
; Using BigInt numeric literal style to avoid integer overflow
(double (rump 77617 33096N))