Метод простых итераций

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
 
namespace TheMethodOfSimpleIterations
{
    class Program
    {
        static void Main(string[] args)
        {
            try
            {
                Console.Write("Точность> ");
                var p = int.Parse(Console.ReadLine());
 
                Console.Write("\nКоличество переменных> ");
                var n = int.Parse(Console.ReadLine());
 
                Console.WriteLine("\nВвод матрицы коэффициентов (разделяя пробелами):\n");
                var a = InputA(n);
 
                Console.WriteLine("\nВвод вектора свободных членов (разделяя пробелами):\n");
                var b = InputB(n);
 
                if (b.Length != n)
                    throw new ArgumentException("b");
 
                GetIterationForm(a, b);
 
                if (!CheckForConvergence(a, b))
                {
                    Console.Write("\nДанная система не удовлетворяет условию сходимости.");
                    Console.ReadKey();
                    return;
                }
 
                var result = GetSolve(a, b, p);
                var s = string.Join(" ", result);
                Console.Write("\nРешение системы равно вектору: {0}", s);
            }
            catch (Exception e)
            {
                Console.Write("\nНекоректный ввод.");
            }
            Console.ReadKey();
        }
 
        static double[,] InputA(int n)
        {
            var a = new double[n, n];
 
            for (int i = 0; i < n; i++)
            {
                Console.Write("[{0}]> ", i);
                var s = Console.ReadLine();
                var numbers = s.Split(' ');
 
                for (int j = 0; j < n; j++)
                    a[i, j] = double.Parse(numbers[j]);
            }
 
            return a;
        }
 
        static double[] InputB(int n)
        {
            Console.Write("> ");
            var s = Console.ReadLine();
            return s.Split(' ').Select(value => double.Parse(value)).ToArray();
        }
 
        static void GetIterationForm(double[,] a, double[] b)
        {
            var c = a.Clone() as double[,];
            var d = b.Clone() as double[];
 
            for (int i = 0; i < a.GetLength(0); i++)
            {
                for (int j = 0; j < a.GetLength(1); j++)
                    a[i, j] = i == j ? 0 : -c[i, j] / c[i, i];
 
                b[i] = d[i] / c[i, i];
            }
        }
 
        static bool CheckForConvergence(double[,] a, double[] b)
        {
            var sum = 0.0;
 
            foreach (var n in a)
                sum += n;
 
            return Math.Sqrt(Math.Abs(sum)) > 1;
        }
 
        static double[] GetSolve(double[,] a, double[] b, int p)
        {
            var min = Math.Pow(10, p);
            var delta = 0.0;
            var x = new double[a.GetLength(0)];
            var x0 = new double[a.GetLength(0)];
            var e = new double[a.GetLength(0)];
 
            do
            {
                for (int i = 0; i < a.GetLength(0); i++)
                {
                    x[i] = 0;
 
                    for (int j = 0; j < a.GetLength(1); j++)
                        x[i] += a[i, j] * x0[j];
 
                    x[i] += b[i];
                    e[i] = Math.Abs(x[i] - x0[i]);
                }
 
                delta = e[0];
 
                for (int i = 0; i < e.GetLength(0); i++)
                    if (delta < e[i])
                        delta = e[i];
 
                x0 = x.Clone() as double[];
            }
            while (delta <= min);
 
            return x;
        }
    }
}