//Elaborado por Leonardo Argüelles

using System;               

public class Program

{

     public static void Main(string[] args)

        {

            double[] voltaje = {50,55,60,61,65,67,69,70,72,73,74,76};

​

            double[] tiempo = { 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12};

​

            double[,] jacobiana = new double[12, 2];

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

            double piv, factor;

            double x1 = 0, x2 = 0, x3 = 0, t = 0;

            double f1, d1, d2, d3;

​

            f1 = x1 * Math.Log(t) + x2;

            d1 = Math.Log(t);

            d2 = 1;

​

            for (int t1 = 0; t1 < 12; t1++)

            {

                jacobiana[t1, 0] = Math.Log(tiempo[t1]);

                jacobiana[t1, 1] = 1;

            }

​

            for (int ren = 0; ren < 2; ren++) 

            {

                for (int col = 0; col < 2; col++)

                {

                    matriz[ren, col] = 0;

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

                    {

                        matriz[ren, col] += jacobiana[i, ren] * jacobiana[i, col];  // matriz transpuesta por original

                        //matriz[ren, col] += matriz[ren, i] * matriz[col, i];      // matriz original por transpuesta

​

                    }

                }

            }

​

            for (int ren = 0; ren < 2; ren++) 

            {

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

                {

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

                    {

                        matriz[ren, 2] += voltaje[i] * jacobiana[i, ren];

​

​

                    }

                }

            }

​

​

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

            {

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

                {

                    Console.Write(matriz[i, j] + " ");

                }

                Console.WriteLine();

            }

​

            //Eliminación gaussiana

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

            {

                piv = matriz[reng, reng];

                for (int colu = 0; colu < 3; colu++)

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

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

                    if (reng_elimi != reng)

                    {

                        factor = matriz[reng_elimi, reng];

                        for (int colu_elimi = 0; colu_elimi < 3; colu_elimi++)

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

                    }

            }

            Console.WriteLine(" ");

​

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

            {

                Console.Write(matriz[i, 2] + " ");

​

            }

            Console.ReadLine();

        }

    }

​