using System;

​

namespace Pavlenko_3PKS_318

{

    class Program

    {

        static void Swap<T>(ref T x, ref T y)

        {

            T t = y;

            y = x;

            x = t;

        }

        static bool checkUsl(ref double[] xk, ref double[] xkp, int n, double eps)

        {

            double norm = 0;

            for (int i = 0; i < n; i++)

                norm += (xk[i] - xkp[i]) * (xk[i] - xkp[i]);

            return (Math.Sqrt(norm) < eps);

        }

        static void matricaOutput(ref double[,] X, ref double[] B, int n)

        {

            Console.WriteLine("Вывод матрицы со стобцом свободных членов(округление значений до 4-x знаков)");

            for (int i = 0; i < n; i++)

            {

                for (int j = 0; j < n; j++)

                {

                    Console.Write("{0}\t", Math.Round(X[i, j], 4));

                }

                Console.Write("{0}\n", Math.Round(B[i], 4));

            }

            Console.WriteLine();

        }

​

        static void change(ref double[,] matrix, ref double[] b, int n, int num)

        {

            //Меняет строку num с другой строкой

            //для выполнения условия преобладания диагоналей

            //num - номер строки

​

            int numOfStr = -1;

            double[] stroka = new double[n];

​

            for (int i = num + 1; i < n; i++)

            {

                for (int j = 0; j < n; j++)

                {

                    stroka[j] = matrix[i, j];

                }

                if (numOfStr == -1)

                {

​

                    if (checkDiag(ref stroka, n, num))

                    {

                        numOfStr = i;

                    }

                }

            }

​

            //Теперь меняем строки местами

            for (int j = 0; j < n; j++)

            {

                Swap(ref matrix[numOfStr, j], ref matrix[num, j]);

            }

            Swap(ref b[num], ref b[numOfStr]);

        }

        static bool checkDiag(ref double[] stroka, int n, int i)

        {

            //проверяет подходит ли строка с диаг.элементом i

            double tmp = Math.Abs(stroka[i]);

            double summ = 0;

            for (int j = 0; j < n; j++)

            {

                if (i == j)

                    break;

                summ += Math.Abs(stroka[j]);

            }

            bool tmpBl;

            if (tmp > summ)

                tmpBl = true;

            else

                tmpBl = false;

            return tmpBl;

        }

        static void Main(string[] args)

        {

            double eps = 0.005;

            Console.Write("Введите размерность матрицы n: ");

            int n = Convert.ToInt32(Console.ReadLine()); ;

            double[,] matrix = new double[n, n];

            double[] b = new double[n];

            Console.WriteLine("Ввод матрицы: ");

            for (int i = 0; i < n; i++)

            {

                for (int j = 0; j < n; j++)

                {

                    Console.Write("\nВведите элемент с индексом i: {0} j: {0} =", i + 1, j + 1);

                    matrix[i, j] = Convert.ToDouble(Console.ReadLine());

                }

​

            }

            Console.WriteLine("\nВвод столбца свободных членов: ");

            for (int i = 0; i < n; i++)

            {

                Console.WriteLine("\n{0}-й элемент: ", i + 1);

                b[i] = Convert.ToDouble(Console.ReadLine());

            }

            Console.WriteLine("Введённая матрица:");

            matricaOutput(ref matrix, ref b, n);

            Console.WriteLine("\n");

​

            //Проверка условия

            for (int i = 0; i < n; i++)

            {

                double tmp = Math.Abs(matrix[i, i]);

                double summ = 0;

                for (int j = 0; j < n; j++)

                {

                    if (i != j)

                        summ += Math.Abs(matrix[i, j]);

​

                }

                if (tmp < summ)

                {

                    change(ref matrix, ref b, n, i);

                }

            }

            Console.WriteLine("Матрица после преобразований: ");

            matricaOutput(ref matrix, ref b, n);

            Console.WriteLine("\n");

​

            double[] Xt = new double[n]; //Текущий массив

            double[] Xpred = new double[n]; //Пред. масс

​

            for (int i = 0; i < n; i++)

            {

                double tmp = Math.Abs(matrix[i, i]);

                double summ = 0;

                for (int j = 0; j < n; j++)

                {

                    if (i != j)

                        summ += Math.Abs(matrix[i, j]);

​

                }

                if (tmp < summ)

                {

                    break;

                }

            }

            for (int i = 0; i < n; i++)

            {

                b[i] /= matrix[i, i];

                double diagElem = matrix[i, i];

                for (int j = 0; j < n; j++)

                {

                    matrix[i, j] /= diagElem;

                }

            }

            for (int i = 0; i < n; i++)

            {

                Xt[i] = b[i]; // В качестве начального приближения возьмём столбец свободных членов

                Xpred[i] = b[i]; //т.к первая итерация

            }

            matricaOutput(ref matrix, ref b, n);

            int m = 0;

            do

            {

                for (int i = 0; i < n; i++)

                    Xpred[i] = Xt[i];

                for (int i = 0; i < n; i++)

                {

                    double var = 0;

                    for (int j = 0; j < i; j++)

                        var += (matrix[i, j] * Xt[j]);

                    for (int j = i + 1; j < n; j++)

                        var += (matrix[i, j] * Xpred[j]);

                    Xt[i] = (b[i] - var) / matrix[i, i];

                }

                m++;

            } while (!checkUsl(ref Xt, ref Xpred, n, eps));

​

            for (int i = 0; i < n; i++)

                Console.WriteLine("X{0} = {1}", i + 1, Math.Round(Xt[i], 4)); //вывод результата

            Console.ReadKey();

        }

​

    }

​

​

}

​