public static double SpectralRadius(MMatrix A)
        {
            Vector Eigenvalues = new Vector();

            if (!A.IsSquared())
                throw new MMatrixException("Matrix must be squared.");
            Eigenvalues = A.QRIterationBasic(GlobalMath.MAXITER).DiagVector();
            Eigenvalues.Sort(true);

            return Math.Abs(Eigenvalues[0]);
        }