using System;

public class Program

{

    public static void Main()

    {

///DATOS DE ENTRADA(VALORES A DISCRETOS)

double pivote, factor;

int datos = 20;

int incognitas = 3;

int renglones = incognitas, columnas = incognitas+1;

double[] voltaje = { -2.4684, 56.8038, 4.8393, 13.2878, -0.1931, -19.5771,

                     55.4146, 74.4829, 62.0632, 37.3616, 14.1876,

                     72.9088, 81.5205, 120.7811, 156.3607, 141.7361,

                     110.898, 129.7614, 192.8516, 253.4778 };

​

double[] tiempo = { 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16,

                    17, 18, 19, 20};

​

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

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

​

////////////////////MATRIZ JACOBIANA ////////////////////////

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

{

  ///PRIMER COLUMNA (DF/DX)

  matrizJ[i,0]=Math.Pow(tiempo[i],2);

  ///SEGUNDA COLUMNA (DF/DY)

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

  ///TERCER COLUMNA (DF/DZ)

  matrizJ[i,2]= Math.Exp(tiempo[i]/10);

}

 //////////INFORMAR MATRIZ CUADRATICA////////////

//////////MULTIPLICACION J * J^-T /////////////////////

for(int j = 0; j < incognitas; j++)/// renglones de la matriz resultante

  {////j = 1 , a = 0 1 2 

    for (int a = 0; a < incognitas; a++)// columnas de la matriz resultante

    {

      ///CICLO RECORRA RENGLONES DE LA JACOBIANA

      for(int e=0; e<datos; e++)//renglones de la jacobiana

       {

            matriz[j,a] += matrizJ[e,j] * matrizJ[e,a];

       }

    }

  } 

 //////////INFORMAR COLUMNA EXTENDIDA////////////

//////////MULTIPLICACION  J^-T * Vector VOLTAJE /////////////////////

for(int col = 0; col < incognitas; col++)///// RECORRE COLUMNAS

{ /// 1

  for(int ren = 0; ren < datos; ren++) ///// RECORRE RENGLONES

  {

    matriz[col , incognitas ] += matrizJ[ren,col] * voltaje[ren];

  }

}

​

​

///////////INICIA ELIMINACION GAUSSIANA ///////////////////

for(int i=0;i<renglones;i++)// RECORRE RENGLONES (i = renglones)

{

  pivote = matriz[i,i]; //ELEMENTOS DE LA DIAGONAL PRINCIPAL

​

  //CONVERTIR A 1 EL RENGLON DE LA DIAGONAL PRINCIPAL

  for(int j=0; j<columnas; j++)//RECORRE COLUMNAS

  {

      matriz[i,j] = matriz[i,j]/pivote;

  }

​

    //CONVERTIR A CERO EL RESTO DE LOS RENGLONES (RECORRER MATRIZ)

    for(int ren_cero=0; ren_cero<renglones; ren_cero++) 

    {

      if(ren_cero!=i)

      {

          factor = matriz[ren_cero,i]; //0  1  2

          for(int col_cero=0; col_cero<columnas; col_cero++)

          {     

//POSICION - FACTOR * MATRIZ[renglon del piviote , columna actual]

          matriz[ren_cero,col_cero]=matriz[ren_cero,col_cero] - (factor * matriz[i,col_cero]);

          }

      }

    }    

​

}

//IMPRIMIR MATRIZ

for(int i=0;i<renglones;i++) //RECORRE RENGLONES

{

​

      Console.WriteLine("VALOR INCOGNITA " + (i+1) + " : " + matriz[i,incognitas]+" ");

}

    }

}