Programa En c Para Resolver Sistemas De Ecuaciones Con El Método De Eliminación De Gauss

#include<stdio.h>
#include<conio.h>
#include<stdlib.h>
#include "matriz.c"
#define TRUE 1
void sustitucion(double **a, double *b, double *x, int n)
{
int i,j;double suma;
x[n-1] = b[n-1] / a[n-1][n-1];
for(i=n-2;i>=0;i--)
{
suma=0;
for(j=i+1;j<n;j++)suma+=a[i][j] * x[j];
x[i]= (b[i] - suma)/ a[i][i];
}
}


void eliminacion(double **a, double *b, int n)
{
inti,j,k;
double cte;
for(i=0;i<n;i++)
for(j=i+1;j<n;j++)
{
cte = a[j][i] /a[i][i];
for(k=i;k<n;k++)a[j][k] -= cte * a[i][k];
b[j] -= cte *b[i];
}
}
main()
{
int n;
double **a,*b,*x;
system("cls");
do{printf("Dame el numero de ecuaciones: ");
scanf("%d",&n);
if(n>1) break;
printf("El numero de ecuaciones debe sermayor que 2\n");
}while(TRUE);
a=obtenMemoria(n,n);
b=(double *)calloc(n,sizeof(double));
x=(double*)calloc(n,sizeof(double));
printf("Matriz de coeficientes: \n");
leeMatriz(a,n,n);
printf("Matriz original : \n");imprimeMatriz(a,n,n);
printf("Terminos independientes: \n ");
leeVector(b,n);
eliminacion (a,b,n);
printf("Matriz triangulada:\n");
imprimeMatriz(a,n,n);
sustitucion (a,b,x,n);
printf("Solucion del sistema : \n");
imprimeVector(x,n);...