triangulacion_con_pivoteo.cpp

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

//Triangulacion con pivoteo
#include<stdio.h>
#include<conio.h>
#include<math.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 n,i,j;
 float E;
 printf ("Dame el margen de error");
 scanf("%f",&E);
 printf("ingrese un numero");    
 scanf("%d",&n);
 float A[n][COLS];
 printf("Ingrese los elementos de la matriz");
 for(i=0; i<n;++i)
 {
   for(j=0; j<n; ++j)
   {
     scanf("%f",&A[i][j]);
   }
 }
 printf("\n\n");
 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("\n");
   }
 }
 printf("\n\n");
getch();
}