Интерполяция заданной функции второй формулой Ньютона (обратная интерполяция)

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
 
double *v;
 
double *Create_Vector(int N);
 
void node(int N, int a, int b);
 
double f(double x);
 
double Polinom(int N, double x);
 
double F3(int i, int n);
 
void Print(int N);
 
void clean(int N);
 
int main()
{
    setlocale(LC_ALL,"RUSSIAN");
    int a,b,N;
    double x,S=0;
    cout << "Функция: y(x) = pow(2,x)" << "\n";
    cout << "Введите границы интервала" << "\n" << "a: "; cin >> a; cout << "b: "; cin >> b;
    cout << "Введите количество узлов: "; cin >> N; cout << "\n";
    v=Create_Vector(N);
    node(N,a,b);
    Print(N);
    cout<<"Введите x: "; cin >> x;
    cout<<"Значение f(x): " << f(x) << "\n";
    cout<<"Значение P(x): " << Polinom(N,x) << "\n";
    clean(N);
    system("pause");
    return 0;
}
 
double *Create_Vector(int N) {
        return (double *)calloc(N,sizeof(double));
}
 
void node(int N, int a, int b) {
    double c=(double)(b-a)/(N-1);
    for(int i=0;i<N;i++)
    {
        v[i]=a+c*i;
    }
}
 
double Polinom(int n, double X) {
    double t=(X-v[n-1])/(v[1]-v[0]);
    double sum=v[n-1];
    double f=1;
    double temp=t;
    for(int i=1; i<n; i++) {
        sum=sum+((temp)/f*i)*F3(i, n-i);
        f=f*i;
        temp= temp*(t+i);
    }
    return sum;
}
 
double F3(int i, int n) {
    return pow(f(v[n])-f(v[n-1]), i);
}
 
double f(double x) {
    return pow(2,x);
}
 
void Print(int N)
{
    cout << "x[i]           ";
    cout << "f(x[i])            ";
    cout << "P(x[i])            \n\n";
    for(int i=0;i<N;i++)
    {
        cout << "x[" << i << "] = " << v[i] << "     ";
        cout << "f(x[" << i << "]) = " << f(v[i]);
        cout << "     P(x[" << i << "]) = " << Polinom(N, v[i]) << "\n";         
    }
    cout << "\n";
}
 
void clean(int N)
{
    free(v);
    v=NULL;
}