Апроксимация линейной функции С++

#include "iostream.h"
#include "windows.h"
 
int main(int argc, char* argv[])
{
 
    double *X, *Y;
    int i,n;
 
    double SX,SY,SX2,SXY,Det,Dk,Db,k,b;
 
    SetConsoleCP(1251);
        SetConsoleOutputCP(1251);
    
    cout << "Введите число точек: ";
 
    cin >> n;
 
    X=new double[n];
    Y=new double[n];
 
    cout << "Вводим массив X" << endl;
 
    for (i=0; i < n; i++)
    {
        cout << "X[" << i << "]=";
        cin >> X[i];
    }
 
    cout << "Вводим массив Y" << endl;
 
    for (i=0; i < n; i++)
    {
        cout << "Y[" << i << "]=";
        cin >> Y[i];
    }
 
    SX=0;
    SY=0;
    SX2=0;
    SXY=0;
 
    for (i=0; i < n; i++)
    {
        SX=SX+X[i];
        SY=SY+Y[i];
        SX2=SX2+X[i]*X[i];
        SXY=SXY+X[i]*Y[i];
    }
 
    Det=SX2 * ((double)n)-SX*SX;
    Dk=SXY *  ((double) n)-SX*SY;
    Db=SX2*SY-SX*SXY;
 
    k=Dk/Det;
    b=Db/Det;
 
    cout << "k=" << k << endl;
    cout << "b=" << b << endl;
 
    delete [] X;
    delete [] Y;
 
    return 0;
}

#include "iostream.h"
#include "windows.h"
 
int main(int argc, char* argv[])
{
 
    double *X, *Y; // указатели на динамические массивы X и Y
    int i,n;
 
    double SX,SY,SX2,SXY,Det,Dk,Db,k,b;
 
    SetConsoleCP(1251);               // установка русских букв в консоли
    SetConsoleOutputCP(1251);
    
    cout << "Введите число точек: ";
 
    cin >> n;  // Ввод размеров массивов
 
    X=new double[n];   // создание массива X
    Y=new double[n];    // создание массива Y
 
    cout << "Вводим массив X" << endl;
 
    for (i=0; i < n; i++)
    {
        cout << "X[" << i << "]=";
        cin >> X[i];
    }
 
    cout << "Вводим массив Y" << endl;
 
    for (i=0; i < n; i++)
    {
        cout << "Y[" << i << "]=";
        cin >> Y[i];
    }
 
    SX=0;   // Здесь будет Сумма X
    SY=0;   // Здесь будет Сумма Y
    SX2=0; // Здесь будет Сумма X*X
    SXY=0; // Здесь будет Сумма X*Y
 
 
    for (i=0; i < n; i++)
    {
        SX=SX+X[i];           // подсчет Суммы X
        SY=SY+Y[i];           // подсчет Суммы Н
        SX2=SX2+X[i]*X[i]; // подсчет Суммы X*X
        SXY=SXY+X[i]*Y[i]; // подсчет Суммы X*Y
    }
 
    Det=SX2 * ((double)n)-SX*SX;  // главный определитель системы лин. уравнений
    Dk=SXY *  ((double) n)-SX*SY; // определитель для k
    Db=SX2*SY-SX*SXY;               // определитель для b
 
    k=Dk/Det;
    b=Db/Det;
 
    cout << "k=" << k << endl;
    cout << "b=" << b << endl;
 
    delete [] X; // удаление массивов
    delete [] Y;
 
    return 0;
}