Hola agradecería su ayuda, aun soy algo nueva en esto. Hasta ahora hice esto pero solo me arroja la ultima raíz, los demás resultados los marca como inf o -inf, no importa cuantas ecuaciones meta solo me arroja el ultimo resultado, todo los demás los marca como inf, seria de gran ayuda para mi si pudieran ayudarme a solucionar el problema porque ya llevo un rato en esto y no me queda.
#include <iomanip>
#include <iostream>
using namespace std;
int main() {
char salir;
int n, i, j, k;
float factor, suma;
float a[30][30], b[30], x[30];
cout << setprecision(4);
cout << fixed;
cout << "\nELIMINACION GAUSSIANA"
<< "\n"; // Eliminacion gaussiana
printf("\n");
do {
cout << "Numero de ecuaciones: ";
cin >> n;
for (i = 0; i < n; i++) {
for (j = 0; j < n; j++) {
printf("Coef. (%i,%i) ", i + 1, j + 1);
scanf("%f", &a[i][j]);
}
printf("Termino independiente Ec. %i: ", i + 1);
scanf("%f", &b[i]);
}
cout << "\nMATRIZ"
<< "\n"
<< endl;
for (i = 0; i < n; i++) {
for (j = 0; j < n; j++) {
cout << " " << a[i][j];
}
cout << " " << b[i];
cout << endl;
}
for (k = 0; k < n - 1; k++) {
for (i = k + 1; i < n; i++) {
factor = a[i][k] / a[k][k];
for (j = k; j < n; j++) {
a[i][j] = a[i][j] - factor * a[k][j];
}
b[i] = b[i] - factor * b[k];
}
}
x[n - 1] = b[n - 1] / a[n - 1][n - 1];
for (i = n - 2; i >= 0; i--) {
suma = b[i];
for (j = i + 1; j < n; j++) {
suma = suma - a[i][j] * x[j];
}
x[i] = suma / a[i][j];
}
cout << "\nMATRIZ ESCALONADA\n";
for (i = 0; i < n; i++) {
for (j = 0; j < n; j++) {
cout << " " << a[i][j];
}
cout << " " << b[i] << "\n";
}
cout << "\nSOLUCION\n";
for (i = 0; i < n; i++) {
cout << x[i] << "\n";
}
cout << "\nDeseas salir s/n: ";
cin >> salir;
} while (salir != 's' && salir != 'S');
return 0;
}