using System;

​

public class Program {

    public static void Main() {

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

        double[] f = {7.2074, 10.5465, 9.7056, 8.2160, 10.2054, 16.6029, 24.2849, 28.9468, 29.0606, 27.2799};

​

        int datos = f.Length, incognitas = 2;

        int ren = incognitas, col = ren + 1;

        double factor, pivote;

​

        double [,] matrizJacobiana = new double[datos, incognitas];

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

​

        //MATRIZ JACO

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

        {

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

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

        }

​

        //MATRIZ CUADRADA GAUSS = JACOBIANA * JACOBIANA TRANSPUESTA 

        for (int i = 0; i < incognitas; i++) //recorre renglones matriz Gauss

        {

            for (int j = 0; j < incognitas; j++) //recorre columnas matriz Gauss

            {

                for (int k = 0; k < datos; k++) //recorre datos/renglones jacobiana

                {

                    //la j se mueve y la i queda fijo

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

                }

            }

        }

​

        //COLUMNA EXTENDIDA: CADA COLUMNA DE LA JACOBIANA * VECTOR DE DATOS

        for (int i = 0; i < incognitas; i++) //recorre columnas matriz Jacobiana

        {

            for (int j = 0; j < datos; j++) //recorre renglones matriz Jacobiana

            {

                matriz[i,incognitas] +=  matrizJacobiana[j,i] * f[j];

            }

        }

​

        //GAUSS

        // Recorrer matriz

        for (int r = 0; r < ren; r++) //RECORRER RENGLONES r = 0 -> 1

        {

            pivote = matriz[ r , r ] ; 

            for (int c = 0; c < col; c++) //RECORRER COLUMNAS C = 0 

            {

                matriz[r,c] = matriz[r,c] / pivote;   

            }

            //VOLVER A RECORRER LA MATRIZ PARA HACER LAS CONVERSIONES A CERO

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

            { 

                if(r != rCero) //BRINCAR EL RENGLON DEL PIVOTE

                {

                        factor = matriz[rCero, r ];

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

                        {

                            //(VALOR ORIGINAL ) – (RENGLON DEL PIVOTE,C)(FACTOR))\

                            matriz[rCero, cCero] = matriz[rCero, cCero] - (factor * matriz[r, cCero]);

            }

        }

    }

}

​

//IMPRESION GAUSS INCOGNITAS

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

{

    Console.WriteLine("Variable " + (r + 1) + ": " + matriz[r, col - 1]);

}

    }

}

​

​