Use the symmetry of the problem . With dexter's notation, notice that whenever s is a solution, so is 1/s. So, we can take his quartic:
z^4 + 2z^3 - 7z^2 + 2z + 1 = 0
And write down its four factors:
(z - s) (z - 1/s) (z - t) (z - 1/t) = 0
If we group them into their symmetric pairs and multiply:
(z^2 + az + 1) (z^2 + bz + 1)
So one quadratic corresponds to the two solutions s and 1/s and the other to t and 1/t. In particular, notice, now, that our unknowns only have two possible values, instead of four, which means we should be able to determine them from a quadratic equation:
Since the polynomials should be equal, we should expand and equate them:
z^4 + (a+b)z^3 + (2+ab)z^2 + (a+b)z + 1
=
z^4 + 2z^3 - 7z^2 + 2z + 1 = 0
Equating coefficients gives us a system of equations to solve.
Once we have a and b, we can then solve for s and t and get the four possible solutions. (Two should be complex, I imagine)
(Incidentally, I don't see how you got your equation, Recon)