Gauss-siedel method to solve system of equations

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();
}

 

Hi Nadr,

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

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