using System;

public class Program

{

    public static void Main()

    {

        //Comenzamos estableciendo las variables 

        double X1,X2,X3, f1, f2,f3, E=0.001;//E es Epsilon es el error aceptado

        double pivote, factor;//Variables de gauss Factor y Pivote 

        //Para este casose requiere una matriz de 3*4

        double [,] matriz = new double [3,4];//La declaramos

        int iteracion = 0;//Con esto establecemos las iteraciones

        X1= -100;//introducimos loS valores iniciales 

        X2= -100;

        X3= -100;

        do//aplicamos un ciclo do while

        {//calculamos  las primeras operaciones y las resolvemos

            f1 = 60.9324975239875*Math.Pow(Math.E,(0.0820527099405327*X1)); 

            f2 = 339.688692801301*Math.Pow(Math.E,(0.133346205242556*X2));

            f3 = 892.945917213366*Math.Pow(Math.E,(0.222060685228468*X3));

//Comenzamos con la matriz jacobiana

                matriz[0, 0] = 4.99951*Math.Pow(Math.E,(0.0820527099405327*X1));

                matriz[0, 1] = 0;

                matriz[0, 2] = -f1;

                matriz[1, 0] =  45.29408*Math.Pow(Math.E,(0.133346205242556*X2));

                matriz[1, 1] = 0;

                matriz[1, 2] = -f2;

                matriz[2, 0] =  198.28756*Math.Pow(Math.E,(0.222060685228468*X3));

                matriz[2, 1] = 0;

                matriz[2, 2] = -f3;

//pasamos a la matriz gauss

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

                            }

                        }

                    }

                }

                //comienza el conteo de iteraciones

                X1 = X1 + matriz[0, 2];

                X2 = X2 + matriz[1, 2];

                X3 = X3 + matriz[2, 2];         

                iteracion++;

        }//se genra la iteración 

        while(Math.Abs(f1) > E | Math.Abs(f2) > E| Math.Abs(f3) > E); //mostramos los resultados obtenidos

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

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

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

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

    }

}