using System;

public class Program

{

    public static void Main()

    {

        //Declarar Variables

double x, y, f1, f2, Epsilon = 0.001;

//Factor y Pivote variables heredadas del método de Gauss Jordan

double pivote, factor;

//Se declara la matriz aumentada para utilizar dependiendo del sistema de ecuaciones

// Si son dos ecuaciones y 2 variables  es una matriz 2 x 3 si es un sistema de 3 ecuaciones con 3 variables es una matriz de 3 x 4

double[,] matriz = new double[2, 3];

//se declara una variable iteracion para ver cuantos pasos realizo

int iteracion = 0;

​

//Se da un valor inicial se deduce.

x = -1;

y = -1.7;

​

​

do

{

    // Se calculan  los valores delas variables en las  ecuaciones o fucniones

    f1 = 2 * x - y - Math.Exp(-x);

    f2 = -x + 2 * y - Math.Exp(-y);

​

    // Se declaran las operaciones en nuestra ecuacion Jacobi y en la columna aumentada multiplicarla por -1

    matriz[0, 0] = 2 * -Math.Exp(-x) ;

    matriz[0, 1] = 1;

    matriz[0, 2] = -f1;

    matriz[1, 0] = -1;

    matriz[1, 1] = 2 - Math.Exp(-y);

    matriz[1, 2] = -f2;

​

    //Eliminación Gaussiana Gauss Jordan//------------------------------

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

    {

        pivote = matriz[reng, reng];

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

        {

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

        }

        for (int reng_elimi = 0; reng_elimi < 2; reng_elimi = reng_elimi + 1)

        {

            if (reng_elimi != reng)

            {

                factor = matriz[reng_elimi, reng];

                for (int colu_elimi = 0; colu_elimi < 3;

                    colu_elimi = colu_elimi + 1)

                {

                    matriz[reng_elimi, colu_elimi] = matriz[reng_elimi,

                       colu_elimi] - factor * matriz[reng, colu_elimi];

                }

            }

        }

    }

    //------------------------------

​

    //se hace el incremento de la raiz en una nueva iteracion. 

    x = x + matriz[0, 2];

    y = y + matriz[1, 2];

    iteracion++;

} while (Math.Abs(f1) > Epsilon | Math.Abs(f2) > Epsilon); //Se toma el Absoluto de la funcion debe ser mayo que el error de tolerancia Epsilon

                                                           // Se imprimen los resultados en de los valores. 

Console.WriteLine("El valor final de x es: " + x);

Console.WriteLine("El valor final de y es: " + y);

Console.WriteLine("Total de Iteraciones es: " + iteracion);

    }

}

​