using System;

​

public class Program

{

    public static double[,] SwapRows(double[,] mat, int row1, int row2)

    {

        double temp;

        for (int j = 0; j < mat.GetLength(1); j++)

        {

            temp = mat[row1, j];

            mat[row1, j] = mat[row2, j];

            mat[row2, j] = temp;

        }

        return mat;

    }

​

    public static double[,] ScaleRow(double[,] mat, int row, double scalar)

    {

        for (int j = 0; j < mat.GetLength(1); j++)

        {

            mat[row, j] *= scalar;

        }

        return mat;

    }

​

    public static double[,] AddScaledRow(double[,] mat, int sourceRow, int targetRow, double scalar)

    {

        for (int j = 0; j < mat.GetLength(1); j++)

        {

            mat[targetRow, j] += scalar * mat[sourceRow, j];

        }

        return mat;

    }

​

    public static void PrintMatrix(double[,] mat)

    {

        for (int i = 0; i < mat.GetLength(0); i++)

        {

            for (int j = 0; j < mat.GetLength(1); j++)

            {

                Console.Write(mat[i, j] + "\t");

            }

            Console.WriteLine();

        }

    }

​

    public static void Main()

    {

        // Definir la matriz m

        double[,] m = new double[4, 5];

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

        {

            for (int j = 0; j < 5; j++)

            {

                m[i, j] = i * 1.0 * j;

            }

        }

​

        // Mostrar la matriz original

        Console.WriteLine("Matriz Original:");

        PrintMatrix(m);

        Console.WriteLine();

​

        // Aplicar operaciones de Gauss-Jordan

        m = SwapRows(m, 1, 2);

        m = ScaleRow(m, 2, 0.5);

        m = AddScaledRow(m, 1, 3, -2);

​

        // Mostrar la matriz después de las operaciones

        Console.WriteLine("Matriz Después de Operaciones de Gauss-Jordan:");

        PrintMatrix(m);

    }

}