using System;

public class Program

{

    public static void Main()

    {       

        //double[,] matriz = { { 3, -1, 2, 4, 30 }, //ojo siempre hay que actualizar el tamaño de la matriz

        //                     { 1, -2, 3, 1, 20 },

        //                     { 3, 1, 3, 2, 10 },

        //                   { 1, 3, -3, 2, 30 }};

        double pivote, factor;

        int nv = 4; // Este es el número de variables

        Console.Write("Cuántas variables tiene el sistema?");

        nv = int.Parse(Console.ReadLine());

        Console.WriteLine();

        double [,] matriz = new double [nv, nv+1];

        for (int reng = 0; reng < nv; reng++)

        {

            Console.WriteLine("Dame los valores de la ecuación #" + (reng+1));

            for(int colu = 0; colu < nv+1; colu++)

            {

                if(colu != nv)

                {

                    Console.Write("Escribe el coeficiente de la variable X"+(colu+1)+": ");

                }

                else

                {

                    Console.Write("Ahora escribe el valor de la constante para esta ecuación: ");

                }

                matriz[reng,colu] = double.Parse(Console.ReadLine());               

            }

            Console.WriteLine();

        }

            //aquí empieza el método de Gauss

            for (int reng = 0; reng < nv; reng = reng + 1)

            {

                //Selecciona el elemento diagonal como pivote

                pivote = matriz[reng, reng];

​

                //Divide todo el renglón entre el pivote

                for(int colu= 0; colu < nv+1; colu++)

                {

                matriz[reng,colu]= matriz[reng,colu]/pivote;

                }

                //

                for(int reng_elimi = 0; reng_elimi < nv; reng_elimi++)

                {

                    if(reng_elimi != reng)

                    {

                        factor = matriz[reng_elimi,reng];

                        for(int colu_elimi = 0; colu_elimi < nv+1; colu_elimi++)

                        {

                            matriz[reng_elimi, colu_elimi] = matriz[reng_elimi,colu_elimi]-

                                factor*matriz[reng,colu_elimi];

                        }

                    }

                }

            }

        //aquí se muestra el resultado en la pantalla

            for (int reng = 0; reng < nv; reng++)

            {

                for(int colu= 0; colu < nv+1; colu++)

                {

                    Console.Write(matriz[reng,colu] + "       ");

                }

                Console.WriteLine();

            }

​

    }

}