esinversiblelamatrz.cpp

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

#include<conio.h>
#include<math.h>
#include<stdio.h>
#define cols 100
void triangular(float A[][cols],int n,float E)
{ float m,aux;
int i,j,k,t;
   for(j=0 ; j<n-1 ; ++j)
   { 
     if(fabs(A[j][j])<E)
      { t=j+1;
        while(t<n && fabs(A[t][j])<E)
           ++t;
        if(t==n)
        {
           printf("La matriz no es inversible");
           return;
        }
        for (k=j; k<=n; ++k)
        { 
            aux =A[t][k];
            A[t][k]=A[j][k];
            A[j][k]=aux;
        }        
     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];
     }
   }
}
}               
main()
{int i,j,n;
 float E;
 printf("Ingresa el margende error");
 scanf("%f",&E);
 printf("Ingresa el numero de filas que queres que tenga la matriz");
 scanf("%d",&n);
 float A[n][cols];
 printf("Ingresa los elementos de la matriz extendida \n "); 
 for(i=0 ; i<n ;++i)
 { 
    for(j=0 ; j<=n ; ++j)
    {     
    scanf("%f",A[i][j]);

}
}
 triangular(A,n,E);
 for(i=0 ; i<n ;++i)
 {
  for(j=0 ; j<=n ; ++j)
  {
  printf("%.0f ",A[i][j]);
  if (j==n-1)
  printf("| ");
  if (j==n)
  printf("\n");
}
}
getch();
}