algoritmo triangulacion[1]..cpp

Topics Practices: Programming and Numerical Methods Practice 1: Introduction to C Practice 2

//Resolver sistemas de ecuaciones..

#include<stdio.h>
#include<conio.h>

 
 main()
{ 
      int i,j,k, suma,n,h,t;
      float m;
      
      printf("Ingrese el tama駉 de la matriz cuadrada ");
      scanf("%d",&n);
      
      float  a[n][n+1];
      printf("Ingrese los elementos de la matriz  ampliada por fila\n ");
      for(i=0;i<n;++i)
        for(j=0;j<n+1;++j)
        {      
        //printf("Ingrese el elemento a[%d][%d]\n",i+1,j+1);
        scanf("%f",&a[i][j]);
        }


//primero triangulamos con algoritmo de triangulacion de gauss.

  
//int a[n][n+1];//matriz aimpliada n+1 significa que le agrego una columna mas no?

for(j=0;j<=n-2;++j)
    //if(a[j][j]==0)//pivotea..hay que intercambiar filas.
//pivotear;
for(i=j+1;i<n;++i)

{                  
    m = -a[i][j]/a[j][j];
    for(k=j;k<=n;++k)
    a[i][k] = a[i][k]+ m*a[j][k];
}

 for(i=0;i<n;++i)
 {
  for(j=0;j<n+1;++j)
    printf("%.2f ", a[i][j]);
  printf("\n");  
 }

//solucion del sistema

     float x[n];
     for(i=n-1;i>=0;i=i-1)

{     suma=0;
      for(j=i+1;j<n;++j)
           suma=suma+a[i][j]*x[j];
      x[i] = (a[i][n] - suma) / a[i][i];
}
      
      for(i=0;i<n;++i)
       printf("%.2f ",x[i]);   
              
getch();



}