Solve (polygamma(0, 1 + 1/x) + x^2 polygamma(0, x) - (1 + x^2) polygamma(0, -1/x + x)) = 0, x