Metoda najmniejszych kwadratow

#include <iostream>
#include <cmath>
using namespace std;

int main(int argc, char** argv) {
    int n;
    cout<<"Podaj n"<<endl;
    cin>>n;
    float tab[2][n];
    for(int i=0;i<n;i++)
    {

        cout<<"Podaj x"<<i<<endl;
        cin>>tab[0][i];
        cout<<"Podaj f"<<i<<endl;
        cin>>tab[1][i];

    }

    float sumy[4]={0,0,0,0};
    for(int i=0;i<n;i++)
    {

        sumy[0]+=tab[0][i];
        sumy[1]+=tab[1][i];
        sumy[2]+=pow(tab[0][i],2);
        sumy[3]+=tab[0][i]*tab[1][i];
    }

     double W,W1,W2;
     W= (n*sumy[2])-(sumy[0]*sumy[0]);
     W1= (sumy[1]*sumy[2])-(sumy[3]*sumy[0]);
     W2=(n*sumy[3])-(sumy[0]*sumy[1]);
     cout<<W<<" "<< W1<< " "<< W2 <<" "<<endl;
     double a,b;
     a=W1/W;
     b=W2/W;
     cout<<"g(x)="<<b<<"x+("<<a<<")"<<endl;
}