вывод

//************Interpoljacija Nutona****************
double func(double x)
{
    return pow(2,x)-log(1+x);
}
void njyton(double N,double p_x,vector<double> &x,vector<double> &c)
{
    double P=0,proizv=0;
    int i,j;
   for(i=0;i<=N;i++)
   {
       proizv=1;
       for(j=0;j<i;j++)
           proizv*=(p_x-x[j]);
       
       P+=c[i]*proizv;
   }
   cout<<"f("<<p_x<<")~ "<<P;
   cout<<"\nOshibka= "<<abs(abs(func(p_x))-abs(P))<<"\n";
}
 
void cc(double N,vector<double> &c,vector<double> &x,vector<double> &f)
{
    int i,k,j;
    long double proizv1=1,proizv2=1;
    long double sum=0;
    c.push_back(f[0]);;
    for(i=1;i<=N;i++)
    {
        sum=0;
      proizv1=1;
      proizv2=1;
      for(j=1;j<i;j++)
      {
          proizv1=1;
          for(k=0;k<j;k++)
          proizv1=proizv1*(x[i]-x[k]);
 
          sum+=c[j]*proizv1;
      }
      for(k=0;k<i;k++)
      {
          proizv2=proizv2*(x[i]-x[k]);
      }
      c.push_back((f[i]-c[0]-sum)/proizv2);
    }
    
}
 
void create_x_and_znach_f(double a,double b,double N,vector<double> &x,vector<double> &f)
{
    double prom=0;
    for(int i=0;i<=N;i++)
    {
        prom=a+i*((b-a)/N);
        x.push_back(prom);
        f.push_back(func(prom));
    }
}
 
int main()
{
double a,b,N,p_x;
    vector <double> c_vect,x_vect,f_vect;
 
cout<<"INTERPOLATION NUYTON:\n";
    cout<<"Vvedite 'a': ";
    cin>>a;
    cout<<"Vvedite 'b': ";
    cin>>b;
    cout<<"Vvedite 'N': ";
    cin>>N;
    create_x_and_znach_f(a,b,N,x_vect,f_vect);
    cc(N,c_vect,x_vect,f_vect);
    int f=1;
    while(f)
    {
        cout<<"Vvedite tochky,znachenie f-ii kotoroy hotite opredelit': ";
        cin>>p_x;
        njyton(N,p_x,x_vect,c_vect);
        cout<<"Prododovgytu?(1-tak,0-ni)"<<endl;
        cin>>f;
    };
    
    return 0;
}