using System;   //Programa 2, actividad 6

        /// Matriz aleatoria o manual           

public class Programa2JC

{

    public static void Main()

    {

    double[,] matriz;       // Es la matriz pero no se conoce el tamaño 

        int r, c;           // Para renglon y columna

        //solicitar tamaño de matriz

        Console.WriteLine("Cuantos renglones deseas? ");

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

        Console.WriteLine("Cuantas columnas deseas? ");

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

        matriz=new double [r,c];        // se declara matriz tipo double con el tamaño de usuario

        int n;

        Console.WriteLine("Deseas una matriz aleatoria o manual, 1=aleatorio o 2=manual");

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

        Random rnd= new Random();   //comando para generar numeros en orden aleatorio

​

        //Se piden valores de la matriz 

        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);

                if(n==0)

                {

                matriz[i,j]=double.Parse(Console.ReadLine());

                }

                else

                {

                matriz[i,j]=rnd.Next(-10,10);       //Comando genera numeros en orden aleatorio para matriz definiendo sus limites  

                }

            }

        }

        ////////Ciclo para visualizar la matriz//////

        Console.Write("\n**Matriz inicial");

            for(int i=0; i<r; i+=1)  // este es para que corra el renglón

             {

             Console.Write("\n");

                 for(int j=0;j<c;j+=1)     //// Este es para que corra las columnas

                 {

                    Console.Write("\t");

                   Console.Write(matriz[i,j]);

                 }

             }

        /////////////////////////////////////

        /////////////// se implementa el método de Gauss para resolver esta matriz

        Console.WriteLine("\n Aquí comienza metodo de Gauss Jordan para resolver matriz");

        double pivote, factor;

​

            for (int renglon = 0; renglon < r; renglon = renglon + 1)

            {

                pivote = matriz[renglon, renglon];                //este es el número de la diagonal y se llama pivote

​

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

                {

                    matriz[renglon, columna] = matriz[renglon, columna] / pivote;      //Divide todo el renglón entre el pivote

                }

                for (int renglon_eliminar = 0; renglon_eliminar < r; renglon_eliminar = renglon_eliminar + 1)    // elimina los elementos de la misma columna

                {

                    if (renglon_eliminar != renglon)

                    {

                        factor = matriz[renglon_eliminar, renglon];               // este es el número por el que se multiplicará el pivote para eliminar el renglón

​

                     for (int columna_eliminar = 0; columna_eliminar < c;

                            columna_eliminar = columna_eliminar + 1)

                        {

                            matriz[renglon_eliminar, columna_eliminar]

                                = matriz[renglon_eliminar, columna_eliminar] - factor              //Resta todo el renglón principal a el renglón a eliminar

                                * matriz[renglon, columna_eliminar];

                        }

                    }

​

        ///////////////////////////////////////////////////////////////////////////

        ////////Ciclo para visualizar la matriz terminada /////

        Console.Write("\n**Matriz resuelta");

            for(int i=0; i<r; i+=1)  // este es para que corra el renglón

             {

             Console.Write("\n");

                 for(int j=0;j<c;j+=1)     //// Este es para que corra las columnas

                 {

                    Console.Write("\t");

                   Console.Write(matriz[i,j]);

                     if(Double.IsNaN(matriz[i,j]))

                       Console.WriteLine("\n Este sistema NO tiene solución");

                 }

             }

        /////////////////////////////////////

        ///////////////////////////////

                }

            }

    }

    }