Gauss-siedel method to solve system of equations

Discussion in 'C++' started by Nadr, May 16, 2009.

Joined:
Oct 16, 2007
Messages:
165
1
Trophy Points:
0

Introduction

Gauss-siedel method to solve system of equations :undecided

The code

Code:
```#include <stdio.h>
#include <conio.h>
#include <math.h>
#define N 100
#define ACC 0.000001

/************************************************************/
/***** Gauss-siedel method to solve system of equations *****/
/************************************************************/

void main()
{
float coeff[N][N+1]={{4,2,1,11},{3,5,3,22},{2,3,4,20},};
/* Co-efficient inputing variables */
float x[N]={0,0,0};	/* values of the variables i.e. x's */
int n=3;			/* Number of equations */
int i,j;			/*	Loop variables */
float acc;			/* Accuracy determining variable */
char ch;			/* choice inputing variable */
do
{
printf("\nEnter the number of variables\n");
scanf("%d",&n);
for(i=0;i<n;i++)
{
for(j=0;j<n+1;j++)
{
printf("\nEnter %d row %d col element\n",i+1,j+1);
scanf("%f",&arr[i][j]);
}
}
do
{
acc=coeff[0][n];
for(i=0;i<n;i++)
{
x[i]=coeff[i][n];
for(j=0;j<n;j++)
{
if(i==j)
continue;
x[i]=x[i]-coeff[i][j]*x[j];
}
x[i]=x[i]/coeff[i][i];
}
for(i=0;i<n;i++)
acc=acc-coeff[0][i]*x[i];
}while(fabs(acc)>=ACC);
for(i=0;i<n;i++)
printf("%g\t",x[i]);
printf("\nDo you wish to continue[y/n]\n");
fflush(stdin);
scanf("%c",&ch);
}while(ch=='Y' || ch=='Y');
getch();
}```

2. SaswatPadhi~ Б0ЯИ Τ0 С0δЭ ~

Joined:
May 5, 2009
Messages:
1,343
55
Trophy Points:
0
Occupation:
STUDENT !
Location:
Orissa, INDIA

You forgot to remove the background, it contains the default text :p

But anyway, the article was very useful, thanx for that.

Joined:
Jul 12, 2004
Messages:
15,293
365
Trophy Points:
83
Done that for him. Thx for pointing it out