Множественная регрессия

void LinReg(double x[], double y[], const int data,
                    double& slope, double& intercept){
    double sum_x=0, sum_x2=0, sum_y=0, sum_xy=0;
    for(int i=0;i<data;i++){
        sum_x+=x[i];
        sum_y+=y[i];
        sum_x2+=x[i]*x[i];
        sum_xy+=x[i]*y[i];
    }
    slope=(sum_x*sum_y - data*sum_xy)
    /(sum_x*sum_x - data*sum_x2 );
    intercept=(sum_y-slope*sum_x)/data; 
}

#include "stdio.h"
#include "math.h"
#include "stdlib.h"
#include "iostream"
int chlk(double a[],int n,int m,double d[])
{
    int i,j,k,u,v ;
    a[0]=sqrt(a[0]);
    for(j=1;j<=n-1;j++)a[j]=a[j]/a[0];
    for(i=1;i<=n-1;i++)
    {
        u=i*n+i ;
        for(j=1;j<=i;j++)
        {
            v=(j-1)*n+i ;
            a[u]=a[u]-a[v]*a[v];
        }
        
        a[u]=sqrt(a[u]);
        if(i!=(n-1))
        {
            for(j=i+1;j<=n-1;j++)
            {
                v=i*n+j ;
                for(k=1;k<=i;k++)
                    a[v]=a[v]-a[(k-1)*n+i]*a[(k-1)*n+j];
                    a[v]=a[v]/a[u];
            }
            
        }
        
    }
    for(j=0;j<=m-1;j++)
    {
        d[j]=d[j]/a[0];
        for(i=1;i<=n-1;i++)
        {
            u=i*n+i ;
            v=i*m+j ;
            for(k=1;k<=i;k++)
                d[v]=d[v]-a[(k-1)*n+i]*d[(k-1)*m+j];
            d[v]=d[v]/a[u];
        }
    }
    for(j=0;j<=m-1;j++)
    {
        u=(n-1)*m+j ;
        d[u]=d[u]/a[n*n-1];
        for(k=n-1;k>=1;k--)
        {
            u=(k-1)*m+j ;
            for(i=k;i<=n-1;i++)
            {
                v=(k-1)*n+i ;
                d[u]=d[u]-a[v]*d[i*m+j];
            }
            v=(k-1)*n+k-1 ;
            d[u]=d[u]/a[v];
        }
    }
    return(2);
}
 
void sqt2(double x[],double y[],int m,int n,double a[])
{
    int i,j,k,mm ;
    double p,*b ;
    b=(double *)malloc((m+1)*(m+1)*sizeof(double));
    mm=m+1 ;
    b[mm*mm-1]=n ;
    for(j=0;j<=m-1;j++)
    {
        p=0.0 ;
        for(i=0;i<=n-1;i++)
            p=p+x[j*n+i];
        b[m*mm+j]=p ;
        b[j*mm+m]=p ;
    }
    for(i=0;i<=m-1;i++)
        for(j=i;j<=m-1;j++)
        {
            p=0.0 ;
            for(k=0;k<=n-1;k++)
                p=p+x[i*n+k]*x[j*n+k];
            b[j*mm+i]=p ;
            b[i*mm+j]=p ;
        }
    a[m]=0.0 ;
    for(i=0;i<=n-1;i++)
        a[m]=a[m]+y[i];
    for(i=0;i<=m-1;i++)
    {
        a[i]=0.0 ;
        for(j=0;j<=n-1;j++)
            a[i]=a[i]+x[i*n+j]*y[j];
    }
    chlk(b,mm,1,a);
    
    free(b);
    return ;
}
                              
int main()
{
    double a[4];
    double x[3][5]={ {548,460,608,642,836},{7,6,7,7,8},{2,2,2,3,3,}};
    double y[5]={208.5,181.5,223.5,140,250};
    sqt2(&x[0][0],y,3,5,a);
    std::cout<<"y = "<<a[0]<<"*x1"<<a[1]<<"*x2"<<a[2]<<"*x3"<<"+"<<a[3];
}