polygamma: fix expand(func=True) (#643)