With solve(2*cot(x)==0, x, to_poly_solve = 'force') you will get all solutions, where $\pi z_{xx}$ should be interpreted as a multiple integer of $\pi$.