Вычислить интерполяционный многочлен Ньютона с равноотстоящими узлами

#include <stdio.h>
#include <iostream>
#include <conio.h>
#include <math.h>
#include <iomanip>
using namespace std;
double func(double x)
  {
return exp(x);
  }
int main()
  {
  setlocale(LC_ALL, "Russian");
int i, j, n, period = 1;
double lnx, pr, a, b;
double x[100], f[100][100];
cout << "Введите количество узлов n = " ;
cin >> n; 
while (n <= 0){
    cout << "Введите n больше 0" <<  endl;
    cout << "n = ";
    cin >> n;
  } 
  cout << "Введите a = ";
  cin >> a; 
  cout << "Введите b = ";
  cin >> b;
while (b < a)
  {
      cout << "Введите b, которое больше a " << endl;
      cout << "Введите b = ";
      cin >> b;
  } 
for (i = 0; i < n; i++)
  {
      cout << "Xi = ";
      cin >> x[i];
      f[i][0] = func(x[i]);
  }
 
for (j = 1; j < n; j++)
for (i = j; i < n; i++)
  f[i][j] = (f[i][j - 1] - f[i - 1][j - 1]);
 
double xtest;
 
int t;
 
  cout << "Количество точек = ";
  cin >> t;
const double h = (b - a) / n;
int iter;
for (iter = 0; iter < t; iter++)
  {
      xtest = h * iter + a;
 
  lnx = f[0][0];
  pr = 1;
for (i = 1; i < n; i++)
  {
  pr = pr * (iter - x[i + 1]);
  lnx = lnx + f[i][i] * pr;
  }
  cout << xtest;
  cout << " " << func(xtest);
  cout << " " << lnx;
  cout << " " << abs(func(xtest) - lnx) << endl;
 
}
 
  system("pause");
  return 0;
}