An algorithm for the numerical solution of an important algebraic matrix equation in control system design.