Código: Seleccionar todo
#include <iostream>
#include <fstream>
#include <iomanip>
using namespace std;
void gaussjordan (int n, double **A);
int main(){
ifstream label1 ("datos//datos.dat");
ofstream label2 ("salida//salida.dat");
int i, j, k;
cout << "Resolucion de sistemas lineales por el metodo de Gauss-Jordan\n\n";
label2 << "Resolucion de sistemas lineales por el metodo de Gauss-Jordan\n\n";
int n;
label1 >> n;
cout << "Numero de incognitas = " << n;
cout << "\n\n";
double **A;
A = new double *[n+1];
for (j=0; j < n+1; j++)
A[j] = new double [n+1];
cout.setf(ios::fixed);
cout.precision(4);
label2.setf(ios::fixed);
label2.precision(4);
// Introducción de datos a partir del archivo correspondiente
for (i = 0; i < n; i++) {
for (j = 0; j < n + 1; j++) {
label1 >> A[i][j];
}
}
// Aplicación del método de Gauss-Jordan
gaussjordan (n, A);
// Impresión de valores de las incógnitas
for (i = 0; i < n; i++) {
cout << "X[" << i+1 << "] = " << setw(10) << setiosflags(ios::right) << A[i][n] << " ";
label2 << "X[" << i+1 << "] = " << setw(10) << setiosflags(ios::right) << A[i][n] << " ";
cout << endl;
label2 << endl;
}
cout << endl;
return 0;
}
void gaussjordan (int n, double **A) {
// Aplicación del método de Gauss-Jordan
int i, j, k;
// Normalización de la fila pivote
double factor;
for (k = 0; k < n; k++) {
factor = A[k][k];
for (j = k; j < n + 1; j++) {
A[k][j] = A[k][j]/factor;
}
// Eliminación de variables por fila
for (i = 0; i < n + 1; i++) {
if (i != k) {
factor = A[i][k];
for (j = 0; j < n + 1; j++) {
A[i][j] = A[i][j] - factor * A[k][j];
}
}
}
}
}
Código: Seleccionar todo
5
4 2 3 -1 1 16
1 1 1 1 1 3
1 -1 1 -1 1 1
1 1 2 1 1 6
1 1 -1 1 2 -5 g++ gaussjordan.c++ -o gaussjordan
y ejecuten con:
./gaussjordan
El resultado será:
Código: Seleccionar todo
Resolucion de sistemas lineales por el metodo de Gauss-Jordan
Numero de incognitas = 5
X[1] = 1.0000
X[2] = 2.0000
X[3] = 3.0000
X[4] = -1.0000
X[5] = -2.0000 
. Con relación a los de los punteros dobles esa es la manera como lo implementaba un físico del cual adapté esa forma de programar en C++. Todavía no me ha dado errores.