public static void Main()
Console.WriteLine("Cuantos renglones deseas? ");
r=int.Parse(Console.ReadLine());
Console.WriteLine("Cuantas columnas deseas? ");
c=int.Parse(Console.ReadLine());
Console.WriteLine("Deseas una matriz aleatoria o manual, 1=aleatorio o 2=manual");
n=int.Parse(Console.ReadLine());
Random rnd= new Random();
Console.WriteLine("Introduce valores de matriz");
for(int i=0 ; i < r ; i++)
for(int j=0; j < c ; j++)
Console.WriteLine("introduce valor de la entrada ({0},{1}", i,j);
matriz[i,j]=double.Parse(Console.ReadLine());
matriz[i,j]=rnd.Next(-10,10);
Console.Write("\n**Matriz inicial");
Console.Write(matriz[i,j]);
Console.WriteLine("\n Aquí comienza metodo de Gauss Jordan para resolver matriz");
for (int renglon = 0; renglon < r; renglon = renglon + 1)
pivote = matriz[renglon, renglon];
for (int columna = 0; columna < c; columna = columna + 1)
matriz[renglon, columna] = matriz[renglon, columna] / pivote;
for (int renglon_eliminar = 0; renglon_eliminar < r; renglon_eliminar = renglon_eliminar + 1)
if (renglon_eliminar != renglon)
factor = matriz[renglon_eliminar, renglon];
for (int columna_eliminar = 0; columna_eliminar < c;
columna_eliminar = columna_eliminar + 1)
matriz[renglon_eliminar, columna_eliminar]
= matriz[renglon_eliminar, columna_eliminar] - factor
* matriz[renglon, columna_eliminar];
Console.Write("\n**Matriz resuelta");
Console.Write(matriz[i,j]);
if(Double.IsNaN(matriz[i,j]))
Console.WriteLine("\n Este sistema NO tiene solución");