Метод интерполяции Эрмита и полинома Лагранжа

 double f(double X)
 {
     return Math.Pow(X, 3) + 3 * Math.Pow(X, 2) - 24 * X +10;
 }
 
 
 //Для Эрмита
 double df(double x)
 {
     return 3 * (x * x) + 6 * x - 24;
 }
 
 //Полином Лагранжа
 double Lag(int i, double[]X, double x)
 {
     int n = X.Length;
     double result = 0;
     for (int k = 0; k < n; k++)
     {
         if (k == i) continue;
         double term = 1;
         for (int j = 0; j < n; j++)
         {
             if (j == i || j == k) continue;
             
                 
                term *= (x - X[j]) / (X[i] - X[j]);
             
         }
         result += term/(X[i] - X[k]);
     }
 
     return result;
 }
 
 
 double Ermith( double[] xar, double[] yar, double[] dyar, double x)
 {
     int n = xar.Length;
     double result = 0;
 
     for (int i = 0; i < n; i++)
     {
         double Li = 1;
         double dLi = 0;
 
         for (int j = 0; j < n; j++)
         {
             if (j != i)
                 Li *= (x - xar[j]) / (xar[i] - xar[j]);
         }
 
         //for (int j = 0; j < n; j++)
         //{
         //    if (j != i)
         //    {
         //        double term = 1 / (xar[i] - xar[j]);
 
         //        for (int k = 0; k < n; k++)
         //        {
         //            if (k != i && k != j)
         //                term *= (x - xar[k]) / (xar[i] - xar[k]);
         //        }
 
         //        dLi += term;
         //    }
         //}
 
         dLi = Lag(i, xar, x);
 
 
         double Hi = (1 - 2 * (x - xar[i]) * dLi) * Li * Li;
         double Ki = (x - xar[i]) * Li * Li;
 
         result += Hi * yar[i] + Ki * dyar[i];
     }
 
     return result;
 }
 
 private void buttonX_Click(object sender, EventArgs e)
 {
         double x = double.Parse(Xpoint.Text, CultureInfo.InvariantCulture); //Current Culture чтобы понимал 1,5 вместо 1.5  //теперь не работает после обновления. Непонятно почему
         double[] Xnodes = Xn.Text.Split(',').Select(v=>double.Parse(v.Trim(),CultureInfo.InvariantCulture)).ToArray();
         double[] Ynodes = Yn.Text.Split(',').Select(v => double.Parse(v.Trim(), CultureInfo.InvariantCulture)).ToArray();
        
        //....
 
 
         ....//
 
//Эрмит
 
     double[] Dy = Xnodes.Select(xi => df(xi)).ToArray();
     double E = Ermith(Xnodes, Ynodes, Dy, x);
     ResultX.Text = $"Сплайн= {Sx}\r\n" + $"Функция = {Fx}\r\n" + $"Погрешность= {Math.Abs(Sx - Fx):F6}\r\n" + $"Эрмит ={E:F6}\r\n" + $"Погрешность Эрмита = {Math.Abs(E - Fx):F6}\r\n";
 }