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

void LU(float A[100][100],float L[100][100],float U[100][100],int n){

for(int i=0;i<=n-1;i++){

if(i==0){
L[0][0]=1;
for(int j=0;j<=n-1;j++)
{U[0][j]=A[0][j];}
for(int j=i+1;j<=n-1;j++)   
{L[j][0]=A[j][0]/U[0][0];}
}
else{

for(int j=i;j<=n-1;j++)
{
float suma=0;
for(int k=0;k<=i-1;k++)
{suma=suma+U[k][j]*L[i][k];}    
U[i][j]=A[i][j]-suma;   
}
}

for(int j=i;j<=n-1;j++)
{
if(j==i)
{L[j][i]=1;}
else
{
float suma=0;
for(int k=0;k<=i-1;k++)
{suma=suma+U[k][i]*L[j][k];}    
L[j][i]=(A[j][i]-suma)/(U[i][i]);
}

}          
}   
}


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

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

}

void trianL(float L[100][100],float B[100],float y[100],int n){
float suma=0;
y[0]=B[0]/L[0][0];
printf("n%fn",y[0]);
for(int k=1;k<=n-1;k++){
suma=B[k];
for(int j=0;j<=k-1;j++){
suma=suma-L[k][j]*y[j];
}  
printf("%fn",y[k]=suma/L[k][k]);
}
}


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

int main(){

float A[100][100],L[100][100],U[100][100],B[100],y[100];
int n;
printf("tttt  Resolucion de Sistemas AX=B (LU) n");
printf("nIngresar numero de incognitas:");
scanf("%d",&n);
printf("Ingresar Atn");
llenarm(A,n);
printf("Ingresar Btn");
llenarv(B,n);
printf("nEl Sistema es AX=B:n");
mostrar(A,B,n);

LU(A,L,U,n);

printf("nnPresionar cualquier tecla para hallar L y U (LU=A)n");
getch();
getch();
printf("nL=n");
mostrar2(L,n);
printf("nU=n");
mostrar2(U,n);
printf("nnPresionar cualquier tecla para hallar las soluciones:n");
getch();
printf("nLo que bota el algoritmo es:nny=");
trianL(L,B,y,n);
printf("nLas soluciones son: X=n");
trianU(U,y,n);

getch();
getch();
}