using System;

public class Program

{

    public static void Main()

    {

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

            double pivote, factor;

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

​

            x = 1;

            y = 4;

            z=7;

            f1 =4*Math.Pow(x,2)+8*Math.Sin(y)+2*z+3;

            f2 = Math.Exp(x)+Math.Pow(7,-y)+3*Math.Sin(z)-1;

            f3= Math.Pow(y,2)-Math.Exp(x)+2*z;

​

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

            {

            f1 =4*Math.Pow(x,2)+8*Math.Sin(y)+2*z+3;

            f2 = Math.Exp(x)+Math.Pow(7,-y)+3*Math.Sin(z)-1;

            f3= Math.Pow(y,2)-Math.Exp(x)+2*z;

                matriz[0, 0] = 8 * x;

                matriz[0, 1] = 8*Math.Cos(y);

                matriz[0, 2] = 2;

                matriz[0, 3] = f1;

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

                matriz[1, 1] = Math.Log(7)*Math.Pow(7,-y);

                matriz[1, 2] = 3*Math.Cos(z);

                matriz[1, 3] = -f2;

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

                matriz[2, 1] = 6*y;

                matriz[2, 2] = 2;

                matriz[2, 3] = -f3;

​

                //Eliminación Gaussiana //------------------------------

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

                {

                    pivote = matriz[reng, reng];

                    for (int colu = 0; colu < 4; 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];

                            }

                        }

                    }

                }

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

​

                x = x + matriz[0, 3];

                y = y + matriz[1, 3];

                 z=z + matriz[2,3];

            }

​

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

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

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

        Console.WriteLine("EL CODIGO FUE CREADO POR:ULISES ALBERTO RODRIGUEZ ORTEGA");

            Console.ReadLine();

        }

}