Программирование помогите пж C++

А что тут нужно делать? Откуда берутся массивы узловых значений? Как поступать с интерполяционным полиномом?
Можно попробовать вводить разные узловые значения координат, упорядочивать их по возрастанию их абсцисс и печатать таблицу равномерной табуляции интерполяционного полинома.
#include <iostream>
#include <iomanip>
using namespace std;
double Lagrange(double x,
double *X, double *Y, int N)
{
double sum = 0., product;
int i, j;
for (i = 0; i <= N; i++)
{
product = 1.;
for (j = 0; j <= N; j++)
if (j == i) continue; else
product *= (x - X[j]) / (X[i] - X[j]);
sum += product * Y[i];
}
return sum;
}
int main()
{
int i, j, n;
double a, h, x, u, v;
cout << "n » ";
cin >> n;
double *X = new double [n + 1],
*Y = new double [n + 1];
cout << "x y(x): " << endl;
for (i = 0; i <= n; i++) cin >> X[i] >> Y[i];
for (i = 1; i <= n; i++) for (j = n; j >= i; j--)
if (X[j - 1] > X[j])
{
u = X[j - 1];
v = Y[j - 1];
X[j - 1] = X[j];
Y[j - 1] = Y[j];
X[j] = u;
Y[j] = v;
}
for (i = 0; i <= n; i++)
cout << setprecision(10) << fixed
<< setw(15) << X[i] << setprecision(15)
<< setw(20) << Y[i] << endl;
cout << endl;
h = (X[n] - X[0]) / n;
for (i = 0; i <= n; i++)
{
x = X[0] + i * h;
cout << setprecision(10) << fixed
<< setw(15) << x << setprecision(15)
<< setw(20) << Lagrange(x, X, Y, n) << endl;
}
}

Я завтра в колетж приду что вы делаете скиньте ответ