using System;

public class Program

{

    public static void Main()

    {

        //Ejercicio 2 elaborado por Omar Juárez y Franco Nuñez

            //Aqui declaramos variables

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

        double e= 2.71828;

            //Establecemos el factor y Pivote variables heredadas del método de Gauss Jordan

            double pivote, factor;

           //Declaramos la matriz aumentada para utilizar dependiendo del sistema de ecuaciones

           // En esta ocasion es un sistema de 3 ecuaciones con 3 variables por lo que es una matriz de 3 x 4

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

          //Declaramos una variable iteracion el numero de pasos

          int iteracion = 0;

            //Aqui deducimos los valores sin importar solo al tanteo.

            x = 35.9;

            y = 36.9;

            z = 6.9;

            do

            {

                 // Calculamos los valores de las variables en las  ecuaciones

                f1 = e*Math.Exp(x) + 2*Math.Exp(-y) + 2*Math.Cos(z) -6;

                f2 = x*Math.Exp(2)+10*Math.Cos(y)+z+3;

                // Declaramos las operaciones en nuestra ecuacion Jacobi y en la columna aumentada que será multiplicada por -1

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

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

                matriz[0, 2] = -2*Math.Sin(z);

                matriz[0, 3] = -f1;

                matriz[1, 0] = 2*x;

                matriz[1, 1] = -10*Math.Sin(y);

                matriz[1, 2] = 1;

                matriz[1, 3] = -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];

                            }

                        }

                    }

                }

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

                //Realizamos el incremento de la raiz en una nueva iteracion. 

                x = x + matriz[0, 3];

                y = y + matriz[1, 3];

​

                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

            // Arrojamos los resultados de los valores de x,y,z 

            Console.WriteLine("El valor final de x que es la edad de la madre es: " + x);

            Console.WriteLine("El valor final de y es la edad del padre es: " + y);

            Console.WriteLine("El valor final de z es la edad del niño es: " + z);

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

        Console.ReadLine();

        }

​

}