using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;

namespace MetodyNumeryczne
{
    public class MetodaGaussa
    {
        public static int Gauss(double [,] A,double [] B, double [] X, double eps )
        {
            //Wybór Elementu Głównego
            double T, ZS, MA;
            int k, i, j, N;
            N = A.GetLength(0) - 1;
            for (i = 1; i <= N; i++)
            {
                T = Math.Abs(A[i, i]);
                k = i;
                for (j = i + 1; j <= N; j++)
                {
                    MA = Math.Abs(A[j, i]);
                    if (MA > T)
                    {
                        T = MA;
                        k = j;
                    }
                }
                ZS = B[k];
                B[k] = B[i];
                B[i] = ZS;
                for (j = N; j >= 1; j--)
                {
                    ZS = A[k, j];
                    A[k, j] = A[i, j];
                    A[i, j] = ZS;
                }

            }

            //Przekształcenie Macierzy
            double m;
            for (i = 1; i <= N; i++)
            {
                for (j = i + 1; j <= N; j++)
                {
                    m = A[j, i] / A[i, i];
                    B[j] -= B[i] * m;
                    for (k = 1; k <= N; k++)
                        A[j, k] -= A[i, k] * m;
                }
            }
            //Metoda Eliminacji Wstecznej
            double Z;
            for (i = N; i >= 1; i--)
            {
                Z = B[i];
                for (j = i + 1; j <= N; j++)
                    Z -= A[i, j] * X[j];
                X[i] = Z / A[i, i];

            }

            for (j = 1; j <= B.GetLength(0) - 1; j++)
            {

            }

            return 0;

        }

    }
}