using System;

public class Program

{

    public static void Main()

    {

    int filas,columnas;

    string sfila,scolumna;

    Console.Write("¿Cuántas filas quieres en la matriz?");

    sfila=Console.ReadLine();

    filas=Convert.ToInt32(sfila);

    Console.Write("¿Cuántas columnas quieres en la matriz?");

    scolumna=Console.ReadLine();

    columnas=Convert.ToInt32(scolumna);

    double[,]matriz=new Double[filas,columnas];

    int[,]matrizR=new int[filas,columnas];

    for(int f=0;f<filas;f++)

  {

      for(int c=0;c<columnas;c++)

      {

        Console.WriteLine("\nEscribe un numero");

        matriz[f, c]=Convert.ToInt32(Console.ReadLine());

      }

    }

    Console.Write("\n La Matriz es : \n");

    for(int f=0;f<filas;f++)

  {

      for(int c=0;c<columnas;c++)

      {

    Console.Write("" +matriz[f, c]);

      }

        Console.WriteLine();

    }

        Console.ReadLine();

    double pivote, factor;

​

    for (int f = 0; f < 3; f = f + 1)

            {

                pivote = matriz[f, f];

                if (pivote==0)

                {

                    double suma=0;

                    for  (int x=f; x<3; x=x+1)

                          {

                              suma=suma+matriz[f,x];

                    }

                          if (suma==0)

                          {

                            Console.Write("Matriz sin solución");

                            Console.ReadLine();

                          }

                          else

                          {

                            Console.Write("La Matriz tiene una cantidad infinita de soluciones");

                            Console.ReadLine();

                          }

                          }

                         for (int c=0; c<4;c=c+1)

                          {

                              matriz[f,c]=matriz[f,c]/pivote;

                          }

                for (int f_elimi = 0; f_elimi < 3; f_elimi = f_elimi + 1)

                {

                   if (f_elimi !=f)

                   {

                       factor=matriz[f_elimi, f];

                       for (int c_elimi = 0; c_elimi < 3; c_elimi = c_elimi + 1)

                           {

                               matriz[f_elimi,c_elimi]=(matriz[f_elimi,c_elimi])-(factor*matriz[f,c_elimi]);

                           }

                   }

                }

                     Console.Write("Resultado Gauss");   

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

                          {

                              for(int c=0;c<4;c++)

                              {  

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

                                                  }

                                                  Console.WriteLine();

                                                  }

                                                  Console.ReadLine();

                                                  }

                                                Console.ReadLine();

    }

}

​

​

​