Интерполяционный полином Ньютона для интерполяции назад

#include <iostream>
#include <conio.h>
using namespace std;
 
double T(double t, int n)
{
    double temp = t;
    for (int i = 1; i < n; i++)
        temp = temp * (t + i);
    return temp;
}
 
int factorial(int n)
{
    int f = 1;
    for (int i = 2; i <= n; i++)
        f *= i;
    return f;
}
 
int main() 
{
    setlocale(LC_ALL, "Russian");
 
    const int n = 6;
 
    double x[] = { 0, 0.31, 0.62, 0.93, 1.24, 1.55 };
    double y[n][n];
    y[0][0] = 2;
    y[1][0] = 1.96;
    y[2][0] = 1.793;
    y[3][0] = 1.408;
    y[4][0] = 0.909;
    y[5][0] = 0.479;
 
    double X;
    cout << "Значение точки х: ";
    cin >> X;
 
    for (int i = 1; i < n; i++) {
        for (int j = n - 1; j >= i; j--)
            y[j][i] = y[j][i - 1] - y[j - 1][i - 1];
    }
 
    double sum = y[n - 1][0];
    double t = (X - x[n - 1]) / (x[1] - x[0]);
    for (int i = 1; i < n; i++) {
        sum = sum + (T(t, i) * y[n - 1][i]) /factorial(i);
    }
 
    cout << "Приближённое значение функции ln(x) = " << sum << endl;
 
    _getch();
    return 0;
}