Resolución de Matrices
Enviado por Christian Arcos • 15 de Marzo de 2020 • Trabajo • 435 Palabras (2 Páginas) • 98 Visitas
using namespace std;
#include <iostream>
#include<stdio.h>
#include<stdlib.h>
int main()
{
int i,j,k,n;
float a[100][100],b[100],x[100],l[100][100],y[100],u[100][100],sum;
cout<<endl;
cout<< "CALCULADORA DE SISTEMAS DE ECUACIONES POR EL METODO LU "<<endl;
cout<< ""<<endl;
cout<< "Ingresa el numero de Ecuaciones: "<<endl;
cin>>n;
cout<<"Ingrese los coeficientes:"<<endl;
for(i=0;i<n;i++){
for(j=0;j<=n;j++){
if(j<n){
cout<<"a["<<i+1<<"]"<<"["<<j+1<<"]"<<endl;
cin>>a[i][j];
}else {
cout<<"b["<<i+1<<"](Termino Independiente)"<<endl;
cin>>b[i];
cout<<endl;
}
}
}
cout<<endl;
cout<<"La matriz ingresada es: "<<endl;
cout<<endl;
for(i=0;i<n;i++){
for(j=0;j<=cout<<a[i][j]<<" ";
}
}cout<<endl;
}
cout<<endl;
cout<<"La matriz de terminos idependientes ingresada es: "<<endl;
cout<<endl;
for(i=0;i<n;i++){
cout<<b[i]<<endl;
}
for(i=0;i<n;i++)
for(j=0;j<n;j++)
if(i>j)
{
u[i][j]=0;
}
else if(i==j)
{
l[i][j]=1;
}
else
{
l[i][j]=0;
}
for(i=0;i<n;i++)
{
{
for(k=0;k<n;k++)
if(k!=i)
sum=sum+l[i][k]*u[k][j];
u[i][j]=a[i][j]-sum;
}
else
{
for(k=0;k<n;k++)
if(k!=j)
sum=sum+l[i][k]*u[k][j];
l[i][j]=(a[i][j]-sum)/u[j][j];
}
}
}
cout<<"La Matriz L es :"<<endl;
cout<<endl;
cout<<"L ="<<endl;
for(i=0;i<n;i++)
{
for(j=0;j<n;j++)
{
cout<<
...