Hi, I'm just learning to use Mathematica and I've defined the following function:
f[a_, b_] = 333.75*b^6 +
a^2*(11*a^2*b^2 - b^6 - 121*b^4 - 2) + 5.5*b^8 + a/(2*b)
f[77617, 33096] is supposed to be -0.8273960599468214 while Mathematica gives 1.1805916207174113e21. I've verified my result with Maple and GNU MP/MPFR with 1024 bit mantissa and both gave the number it is supposed to be.
Neither of single, double or extended precision floating point numbers was able to calculate this expression.
What could I be doing wrong?