MATH Lab- Gauss-Elimination method (Solving system of equation);

#include<stdio.h>

main()
{
    printf("\nEnter the size of matrix:");
    int size,i,j,n;

    scanf("%d",&size);
    float mat[size][size+1],x[size];

    printf("\nEnter the augmented matrix?\n");

    for(i=0;i<size;i++)
    {
        for(j=0;j<size+1;j++)
        {
            scanf("%f",&mat[i][j]);

        }
    }//for-i

    //execution;
    double k;
    for(n=0;n<size-1;n++)
    {
        for(i=n+1;i<size;i++)
        {
            for(j=n;j<size+1;j++)
            {
                if(j==n)
                {
                    k=(mat[i][j]*1.0)/(1.0*mat[n][j]);
                }//if

                mat[i][j]-=k*mat[n][j];
            }//for-j
        }//for-i

    }//for-n

    //termination;
    for(i=size-1;i>=0;i--)
    {
        k=0;
        for(j=i+1;j<size;j++)
        {
            k+=mat[i][j]*x[j];
        }

        x[i]=(mat[i][size]-k)/mat[i][i];
    }//for i

    printf("\nAnswer:-\n");
    for(i=0;i<size;i++)
        printf("x%d=%f, \n",i+1,x[i]);

}//end

OUTPUT:-

-kkmohanta