fix issue #8323

	harmonic(n).evalf() and catalan(n).evalf()
Yes, though maybe I'm barking up the wrong tree to fix the printing issue.

Should this be documented?  `e.as_coeff_mul()` does not behavie as I expected it would.