C PROGRAM FOR GAUSS ELIMINATION METHOD

/* PROGRAM FOR GAUSS ELIMINATION METHOD*/
#include<stdio.h>
#include<conio.h>
void main() {
       float a[20][20],ratio,x[20];
       int i,j,k,n;
       clrscr();
       printf("\n Enter order of matrix:");
       scanf("%d", &n);
       printf("Enter the coefficients & RHS \n");
       for(i=1;i<=n;i++){
          for(j=1;j<=n+1;j++)
          scanf("%f", &a[i][j]);
          printf("\n"); }
       for(k=1;k<=n-1;k++){
          for(i=k+1;i<=n;i++){
                  ratio=a[i][k]/a[k][k];
                  for(j=1;j<=n+1;j++)
                  a[i][j] = a[i][j] - ratio * a[k][j];}}
                  x[n]=a[n][n+1]/a[n][n];                            
       for(k=n-1;k>=1;k--){
         x[k] = a[k][n+1];
         for(j = k+1;j<=n;j++)
         x[k] = x[k]-a[k][j] * x[j];
         x[k]= x[k]/a[k][k];}
       for(i=1;i<=n;i++)
         printf("\n x(%d)=%g",i,x[i]);                
       getch();
       }

/*OUTPUT:-
Enter order of matrix:4
Enter the coefficients & RHS
2       2       1       2       7

-1      2       0       1       -2

-3      1       2       1       -3

-1      0       0       2       0


 x(1)=2.16216
 x(2)=-0.45946
 x(3)=1.43243
 x(4)=1.08108 */