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

int n;
float mat[10][10];
float ans[10];

void swaprows(int r1,int r2)
{
	int i;
	float temprow[10];
	for(i=0;i<=n;i++) temprow[i]=mat[r1][i];
	for(i=0;i<=n;i++) mat[r1][i]=mat[r2][i];
	for(i=0;i<=n;i++) mat[r2][i]=temprow[i];
}

void checkmat(int c)
{
	int i;
	if(mat[c][c]==0)
	{
		if(c==n-1)
		{
			printf("\nThe system of linear equations is inconsistent");
			exit(1);
		}
		for(i=c+1;i<=n-1;i++)
		{
			if(mat[i][c]!=0)
			{
				swaprows(c,i);
				break;
			}
		}
		checkmat(c);
	}
}

void getmat()
{
	int i,j;
	printf("Enter the number of variables: ");
	scanf("%d",&n);
	printf("\nEnter the coefficients along with the constant at the end:\n\n");
	for(i=0;i<n;i++)
	{
		printf("Enter the coeffs of eqn %d : ",i+1);
		for(j=0;j<n+1;j++)
		{
			scanf("%f",&mat[i][j]);
		}
	}
}

void showmat()
{
	int i,j;
	printf("\nThe augmented form of the matrix is:\n");
	for(i=0;i<n;i++)
	{
		for(j=0;j<n+1;j++)
		{
			printf("%.2f\t",mat[i][j]);
		}
		printf("\n");
	}
}

void makezero(int r,int c)
{
	int i;
	float pivot=mat[c][c];
	float factor=mat[r][c]/pivot;
	for(i=0;i<=n;i++)
	{
		mat[r][i]-=factor*mat[c][i];
	}
}

void getans()
{
	int i,j;
	float term=0;
	for(i=n-1;i>=0;i--)
	{
		for(j=i+1;j<n;j++)
		{
			term+=ans[j]*mat[i][j];
		}
		ans[i]=(mat[i][n]-term)/mat[i][i];
		term=0;
	}
}

void showans()
{
	int i;
	printf("\n\nThe solution to the system of linear equations is:\n\n");
	for(i=0;i<n;i++) printf("x%d=\t%f\n",i+1,ans[i]);
}

int main()
{
	int i,j;
	clrscr();
	getmat();
	showmat();
	for(j=0;j<=n-1;j++)
	{
		checkmat(j);
		for(i=j+1;i<n;i++)
		{
			if(i==j) continue;
			makezero(i,j);
			showmat();
			getch();
		}
	}
	showmat();
	getans();
	showans();
	return 0;
}