I want to understand the polygamma reflection. I can get the polynomials in cos^2. Go into test_polygamma and push it through.