Adding hybrid analytic / companion matrix polynomial root solver adapted from David Eberly's Wild Magic code