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("%d",&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];
}

//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("%d",x[i]);   
              
getch();



}