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 = x * x + y * y - 4;

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

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

                matriz[0, 0] = 2 * x;

                matriz[0, 1] = 2 * y;

                matriz[0, 2] = -f1;

                matriz[1, 0] = Math.Exp(x);

                matriz[1, 1] = 1;

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

        }

}