GAUSS ELIMINATION METHOD IN C PROGRAMMING

//gauss elimination
#include<stdio.h>
void main()
{
    int i,j,k,n;
    float r,a[20][20],x[20],sum;
    printf("Enter matrix size : ");
    scanf("%d",&n);
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=n+1;j++)
        {
            printf("Enter a%d%d : ",i,j);
            scanf("%f",&a[i][j]);
        }
    }
    for(j=1;j<=n-1;j++)
    {
        for(i=j+1;i<=n;i++)
        {
            r=a[i][j]/a[j][j];
            for(k=1;k<=n+1;k++)
            {
                a[i][k]=a[i][k]-r*a[j][k];
            }
        }
    }
    x[n]=a[n][n+1]/a[n][n];
    for(i=n-1;i>=1;i--)
    {
        sum=0;
        for(j=i+1;j<=n;j++)
        {
            sum=sum+a[i][j]*x[j];
        }
        x[i]=(a[i][n+1]-sum)/a[i][i];
    }
    for(i=1;i<=n;i++)
    {
        printf("x%d = %.4f\t",i,x[i]);
    }
}