using System;

public class Program

{

    public static void Main()

    {

            //Variables de Entrada 

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

        // pivote y factores variables heredadas 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];

​

        // Paso 1 se darán los valores iniciales a ojo de buen cubero

            x = 1;

            y = -1.7;

        // Se calculan las funciones con  los valores inciales de las variables

            f1 = x * x + y * y - 4;

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

​

          //Mientras las ecuaciones iniciaes sigan siendo mayores a un numero muy aproximado a cero

            while (Math.Abs(f1) > epsilon | Math.Abs(f2) > epsilon) 

            {

                // Se vuelve a calcular las variables ya que estas cambiaran 

                f1 = x * x + y * y - 4;

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

                //Se colocan la matriz jacobi con la columna aumentada multiplicada 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 //------------------------------

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

                            }

                        }

                    }

                }

                //------------------------------ Termina Gauss Jordan ------------------------------

                // Se aplica el incremento de los nuevos valores obtenidos  y se evaluan en la nueva matriz

                x = x + matriz[0, 2];

                y = y + matriz[1, 2];

​

            }

            //Se Imprimen las Raices

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

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

​

            Console.ReadLine();

        }

    }

​