Gauss Seidel

Descrição

Esconder código-fonte

/*Matematica numerica
Claudir Oliveira
Método de Gauss Seidel
*/
/*
<1º passo-Criação de ponteiros e alocação de ponteiros>
<2º passo-Criar ponteiro para arquivos e escrever aquivos>
############################################################
compilar:gcc SeidelModificado.c -o <name>
*/
#include<stdio.h>
#include<math.h>
#include<malloc.h>
double **a, *b, *x;
int i, j, n;

main()
{
 void gauss_Seidel(), aloc_var();
 int s, m, n;
 s = 0;
 while(s==0)
 {
  printf("\033[036m");
  printf("\033[04m");
  printf("Resolução de SELAs - Gauss-Seidel \n");
  printf("\033[00m");
  printf("\033[036m");
  printf("\n Forneça as dimensões da matriz A: m n = ");
  printf("\033[00m");
  scanf("%d %d", &m, &n);
  aloc_var(m,n);
  gauss_Seidel(m,n);
  printf("\033[036m");
  printf("\033[05m");
  printf("\n Deseja sair do progrma? (1-Sim, 0-Não)");
  printf("\033[00m");
  scanf("%d", &s);
 }
}


void gauss_Seidel(m,n)
 int m, n;
{
 double somat();
 int k, tol;
 printf("\033[036m");
 printf("\n Defina a matriz A[%dx%d]\n",m,n);
 printf("\033[00m");
 for(i=0; i<= m-1; i++)
 {
  for(j=0; j<=n-1; j++)
  {
   printf("\n a[%d,%d]=  ",i+1,j+1);
   scanf("%lf",&a[i][j]);
  }
 }
 printf("\033[036m");
 printf("\n Defina o vetor B(%d)\n",m);
 printf("\033[00m");
 for(i=0; i<=m-1; i++)
 {
   printf("\n b[%d]=  ",i+1);
   scanf("%lf",&b[i]);
 }
 printf("\033[036m");
 printf("\n Defina o vetor x(%d)\n",n);
 printf("\033[00m");
 for(i=0; i<=n-1; i++)
 {
   printf("\n x[%d]=  ",i+1);
   scanf("%lf",&x[i]);
 }
 printf("\033[036m");
 printf("\n Defina o número de iterações = ");
 printf("\033[00m");
 scanf("%d", &tol);
 k = 0;
 while (k < tol)
 {
  for(i=0; i<=n-1; i++)
  {
   x[i] = (b[i] - somat(i,n) + a[i][i]*x[i])/a[i][i];
  }
  k = k+1;
 }
 printf("\033[036m");
 printf("\n O vetor x é = \n");
 printf("\033[00m");
 for (i = 0; i <= n-1; i++)
 {
  printf("%g \n", x[i]);
 }
}

double somat(i,n)
 int i, n;
{
 double soma;
 soma = 0.0;
 for(j=0; j<=n-1; j++)
 {
  soma = soma + a[i][j]*x[j];
 }
 return soma;
}
void aloc_var(m,n)
 int m,n;
{
 int i;
printf("\n Valor de m , n = %d,%d\n",m,n);
a=(double **) malloc((m+1)*sizeof(double *));
b=(double *) malloc((m+1)*sizeof(double));
x=(double *) malloc((n+1)*sizeof(double));
 for(i=0;i < n+1;i++)
 {
  a[i]=(double *)malloc((n+1)*sizeof(double));
 }
}