Нахождение формулы аппроксимации

double s = 0;
double m = 1;
for (int i = 0; i < n; i++)
{
    s += a[i] * m;
    m *= x;
}

    for (i = 0; i < N + 1; i++)
    {
        y1 += a[i] * pow(x[i], i);
        ShowMessage(FloatToStr(y1));
        Form1->ListBox1->Items->Add(FloatToStr(y1));
        Form1->Series2->AddY(y1);
    }

double a[3] = {1, 2, 3};
int n = 3;
double x = 8;
 
double s = 0;
double m = 1;
for (int i = 0; i < n; i++)
{
    s += a[i] * m;
    m *= x;
}

const int N = 3, M = 4;
double a[N] = {1, 2, 3};
double x[M] = {3, 2, 1, 0}; 
double s[M] = {0};
 
for ( int j = 0, xj = 0, i; j < M; xj = x[++j];) {  
    for ( i = 0; i < N; i++ ) {
        s[j] += a[i] * xj;
        xj *= x[j];
    }
}

for(j = 0; j < K + 1; xj = x[j++])
{
    for(i = 0; i < N; i++)
    {
        func[j] += a[i] * xj;
        xj *= x[i];
    }
 
}
for (int i = 0; i < N; i++) {
Form1->ListBox1->Items->Add(FloatToStr(func[i]));
}

func[0] = a[0];
func[1] = a[0] + a[1] * x[0];
func[2] = a[0] + (a[1] * x[1]) + (a[2] * pow(x[1], 2));
func[3] = a[0] + (a[1] * x[2]) + (a[2] * pow(x[2], 2)) + (a[3] * pow(x[2], 3));
func[4] = a[0] + (a[1] * x[3]) + (a[2] * pow(x[3], 2)) + (a[3] * pow(x[3], 3)) + (a[4] * pow(x[4], 4));

    const int N = 5;
    double a[N] = {1, 2, 3, 4, 5};
    double x[N] = {3, 2, 1, 5};
    double s[N] = {0};
 
    for ( int j = 0, xj = 0, i; j < N; j++ ) {
        for ( i = 1, ( j > 0 ? xj = x[j - 1] : xj ), s[j] = a[0]; i <= j; i++ ) {
            s[j] += a[i] * xj;
            xj *= x[j-1];
        }
    }