GAUSS JODON INVERSE IN C PROGRAMMING

#include<stdio.h>
#include<stdlib.h>
void main()
{
    int i,j,k,n;
    float a[15][15],d,r;
    printf("Enter order of matrix :");
    scanf("%d",&n);
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=n;j++)
        {
            printf("Enter a[%d][%d] : ",i,j);
            scanf("%f",&a[i][j]);
        }
    }
    for(i=1;i<=n;i++)
    {
        for(j=n+1;j<=2*n;j++)
        {
            if(j==n+i)
            {
                a[i][j]=1;
            }
            else
            {
                a[i][j]=0;
            }
        }
    }
    for(j=1;j<=n;j++)
    {
        if(a[j][j]==0)
        {
            printf("Error");
            exit(0);
        }
        for(i=1;i<=n;i++)
        {
            if(i!=j)
            {
                r=(a[i][j])/(a[j][j]);
                for(k=1;k<=2*n;k++)
                {
                    a[i][k]=a[i][k]-r*a[j][k];
                }
            }
        }
    }
    for(i=1;i<=n;i++)
    {
        d=a[i][i];
        for(j=1;j<=2*n;j++)
        {
            a[i][j]=a[i][j]/d;
        }
    }
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=2*n;j++)
        {
            printf("  %f",a[i][j]);
        }
        printf("\n");
    }
}