GAUSS JORDON

#include<stdio.h>
int main()
{
    int i,j,k,n;
    float a[50][50],x[50],r;
    printf("Enter the number of equations: ");
    scanf("%d",&n);
    printf("\nEnter co-efficient matrix [A] row wise:\n");
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        {
            printf("\na[%d][%d]= ",i+1,j+1);
            scanf("%f",&a[i][j]);
        }
        printf("\n");
    }
    printf("\n\nEnter column matrix [B] row wise:\n");
    for(i=0,j=n;i<n;i++)
    {
        printf("\nb[%d]= ",i+1);
        scanf("%f",&a[i][j]);
    }
    printf("\n");
    for(k=0;k<n;k++)
    {
        for(i=0;i<n;i++)
        {
            if(i==k)
                continue;
            r=a[i][k]/a[k][k];
            for(j=0;j<=n;j++)
            {
                a[i][j]=a[i][j]-r*a[k][j];
            }
        }
        printf("\nStep %d:-\n\n",k+1);
        for(i=0;i<n;i++)
        {    
            printf("|");
            for(j=0;j<n;j++)
            {
                 printf("%5.2f",a[i][j]);
            }
            printf(" |\t| %5.2f |\n",a[i][j]);
        }
        printf("\n");
    }
    for(i=0;i<n;i++)
    {
        x[i]=a[i][n]/a[i][i];
        a[i][i]=a[i][i]/a[i][i];
    }
    printf("\nStep %d:-\n\n",k+1);
    for(i=0;i<n;i++)
    {    
        printf("|");
        for(j=0;j<n;j++)
        {
             printf("%5.2f",a[i][j]);
        }
        printf(" |\t| %5.2f |\n",x[i]);
    }
    printf("\n");
    printf("\nThe solution is:\n");
    for(i=0;i<n;i++)
    {
       printf("x[%d]=%5.2f\n",i+1,x[i]);
    }
    return 1;
}