I solved it by giving 2 complex equations for 2 unknowns zeta and omegaN, as there are 2 complex conjugate poles. And now Mathematica immediately finds the solution. This method also works for 2 real poles: Full Mathematica 13.2 notebook: