Алгоритм к программе (интерполяция методом Лагранжа)

#include <iostream.h>
#include <DBCtrls.hpp>
#include <DBTables.hpp>
#include <ExtCtrls.hpp>
#include <math.h>
#include <conio.h>
 
double divided_differences(const double* x, const double* y, const int xstart, const int xend)
{
    double dn = x[xend] - x[xstart];
    if (dn == 0) return 0;
    if(xend - xstart == 1)
    {
        return (y[xend] - y[xstart]) / dn;
    } else if(xend - xstart > 1) {
        return (divided_differences(x, y, xstart + 1, xend) - divided_differences(x, y, xstart, xend - 1)) / dn;
    } else {
        return 0;
    }
}
 
 
double newtone_polynomial(const double* x, const double* y, const double xk, const int degree)
{
    if (degree < 1) return 0;
    int i, k; double P = 0;
    for(i = 1; i <= degree; i++)
    {
        double t = 1;
        for (k = 0; k < i; k++) t *= xk - x[k];
        P += divided_differences(x, y, 0, i) * t;
    }
    return y[0] + P;
}
 
int main(int argc, char* argv[])
{
    cout << "Vedit kol-vo uzlov: "; int n; cin >> n; int i;
    double* y = new double[n]; double* x = new double[n];
    cout << "Vedit " << n << " uzlov interpol\n";
    for(i = 0; i < n; i++) cin >> x[i];
    cout << "Vedit " << n << " znach funkc\n";
    for(i = 0; i < n; i++) cin >> y[i];
 
    cout << "\nVedit x: "; double xr; cin >> xr;
 
    cout << newtone_polynomial(x, y, xr, n - 1) << endl;
    delete []x; delete []y;
     getch( ) ;
    return 0;
}