solve(x*(csc(x) + 1) = pi) for x