Решение систем линейных уравнений методами Гаусса и Крамера

using System;
using System.Collections.Generic;
using System.Text;
namespace ConsoleApplication_Kramer
{
class Program
{
static void Main(string[] args)
{
int n; /* количество уравнений */
double [,] A = new double [3,3]; /* матрица системы */
double [] b = new double [3]; /* вектор правых частей */
double [] x = new double [3]; /* вектор решения */
char qq;
Console.Write("Введите количество уравнений(<=3) n -> ");
n = Convert.ToInt32(Console.ReadLine());
if (n > 3 || n <= 1)
{
Console.WriteLine("Ошибка в размерности системы (n=2,3)");
Convert.ToInt32(Console.ReadLine());
return;
}
for(int i=0; i<n; i++)
for(int j=0; j<n; j++)
{
Console.Write("A{0}{1} -> ",i,j);
A[i,j] = Convert.ToDouble(Console.ReadLine());
}
for(int i=0; i<n; i++)
{
Console.Write("b{0} -> ", i);
b[i] = Convert.ToDouble(Console.ReadLine());
}
if(SLAU_kramer(n,A, b, x)==1)
{
Console.WriteLine("Система не имеет решение");
Convert.ToInt32(Console.ReadLine());
return;
}
else
for(int i=0; i<n; i++)
Console.WriteLine("x"+i+" = "+x[i]);
Console.ReadLine();
}
private
static double det(int n, double [,]B)
{
if (n == 2)
return B[0,0] * B[1,1] - B[0,1] * B[1,0];
return B[0,0] * (B[1,1] * B[2,2] - B[1,2] * B[2,1]) - B[0,1] * (B[1,0] * B[2,2] - B[1,2] * B[2,0])+
B[0,2]*(B[1,0] * B[2,1] - B[1,1] * B[2,0]);
}
static void equal(int n, double [,]A, double [,]B)
{
for(int i=0; i<n; i++)
for(int j=0; j<n; j++)
A[i,j]=B[i,j];
}
static void change(int n, int N, double[,] A, double[] b)
{
for(int i=0; i<n; i++)
A[i,N]=b[i];
}
public
static int SLAU_kramer(int n, double[,] A, double[] b, double[] x)
{
double [,]An = new double [3,3];
double det1 = det(n, A);
if (det1 == 0) return 1;
for (int i = 0; i < n; i++)
{
equal(n, An, A);
change(n, i, An, b);
x[i] = det(n, An) / det1;
}
return 0;
}
}
}

using System;
 
class Gauss
{
    public static int Main(string[] args)
    {
        try
        {
            Matrix a = (args.Length == 0) ? SLAU.In() : SLAU.In(args[0]);
 
            for (int k = 1; k <= a.AmountOfRows - 1; k++)
            {
                if (a[k, k] == 0)
                    for (int i = k + 1; i <= a.AmountOfRows; i++)
                        if (a[i, k] != 0)
                        {
                            a.SwapRows(k, i);
                            break;
                        }
 
                double delElem = a[k, k];
                if (delElem == 0)
                {
                    Console.WriteLine("Система не совместна");
                    return 0;
                }
                a.MulRow(k, 1 / delElem);
                for (int i = k + 1; i <= a.AmountOfRows; i++)
                    a.SumRows(k, i, -a[i, k]);
            }
 
            if (a[a.AmountOfRows, a.AmountOfColumns - 1] == 0)
            {
                Console.WriteLine("Система не совместна");
                return 0;
            }
            a.MulRow(a.AmountOfRows, 1 / a[a.AmountOfRows, a.AmountOfColumns - 1]);
 
            Matrix res = new Matrix(a.AmountOfRows, 1);
 
            for (int i = a.AmountOfRows; i >= 1; i--)
            {
                double sum = 0;
                for (int j = a.AmountOfRows; j >= i + 1; j--)
                    sum += res[j, 1] * a[i, j];
                res[i, 1] = a[i, a.AmountOfColumns] - sum;
            }
            res.Out();
        }
        catch (Exception exc)
        {
            Console.WriteLine("Произошла критическая ошибка. 
 Дальнейшее выполение программы нецелесобразно!!!");
            Console.WriteLine(exc.StackTrace);
            Console.WriteLine(exc.TargetSite);
            return 1;
        }
        finally
        {
            Console.WriteLine();
            Console.Write("Для продолжения нажмите Enter ... ");
            Console.ReadLine();
        }
        return 0;
    }
}