using System;

public class Program

{

    public static void Main()

    {

        Console.WriteLine("RETO 2-C | Xiao Ping");

        Console.WriteLine("Guerreros:");

        Console.WriteLine("Angy y Kaleb");

        Console.WriteLine();

        Console.WriteLine("El primer programa debe obtener el valor cuando x = 6.574, según los siguientes datos:");

        //Arreglos para la dupla de valores x,y

        double[] x = {-5, -3, -0.7, 0.25, 2.1, 6, 7.46, 19.1, 15.5};

        double[] y = {6, 5.3, 1.53, -2.7, 4, 9.1, 2.2, 3.5, 6.2};

        Console.WriteLine();

        Console.WriteLine("[x: " + x[0] + " | " + x[1] + "  | " + x[2] + " | " + x[3] + " | " + x[4] + " | " + x[5] + "   | " +  x[6] + " | " + x[7] + " | " + x[8] + " ]");

        Console.WriteLine("[y: " + y[0] + "  | " + y[1] + " | " + y[2] + " | " + y[3] + " | " + y[4] + "   | " + y[5] + " | " +  y[6] + "  | " + y[7] + "  | " + y[8] + "  ]");

        Console.WriteLine();

        Console.WriteLine("c.El tercer programa deberá arrojar todos los valores intermedios entre x = -5 y x = 32,");

        Console.WriteLine("cada 0.4. Es decir, x = -5, x = -4.6, x = -4.2, x = -3.8, etc. Esto se hará con el método que se desee.");

        Console.WriteLine();

        // se necesita una matriz para crear el gauss jordan.

        double[,] matriz = new double[9,10];

        // variables auxiliares para el pivote y factor parte del algoritmo de Gauss Jordan

        double pivote, factor;

        // La variable x_0 es el valor buscado de la interpolacion. 

        double fx = 0;

        //Se crea la matriz de n ecuaciones para crear nuestro sistema de ecuaciones las primera n-1 columnas 

        for (int i = 0; i < 9; i = i + 1)

        {

            for (int j = 0; j < 9; j = j + 1)

            {

                matriz[i, j] = Math.Pow(x[i],j);

            }

        }

        //Se procede a colocar la ultima columna para tener nuestra matriz aumetada y poder realizar el Gauss Jordan

        for (int i = 0; i < 9; i = i + 1)

        {

            matriz[i, 9] = y[i];

        }       

        //Eliminación Gaussiana

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

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

        {

            pivote = matriz[reng, reng];

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

            {

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

            }

            for (int reng_elimi = 0; reng_elimi < 9; reng_elimi = reng_elimi + 1)

            {

                if (reng_elimi != reng)

                {

                    factor = matriz[reng_elimi, reng];

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

                    {

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

                            - factor * matriz[reng, colu_elimi];

                    }

                }

            }

        }

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

       // Se realiza el calculo de la ecuacion para obtener el valor de y (fx) buscado.

        double [] x_0 = new double [1000];

        int a = 0;

        for (double i = -5; i <= 32; i = i + 0.4)

        {

            x_0[a] = i;

            a++;

        }

        int b = 0;

        while (b < a)

        {

            for (int i = 0; i < 9; i++)

            {

            fx = fx + matriz[i, 9]*Math.Pow(x_0[b],i);

            }

            Console.WriteLine("El valor de F(X) = " + x_0[b] + " es: " + fx);

            b++;

        }

    }

}

​