practica7_1.cpp

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

/*PACTICA 7 EJERCICIOS 1
Escribir una funci髇 que tome como par醡etro una matriz cuadrada y que triangule
la matriz seg鷑 el m閠odo de Gauss-Jordan sin pivoteo.
*/

#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();
}