Апроксимация функции МНК

#include "stdafx.h"
#include "stdio.h"
#include "math.h"
#define N 6 
#define M 4
main()
{ 
    double x[N]={0.115,0.120,0.125,0.130,0.135,0.140},
        y[N]={8.65729,8.29329,7.95829 ,7.64893,7.36235,7.09613};
    double k,h[M],b[M],a[M][M];
    int i,j,s;
        for (i=0;i<M;i++)
        {
            for(j=0;j<M;j++)
            {
                a[i][j]=0;
                    for (s=0;s<M;s++)
                    {a[i][j]=a [i][j]+pow(x[s],i);}}
                        b[i]=0;
                            for (s=0;s<M;s++)
                            {b[i]=b[i]+y[s]*pow(x[s], i+j);
                            }
                    }
            for(i=0;i<M;j++)
            {for(j=0;j<M;j++)
            printf(" %f " , a[i][j]);
            printf("\n"); }
            printf("\n");
            for (i=0;i<M;i++)
                printf(" %f ",b[i]);
            printf("\n");
 
            for (s=1; s<M-1;s++)
            {for (i=s+1; i<M; i++)
            {k=a[i][s]/a[s][s];
            b[i]=b[i]-k*b[s];
            for( j=s+1;j<M;j++)
            {a[i][j]=a[i][j]-k*a[s][j];}}}
 
            for (i=0;i<M;i++)
            {for(j=0;j<M;j++)
            printf(" %f ", a[i][j]);
            printf("\n");}
            printf("\n");
            for(i=0;i<M;i++)
                printf(" %f/n ",b[i]);
            printf("\n");
 
            h[M]=b[M]/a[M][M];
            for (i=M-1;i>1;i--)
            {k=b[i];
            for(j=i+1;i<N;i++)
            {k=k-a[i][j]*h[i];}
            h[i]=k/a[i][i];}
 
            for (i=0; i<M; i++)
                printf("x%d=f\n",i,h[i]);
                    return 0;
}