ESTE PROGRAMA ES ELAVORADO EN C++
#include<stdio.h>
#include<conio.h>
#include<math.h>
void llenarm(float matriz[100][100],int n){                                        
for(int i=0;i<=n-1;i++){
for(int j=0;j<=n-1;j++){
printf("Ingresar A[%d][%d]=",i+1,j+1);
scanf("%f1.1",&matriz[i][j]);}
}
}   

void mostrar(float matriz[100][100],float B[100],int n){                         
for(int i=0;i<=n-1;i++){
printf("nn");         
for(int j=0;j<=n-1;j++){
printf("t%f",matriz[i][j]);}
printf("t%f",B[i]);
}               
}

void eliminacion(float A[100][100],float B[100],int n)
{
float m=0;
for(int k=0;k<=n-2;k++){
for(int i=k+1;i<=n-1;i++){
m=A[i][k]/A[k][k];
B[i]=B[i]-m*B[k];
for(int j=k;j<=n-1;j++){
A[i][j]=A[i][j]-m*A[k][j];
}
}               
}
}

void llenarv(float B[100],int n){
for(int j=0;j<=n-1;j++){
printf("Ingresar B[%d]=",j+1);
scanf("%f1.1",&B[j]);}

}

void trii(float A[100][100],float B[100],int n){
float x[100];
float suma=0;
x[n-1]=B[n-1]/A[n-1][n-1];
printf("n%fn",x[n-1]);
for(int k=n-2;k>=0;k--){
suma=B[k];
for(int j=k+1;j<=n-1;j++){  
suma=suma-A[k][j]*x[j];
}  
printf("%fn",x[k]=suma/A[k][k]);
}                     

}  
int main(){

float A[100][100],B[100];
int n;
printf("ttResolucion de Sistemas AX=Bnt  nnIngresar numero de incognitas:n");
scanf("%d",&n);
llenarm(A,n);
printf("Ingresar Bn");
llenarv(B,n);
printf("nEl Sistema ingresado es:n");
mostrar(A,B,n);
eliminacion(A,B,n);
printf("nnEl Sistema reducido a triangular es:n");
mostrar(A,B,n);
printf("nn****Las Soluciones del Sistema son:****n");
trii(A,B,n);


getch();
getch();
}