Gauss-Seidel
quien tenga el programa de Gauss-Seidel por favor me lo mandan a mi correo [email protected]............. se los agradeceria...............gracias
Esto lo encontre en un sitio llamado "El rincon del C'
#include <stdio.h>
#include <math.h>
#include <conio.h>
#define L 10
#define A L
main()
{
float M[L][A],X[L],CX[L],C[L];
float a, b, c, d, e, f;
int it,ini,n,x,y,z;
clrscr();
printf("nnnSolucion de ecuaciones simultaneasnnnpor el metodo de Gauss-Seidel");
printf("ntCuantas incognitas tendra el sistema: ");
scanf("%d",&n);
for(x=0;x<n;x++)
{
CX[x]=0;
X[x]=0;
}
for(y=0;y<n;y++)
{
for(x=0;x<n;x++)
{
printf("A[%d][%d] = ",y,x);
scanf("%f",&e);
M[y][x]=e;
}
printf("Y la constante C[%d] = ",y);
scanf("%f",&f);
C[y]=f;
}
for(it=0;it<20;it++)
{
for(y=0;y<n;y++)
{
for(x=0;x<n;x++)
CX[y]-=(M[y][x]*X[x])/M[y][y];
CX[y]+=(C[y]/M[y][y]);
X[y]=CX[y];
}
}
printf("El resultado es ");
for(x=0;x<n;x++)
printf("nX[%d]=%f",x+1,X[x]);
getch();
}
#include <stdio.h>
#include <math.h>
#include <conio.h>
#define L 10
#define A L
main()
{
float M[L][A],X[L],CX[L],C[L];
float a, b, c, d, e, f;
int it,ini,n,x,y,z;
clrscr();
printf("nnnSolucion de ecuaciones simultaneasnnnpor el metodo de Gauss-Seidel");
printf("ntCuantas incognitas tendra el sistema: ");
scanf("%d",&n);
for(x=0;x<n;x++)
{
CX[x]=0;
X[x]=0;
}
for(y=0;y<n;y++)
{
for(x=0;x<n;x++)
{
printf("A[%d][%d] = ",y,x);
scanf("%f",&e);
M[y][x]=e;
}
printf("Y la constante C[%d] = ",y);
scanf("%f",&f);
C[y]=f;
}
for(it=0;it<20;it++)
{
for(y=0;y<n;y++)
{
for(x=0;x<n;x++)
CX[y]-=(M[y][x]*X[x])/M[y][y];
CX[y]+=(C[y]/M[y][y]);
X[y]=CX[y];
}
}
printf("El resultado es ");
for(x=0;x<n;x++)
printf("nX[%d]=%f",x+1,X[x]);
getch();
}
