Построение полинома методом Ньютона

В чем у меня ошибка? В узлах интерполирования значения полинома должны совпадать со значениями функции. #c++

#include <iostream>
#include <cmath>
using namespace std;

#define c 1000
#define n 5
double function_1(double l)
{
double m;
m = pow(l,5);
return m;
}

double no(double x, int r, double mas_x[n+1], double mas_y[n+1])
{
int p;
double m;
double l0 = 1;
double l_0[n];
int oi, ok;
int jj = r-2;
int x_i = 2;
if (r==1){
l0 = mas_y[r];
}
else if(r==2){
l0 = (mas_y[2] - mas_y[1])/(mas_x[2] - mas_x[1]);
}
else{
for (oi = 1; oi<r-1; oi++){
l_0[oi] = (mas_y[oi+1] - mas_y[oi])/(mas_x[oi+1] - mas_x[oi]);
}
for (oi = 1;oi<= jj; oi++){
for (ok = 1; ok<=jj; ok++){
l_0[ok] = (l_0[ok+1] - l_0[ok])/(mas_x[ok+x_i] - mas_x[ok]);
}
jj = jj-1;
x_i = x_i+1;
}
l0 = l_0[1];

}
return l0;
}
double poly (double y[n+1], double x1, double mas_x1[n+1]) {
int kl;
double ko = 1;
double f_meaning = no(x1,1,mas_x1,y);
for (kl=2; kl<=n+1; kl++) {
ko = (x1 - mas_x1[kl-1])*ko;
f_meaning = ko*no(x1, kl, mas_x1,y) + f_meaning;
}
return f_meaning;
}
int main() {
const int pj=100;
double a=0;
double b=1;
double x [c];
double f[c];
double i;
int j =1;
int m;
double max[n+1];
double step = (b-a)/n;
double p_y[n+1];
cout<<"x meanings\n";
for (i=a; i<=b {
x[j] = i;
i=i+step;
cout<<x[j]<<"\n";
j++;
};
int r;
cout<<"\n";
cout<<"y meanings\n";
for (r=1; r<j;r++){
f[r] = function_1(x[r]);
cout<<f[r]<<"\n";
}
int pl;
cout<<"\n";
cout<<"polynom meanings\n";
for (pl=1; pl<j;pl++){
cout<<poly(f, x [pl], x);
cout<<"\n";}
system("pause");
return 0;
}

0