Полином Ньютона не выдаёт нужного результата

namespace Newton
{
    class Program
    {
        private static double h;//шаг
        private static double[] X = { 0.0, 1.0, 2.0, 3.0, 4.0};//узлы
        private static double[] Y = { 1.0, -1.256, -9.361, -22.863, -57.551};//значения функции в узлах
        public static int n = X.Length;//кол-во узлов
        // Function y = sin(x^3) + 3cos(x^2) - e^x - 1
        static double Function(double x)
        {
             double result = Math.Sin(Math.Pow(x, 3)) + 3 * Math.Cos(Math.Pow(x, 2)) - Math.Exp(x) - 1;
            return result;
        }
        //Полином Ньютона 
        static double Polinom(double x)
        {
            double result = Cj(0);
            for (int j = 1; j < n; j++)
            {
                double something = Cj(j);
                for (int k = 0; k < j; k++) 
                        something *= x - X[k];
                result += something;
            }
            return result;
        }
 
        static double Cj(int j)
        {
                if (j == 0)
                    return Y[0];
                    double delta = Delta(j, j);
            return delta / (Factorial(j) * Math.Pow(h, j));
        }
 
        static double Delta(int p, int i)
        {
            if (p == 1) 
            return Y[i] - Y[i - 1];
            return Delta(p - 1, i) - Delta(p - 1, i - 1);
        }
 
        static double Factorial(double value)
        {
            if (value == 0)
                 return  1;
            return value * Factorial(value - 1);        
        }
        static void Main(string[] args)
        {Link:
            h = X[1] - X[0];
            Console.Write("Точка = ");
            double point = double.Parse(Console.ReadLine());
            Console.WriteLine("Исходная функция в точке " + point + "  " + Function(point));
            Console.WriteLine("Полином Ньютона в точке " + point + "   " + Polinom(point));
            Console.ReadKey();
            goto Link;
        }
    }
}

using System;
using System.Collections.Generic;
 
class NewtonApproximator {
    
    private List<(double b, double x)> records = new List<(double, double)>();
    
    public void AddPoint(double x, double y) {
        var t = DoInterpolate(x);
        double b = (y - t.result) / t.product;
        records.Add((b, x));
    }
    
    public double Approximate(double x) {
        return DoInterpolate(x).result;
    }
    
    private (double result, double product) DoInterpolate(double x) {
        double result = 0.0;
        double product = 1.0;
        foreach (var r in records) {
            result += r.b * product;
            product *= x - r.x;
        }
        return (result, product);
    }
}
 
class Program {
    
    public static void Main(string[] args) {
        //Данные для аппроксимации
        Func<double, double> function = (double x) => Math.Sin(Math.Pow(x, 3)) + 3 * Math.Cos(Math.Pow(x, 2)) - Math.Exp(x) - 1;
        double[] valuesOfX = new[]{0.0, 1.0, 2.0, 3.0, 4.0};
        
        //Аппроксимация
        var approximator = new NewtonApproximator();
        foreach (var x in valuesOfX){
            approximator.AddPoint(x, function.Invoke(x));
        }
        
        //Вывод результата
        for(double x = -1.0; x <= 5.0; x += 0.5) {
            Console.WriteLine("x = " + x + ", approximated = " + approximator.Approximate(x) + ", real = " + function.Invoke(x));
        }
    }
}