En este programa se intenta realizar el procedimiento del metodo de Gauss para resolver una matriz. La problematica se da cuando empiezan las iteraciones para convertir en cero los valores de la parte inferior de la matriz el programa calcula un valor, pero al momento de imprimirlo imprime uno totalmente diferente. Dicho valor erroneo es el que se toma para la siguiente iteración como valor verdadero, es decir, que esto termina afectando el resultado final de la matriz.
En este programa se intenta realizar el procedimiento del metodo de Gauss
para resolver una matriz. La problematica se da cuando empiezan las
iteraciones para convertir en cero los valores de la parte inferior de la
matriz el programa calcula un valor, pero al momento de imprimirlo imprime
uno totalmente diferente. Dicho valor erroneo es el que se toma para la
siguiente iteracion como valor verdadero, es decir, que esto termina
afectando el resultado final de la matriz
#include<iostream>
#include<windows.h>
using namespace std;
int main(){
int u=0,N,i,j;
float aux,aux1,V[30][30]={0},C[30];
cout<<"DIMENSION: ";
cin>>N;
system("cls");
for(i=0;i<N;i++){
cout<<"INGRESE VALORES MATRIZ:"<<endl;
for(j=0;j<N;j++){
cout<<i+1<<j+1<<": ";
cin>>V[i][j];
}
cout<<"INGRESE CONSTANTE: ";
cin>>C[i];
system("cls");
}
system("cls");
for(i=0;i<N;i++){
for(j=0;j<N;j++){
cout<<V[i][j]<<"\t";
}
cout<<"|\t"<<C[i];
cout<<endl;
}
system("pause");
do{
system("cls");
// for(i=0;i<N;i++){
// for(j=0;j<N;j++){
// if(i==j&&V[i][j]==0){
// cout<<"LA MATRIZ NO TIENE SOLUCION";
// return 0;
// }
// }
// }
aux1=0;
for(i=0;i<N;i++){
for(j=0;j<N;j++){
if(i==u){
aux=V[u][u];
V[i][j]=(V[i][j]/aux);
C[i]=(C[i]/aux);
}
}
}
for(i=0;i<N;i++){
for(j=0;j<N;j++){
cout<<V[i][j]<<"\t";
}
cout<<"|\t"<<C[i];
cout<<endl;
}
system("pause");
system("cls");
//----------------------------------------------------------------------------
// Aplicando el Metodo Gauss para Resolucion de Ecuaciones Lineales hemos
observado
// que los valores que corresponden a las coordenadas, al momento de
imprimir
// no coinciden. Ej: en V[2][2] el valor registrado que debe guardar es -2
// pero al momento de imprimir mu
cout<<"ENTRANDO"<<endl;
for(i=u+1;i<N;i++){
float aux3;
aux3=V[i][u];
for(j=u;j<N;j++){
cout<<i<<endl;
cout<<j<<endl;
cout<<"b: "<<V[i][j]<<endl;
V[i][j]=((-aux3*V[u][j])+V[i][j]);
cout<<"a: "<<V[i][j]<<endl;
cout<<"aux: "<<aux3<<endl;
system("pause");
}
C[i]=((-aux3*C[u])+C[i]);
cout<<C[i]<<endl;
system("pause");
}
//----------------------------------------------------------------------------
for(i=0;i<N;i++){
for(j=0;j<N;j++){
cout<<V[i][j]<<"\t";
}
cout<<"|\t"<<C[i];
cout<<endl;
}
for(i=0;i<N;i++){
for(j=0;j<N;j++){
if(i==j){
aux1=aux1+V[i][j];
}
}
}
u++;
}while(aux1!=N);
return 0;
}