ECUACIONES-LINEALES-POR-GAUSS

[caption id="attachment_1102" align="aligncenter" width="671"] Ecuaciones lineales por Gauss[/caption]

ESTE PROGRAMA DESARROLLARA CUALQUIER SISTEMA DE ECUACIONES POR EL METODO DE GAUSS QUE LO MOSTRARA TODO DESARROLADO POR PANTALLA CON SU RESPECTIVA RESPUESTA DEL VALOR DE CADA UNA DE LAS VARIABLES AL FINAL DEL PROGRAMA



[cpp]

# include < stdio.h >
# include < conio.h >
# include < stdlib.h >

main()
{
float s[100][100],t[100][100],jt[100][100],k,rt[100][100];
int m,n,pp,nn,cc,bb,a,i,j,l,tj,vv,zz;
char f[100];

printf("\n\t PROGRAMA QUE RESUELVE SISTEMA DE ECUACIONES POR GAUSS");
printf("\n\n\n");
system("color 04f");

printf("DEBE HABER EL MISMO NUMERO DE ECUACIONES CON EL MISMO NUMERO DE INCOGNITAS: \n\n");
printf("CANTIDAD DE ECUACIONES: ");
scanf("%d",&l);

pp=1;

for(m=1;m< =l;m++)
{
printf("\n ECUACION %d: \n",m);

nn=1;

for(n=1;n< =l+1;n++) //se inserta cada una de las ecuaciones que conforma la matriz
{
printf("[%d][%d]: ",m,n);
scanf("%f",&s[m][n]);

jt[pp][nn]=s[m][n];
nn++;
}
pp++;
printf("\n");
}

printf("MATRIZ= \n");
for(m=1;m< =l;m++)
{
printf("\t");
for(n=1;n< =l+1;n++) //se inserta cada una de las ecuaciones que conforma la matriz
{
printf("%.2f ",s[m][n]);
printf("\t");
}
printf("\n");
}

printf("\n");

cc=1;
m=1;
printf("PRESIONA ENTER PARA EMPEZAR\n\n");
getch();

while(m< =l)
{
tj=m;
for(n=1;n< =l+1;n++)
{
if(m==n)
{
k=s[m][n];
a=m;

nn=n; j=n;

for(pp=1,i=1;pp< =l;pp++,i++)
{
if(pp!=m)
t[i][j]=jt[pp][nn];
else
t[i][j]=1;

}
}
}

for(n=1,bb=1;n< =l+1;n++,bb++)
{
s[m][n]=s[m][n]/k;

if(s[m][n]==-0)
{
s[m][n]=0;
}
rt[cc][bb]=s[m][n];

}

for(i=1;i< =l;i++)
{
if(i!=j)
{
if(t[i][j] >0)
t[i][j]=-1*t[i][j];

else if(t[i][j]< 0)
t[i][j]=-1*t[i][j];

else
t[i][j]=0;
}
}

for(m=1,i=1;m< =l;m++,i++)
{
if(m!=a)
{
for(n=1,bb=1;n< =l+1;n++,bb++)
{
s[m][n]=rt[cc][bb]*t[i][j]+ s[m][n];
}
}

printf("\t");
for(n=1;n< =l+1;n++)
{printf("%.3f",s[m][n]);
printf("\t");
}
printf("\n");

}
printf("");

printf("\n\n");

for(m=1,pp=1;m< =l;m++,pp++)
{
for(n=1,nn=1;n< =l+1;n++,nn++)
{
jt[pp][nn]=s[m][n];
}
}
cc++; m=tj;
m++;
}
n-=1;

printf("\n EL RESULTADO ES: ");
printf("\n\n");

for(m=1,i=65;m< =l;m++,i++) {
printf("\n\t%c= %0.3f",f[i]=i,s[m][n]); }

getch();
}

[/cpp]