This function stabilizes the polynomial transfer function.
b = polystab(a)
This is an Octave function. This function stabilizes the polynomial transfer function by replacing all roots outside the unit circle with their reflection inside the unit circle.