using System;

public class Program

{

    public static void Main()

    {       

        double x, y, z, f1, f2, f3, error = 0.001;

        double pivote, factor;

        int nv = 3;

        double [,] matriz = new double [nv, nv+1];

        x = 1;

        y = 1;

        z = 1;

        f1 = -x + z * Math.Sin(x + y);

        f2 = -x + Math.Exp(x + z);

        f3 = x*z - y*z - x*y;

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

        {

            f1 = -x + z * Math.Sin(x + y);

            f2 = -x + Math.Exp(x + z);

            f3 = x*z - y*z - x*y;

            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 //------------------------------

            //aquí empieza el método de Gauss

            for (int reng = 0; reng < nv; reng = reng + 1)

            {

                //Selecciona el elemento diagonal como pivote

                pivote = matriz[reng, reng];

​

                //Divide todo el renglón entre el pivote

                for(int colu= 0; colu < nv+1; colu++)

                {

                matriz[reng,colu]= matriz[reng,colu]/pivote;

                }

                //

                for(int reng_elimi = 0; reng_elimi < nv; reng_elimi++)

                {

                    if(reng_elimi != reng)

                    {

                        factor = matriz[reng_elimi,reng];

                        for(int colu_elimi = 0; colu_elimi < nv+1; colu_elimi++)

                        {

                            matriz[reng_elimi, colu_elimi] = matriz[reng_elimi,colu_elimi]-

                                factor*matriz[reng,colu_elimi];

                        }

                    }

                }

            }

        //aquí se muestra el resultado en la pantalla

            for (int reng = 0; reng < nv; reng++)

            {

                for(int colu= 0; colu < nv+1; colu++)

                {

                    Console.Write(matriz[reng,colu] + "       ");

                }

                Console.WriteLine();

            }

            Console.WriteLine();

​

                x = x + matriz[0, 2];

                y = y + matriz[1, 2];

            }

​

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

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

    }

}