Написать программу по методу Брауна-Робинсона

private static void Braun(double[,] matrix, double eps)
        {
            int m = matrix.GetLength(0);
            int n = matrix.GetLength(1);
            List<double> B = new List<double>();
            List<double> cB = new List<double>();
            List<double> A = new List<double>();
            List<double> cA = new List<double>();
            for (int i = 0; i < n; i++) { B.Add(0.0); cB.Add(0.0); }
            for (int j = 0; j < m; j++) { A.Add(0.0); cA.Add(0.0); }
            Random rand = new Random();
            int f = rand.Next(m);            
            int min = 0;
            double Vn = 0.0, Vv = 0.0, Vs = 0.0, Vold = 0.0; 
            int count=0;
            int width=6;
            Console.Write("┌─────┬───┬"); for (int i = 0; i < B.Count; i++) Console.Write("──────┬"); Console.Write("───┬"); for (int i = 0; i < A.Count + 2; i++) Console.Write("──────┬"); Console.Write("──────┐");
            Console.WriteLine();
            Console.Write("│K".PadRight(width) + "│" + "i  │"); for (int i = 0; i < B.Count; i++) Console.Write("B:"+i.ToString().PadRight(width-2) + "│");
            Console.Write("j  │"); for (int i = 0; i < A.Count; i++) Console.Write("A:"+i.ToString().PadRight(width-2) + "│");
            Console.Write("Vn".PadRight(width)+"│"+"Vv".PadRight(width)+"│"+"Vs".PadRight(width)+"│");
            Console.WriteLine();
            while (true)
            {                
                Console.Write("├─────┼───┼"); for (int i = 0; i < B.Count; i++) Console.Write("──────┼"); Console.Write("───┼"); for (int i = 0; i < A.Count + 2; i++) Console.Write("──────┼"); Console.Write("──────┤");
                Console.WriteLine();
                count++;
                cA[f]++;
                for (int i = 0; i < B.Count; i++) B[i] += matrix[f, i];
                min = B.IndexOf(B.Min());
                cB[min]++;
                for (int j = 0; j < A.Count; j++) A[j] += matrix[j, min];
                Console.Write("│" + count.ToString().PadRight(width-1) + "│" + (f + 1).ToString().PadRight(3) + "│");
                for (int i = 0; i < B.Count; i++) Console.Write(Math.Round(B[i], 3).ToString().PadRight(width) + "│");
                Console.Write((min + 1).ToString().PadRight(3) + "│");
                for (int j = 0; j < A.Count; j++) Console.Write(Math.Round(A[j], 3).ToString().PadRight(width) + "│");
                Console.Write(Math.Round(Vn,3).ToString().PadRight(width)+"│"+Math.Round(Vv,3).ToString().PadRight(width)+"│"+Math.Round(Vs,3).ToString().PadRight(width)+"│");
                Console.WriteLine();
                Vn = B.Min() / count;
                Vv = A.Max() / count;
                Vs = (Vv + Vn)/2;                
                if ((Math.Abs(Vs - Vv)) < eps) break;
                //if ((count%20)==0) Console.ReadLine();
                Vold = Vs;
                f = A.IndexOf(A.Max());
            }            
            Console.Write("└─────┴───┴"); for (int i = 0; i < B.Count; i++) Console.Write("──────┴"); Console.Write("───┴"); for (int i = 0; i < A.Count + 2; i++) Console.Write("──────┴"); Console.Write("──────┘");
            Console.WriteLine();
            Console.WriteLine("V*= " + Math.Round(Vs, 3));
            double[] p = new double[A.Count];
            for (int i = 0; i < p.Length; i++) p[i] = cA[i] / count;
            Console.Write("p*= ( "); for (int i = 0; i < p.Length; i++) Console.Write(Math.Round(p[i],3)+"; "); Console.WriteLine(")");
            double[] q = new double[B.Count];
            for (int i = 0; i < q.Length; i++) q[i] = cB[i] / count;
            Console.Write("q*= ( "); for (int i = 0; i < q.Length; i++) Console.Write(Math.Round(q[i],3)+"; "); Console.WriteLine(")");            
        }