Интерполяция Лагранжа

#include <iostream.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 << "Введите количество узлов: "; int n; cin >> n; int i;
    double* y = new double[n]; double* x = new double[n];
    cout << "Введите " << n << " узлов интерполяции\n";
    for(i = 0; i < n; i++) cin >> x[i];
    cout << "Введите " << n << " значений функции\n";
    for(i = 0; i < n; i++) cin >> y[i];
 
    cout << "\nВведите x: "; double xr; cin >> xr;
 
    cout << newtone_polynomial(x, y, xr, n - 1) << endl;
    delete []x; delete []y;
    return 0;
}

int main(int argc, char* argv[])

void main (){}