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

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 pivote(float A[100][100],float B[100],int n){
for (int i=0;i<=n-2;i++){
float piv=fabs(A[i][i]);
int ipiv=i;
for (int j=i;j<n;j++){
if(fabs(A[j][i])>piv){
ipiv=j;
piv=fabs(A[j][i]);}
}

for (int j=i;j<=n-1;j++) {
float temp=A[i][j];
A[i][j]=A[ipiv][j];
A[ipiv][j]=temp;
} 
float temp = B[i];
B[i]=B[ipiv];
B[ipiv]=temp;

int per[i];
int itemp=per[i];          
per[i]=per[ipiv];
per[ipiv]=itemp;

for (int j=i+1; j<n; j++){
float   piv=A[j][i]/A[i][i];
for (int k=i; k<n; k++){    
A[j][k]=A[j][k]- piv*A[i][k];}
B[j]=B[j]- piv*B[i];}
}
}

void trian(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=Bn");
printf("ttALGORITMO ESTRATEGIA DE PIVOTAMIENTOn");
printf("ttIngresar numero de incognitas:");
scanf("%d",&n);
printf("Ingresar A:n");
llenarm(A,n);
printf("Ingresar Bn");
llenarv(B,n);
printf("nEL SISTEMA INGRESADO ES:n");
mostrar(A,B,n);
pivote(A,B,n);
printf("n EL SISTEMA QUEDA:n");
mostrar(A,B,n);
printf("n LA SOLUCION DEL SISTEMA ES:n");
trian(A,B,n);
getch();
getch();  

}