RosettaCodeData/Task/Pathological-floating-point-problems/PARI-GP/pathological-floating-point-problems-3.parigp

f(a,b)=333.75*b^6+a*a*(11*a*a*b*b-b^6-121*b^4-2)+5.5*b^8+a/(2*b)