fcn rump(a,b){ b=BN(b);
   b2,b4,b6,b8:=b.pow(2),b.pow(4),b.pow(6),b.pow(8);
   a2:=BN(a).pow(2);
   r:=( b6*33375 + a2*(a2*b2*11 - b6 - b4*121 - 2)*100 + b8*550 )*ten2n;
   r+=BN(a)*ten2n*100/(2*b);
   tostr(r,66,32)
}
println("\n3rd: Rump's example: f(77617.0, 33096.0) = ",rump(77617,33096));