Signal Toolkit - polystab

 b = polystab(a)

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