using System;

public class Program

{

    public static void Main()

    {

        Console.WriteLine("Bienvenido amigo");

        //Declarar Variables

            double x1, x2, x3, f1, f2,f3, 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, que en este caso es una matriz de 3X4

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

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

          int iteracion = 0;

        //Se da un valor inicial se deduce.

            x1 = 100;

            x2 = -100;

            x3 = 100;

            do

            {

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

                f1 = x2 * Math.Cos(x1 + x2);

                f2 = x1 * Math.Exp(x1 + x2);

                f3 = (x1 * x3 )+ (x2 * x3);

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

                matriz[0, 0] = (2 * Math.Cos(3*x1));

                matriz[0, 1] = -( 3* x1) * Math.Sin(3*x1);

                matriz[0, 2] = -f1;

                matriz[1, 0] = Math.Exp(3 * x2);

                matriz[1, 1] = 3 * Math.Exp(3*x2) * x2;

                matriz[1, 2] = -f2;

                matriz[2, 0] = 4 * x3;

                matriz[2, 1] = 0;

                matriz[2, 2] = -f3;

​

                //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] = 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. 

                x1 = x1 + matriz[0, 2];

                x2 = x2 + matriz[1, 2];

                x3 = x3 + matriz[2, 2];

                iteracion++;

            } while (Math.Abs(f1) > Epsilon | Math.Abs(f2) > Epsilon | Math.Abs(f3) > 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 x1 es: " + x1);

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

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

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

    }

}