using System;

​

public class Program

{

    public static void Main()

    {

        double pivote, factor;

        int N = 10; // N es el total de elementos 

        int variables = 2; // variables es el número de variables xi en la ecuación 

        double[, ] jacobiana = new double[N, variables];

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

        //Temperaturas de cada mes

        double[] voltaje =

        {

            7.2074,

            10.5465,

            9.7056,

            8.2160,

            10.2054,

            16.6029,

            24.2849,

            28.9468,

            29.0606,

            27.2799

        };

        double[] tiempo =

        {

            1,

            2,

            3,

            4,

            5,

            6,

            7,

            8,

            9,

            10

        };

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

        {

            jacobiana[i, 0] = tiempo[i];

            jacobiana[i, 1] = Math.Sin(tiempo[i]);

        }

​

        //Multiplicación de matrices

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

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

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

                for (int k = 0; k < N; k = k + 1)

                    matriz[i, j] = matriz[i, j] + jacobiana[k, i] * jacobiana[k, j];

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

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

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

                    matriz[i, variables] = matriz[i, variables] - voltaje[k] * jacobiana[k, i];

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

        //Eliminación Gaussiana

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

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

        {

            pivote = matriz[reng, reng];

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

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

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

                if (reng_elimi != reng)

                {

                    factor = matriz[reng_elimi, reng];

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

                        matriz[reng_elimi, colu_elimi] = matriz[reng_elimi, colu_elimi] - factor * matriz[reng, colu_elimi];

                }

        }

​

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

        //Imprime los valores de las variables

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

        {

            Console.WriteLine("Coeficiente " + (i + 1) + " = " + matriz[i, variables]);

        }

    }

}