b = polystab(a)

Stabalize the polynomial transfer function by replacing all roots
outside the unit circle with their reflection inside the unit circle.