SECOND DEGREE CURVE FITTING IN C PROGRAMMING

//second degree curve-fitting
#include<stdio.h>
void main()
{
    int n,i,j,k;
    float a[10][10],x[10],y[10],r;
    float sumx=0,sumx2=0,sumx3=0,sumx4=0,sumxy=0,sumx2y=0,sumy=0;
    printf("Enter no. of data : ");
    scanf("%d",&n);
    for(i=1;i<=n;i++)
    {
        printf("Enter x%d,y%d : ",i,i);
        scanf("%f%f",&x[i],&y[i]);
    }
    for(i=1;i<=n;i++)
    {
        sumx+=x[i];
        sumx2+=x[i]*x[i];
        sumx3+=x[i]*x[i]*x[i];
        sumx4+=x[i]*x[i]*x[i]*x[i];
        sumxy+=x[i]*y[i];
        sumx2y+=x[i]*x[i]*y[i];
        sumy+=y[i];
    }
    a[1][1]=a[2][2]=a[3][3]=sumx2;
    a[1][2]=a[2][3]=sumx;
    a[2][4]=sumxy;
    a[3][4]=sumx2y;
    a[3][1]=sumx4;
    a[2][1]=a[3][2]=sumx3;
    a[1][4]=sumy;
    a[1][3]=n;
    for(j=1;j<=3;j++)
    {
        for(i=1;i<=3;i++)
        {
            if(i!=j)
            {
                r=a[i][j]/a[j][j];

            for(k=1;k<=4;k++)
            {
                a[i][k]=a[i][k]-r*a[j][k];
            }
            }
        }
    }
    for(i=1;i<=3;i++)
    {
        x[i]=a[i][4]/a[i][i];
    }
    printf("a = %f\t b = %f\tc = %f",x[1],x[2],x[3]);
}