fix #523: finding z for |H(z)|=1 computed the wrong polynomials by bnavigator · Pull Request #525 · python-control/python-control