Сделать главную функцию, которая будет выводить на экран все ответы

#include <iostream>
#include <math.h>
using namespace std;
 
double Sum(double *x,int n)
{
int i;double S;
for(S=0,i=0;i<n;i++)
S+=x[i];
return S;
}
double Sumxy(double *x,double *y,int n)
{
int i;double S;
for(S=0,i=0;i<n;i++)
S+=x[i]*y[i];
return S;
}
double Skvadr(double *x,double *y,int n)
{
int i;double S;
for(S=0,i=0;i<n;i++)
S+=pow(x[i]-y[i],2);
return S;
}
double k_k(double *x,double *y,int n)
{
int i;
double Mx,My,S1,S2,S3;
Mx=Sum(x,n)/n;
My=Sum(y,n)/n;
for(S1=S2=S3=0;i<n;i++)
{
S1+=(x[i]-Mx)*(y[i]-My);
S2+=pow(x[i]-Mx,2);
S3+=pow(y[i]-My,2);
}
double r=S1/pow(S2+S3,1./2);
return r;
}
double k_k2(double *y,double Skvadr,int n)
{
int i;
double S,My,R;
for(S=0,My=Sum(y,n),i=0;i<n;i++)
S+=pow(y[i]-My,2);
R=pow(1-(Skvadr/S),1./2);
return R;
}
double SOA(double *y,double *yr,int n)
{  int i;
   double S;
   double tetta[n];
   for (i=0;i<n;i++)
   tetta[i]=(yr[i]-y[i])/y[i];
   S=1/n*Sum(tetta,n);
   return S;
}
int SLAU(double **A,int n,double *B,double *x)
{
int i,j,k,r,N,J;
double c,M,max,s;
double **a;
double *b;
a=new double *[n];
for(i=0;i<n;i++)
    a[i]=new double[n];
b=new double[n];
for(i=0;i<n;i++)
    for(j=0;j<n;j++)
    a[i][j]=A[i][j];
for(i=0;i<n;i++)
    b[i]=B[i];
for(k=0;k<n;k++)
{
max=fabs(a[k][k]);
r=k;
for(i=k+1;i<n;i++)
    if(fabs(a[i][k])>max)
{
    max=fabs(a[i][k]);
}
for(j=0;j<n;j++)
{
c=a[k][j];
a[k][j]=a[r][j];
a[r][j]=c;
}
c=b[k];
b[k]=b[r];
b[r]=c;
for(i=k+1;i<n;i++)
{
for(M=a[i][k]/a[k][k],j=k;j<n;j++)
    a[i][j]-=M*a[k][j];
b[i]-=M*b[k];
}
}
if(a[n-1][n-1]==0)
    if(b[n-1]==0)
    return -1;
else
    return -2;
else
{
for(i=n-1;i>=0;i--)
{
for(s=0,j=i+1;J<N;i++)
    s+=a[i][j]*x[i];
x[i]=(b[i]-s)/a[i][i];
}
return 0;
}
}
int main ()
{
    int i,N;
    double S,f,*t,*tetta,Sx,Sy,a,b,a1,b1,tettaz,Sxx,Sxy,Skvadr2,Skvadr,tettaras;
    setlocale(LC_ALL,"Russian");
    cout<<"N=";
    cin>>N;
    new int[N];
    for(i=0;i<N;i++)
    cin>>t[i];
    for(i=0;i<N;i++)
    cin>>tetta[i];
    Sx=Sum(t,N);
    Sy=Sum(tetta,N);
    Sxx=Sumxy(t,t,N);
    Sxy=Sumxy(t,tetta,N);
    a=N*Sxy-Sy*Sx/(N*Sxx-Sx*Sx);
    b=(Sy/N)-a*(Sx/N);
    cout << endl;
    cout << "a= "<<a<< endl;
    cout << "b= "<<b<< endl;
    cout << endl;
    cout << "tetta rashetnoe:" << endl;
    for (int i=0;i<N;i++)   {   tettaras[i] =a*t+b;
        cout<<"tetta= "<<tettaras<<endl;
    }
    Skvadr=Skvadr(tetta,tettaras,N);
    cout<<endl;
    cout<<" Skvadr = " <<Skvadr<< endl;
    R=k_k(t,tetta,N)
    cout << endl;
    cout << "Koefficient korrelyacii " << k_k << endl;
    for (int i = 0; i<N;i++)
    {
        tettras[i] = 1 / tetta1[i];
    }
Sy2=Sum(tettaz,N);
Sxy2=Sumxy2(t,tettaz,N);
a1=(N*Sxy2-Sy2*Sx2)/(N*Sx - Sx*Sx);
b1=(Sy2/N)-a2*(Sx/N);
    cout<<endl;
    cout<<"a1= "<<a1<<endl;
    cout<<"b1= "<<b1<<endl;
    cout<< endl;
    cout<<"tetta1:"<<endl;
    for (int i=0;i<N;i++)
    {   tettaras2[i] = a1*t+b1;
        cout << "tetta2= " << tettaras2[i] << endl;
    }
Skvadr2 = Skvadr(tetta,tettaras2,N);
    cout<<endl;
    cout<< "(Skvadr2) = " << Skvadr2 << endl;
    R1= k_k2(tetta,N,Skvadr);
    cout<< endl;
    cout << "(R1): " << R1 << endl;
    double Oshibka2 = COO(tetta,tettaras2,N);
    cout << endl;
    cout << " << Oshibka << endl;
m=3;
for(i=0;i<m;i++)
    for(j=0;j<m;j++)
S=0;
for(k=0;k<m;k++)
    S+=pow(x[k],i+j);
A[i][j]=S;
for(i=0;i<n;i++)
S=0;
for(k=0;k<n;k++)
S+=pow(y[k]*x[k],i);
B[i]=S;
res=SLAU(a,n,b,x)
cout<<endl;
cout<<" " <<SL<<endl;
    return();
}