Нужно сделать программу для интерполяции по формуле Лагранжа, y=Sin(x^3), границы от 0-2 с шагом 0,08. Почему не работае

#include <math.h>
#include <stdio.h>
#include <process.h>
#include <iostream>
#include <iomanip>
#include <cmath>
#include <fstream>
 
using namespace std;
 
double lagrange(double *x1, double *y, short n, double _x)
{
   double result = 0;
   for (short i = 0; i < n; i++)
  {
     double P = 1;
     for (short j = 0; j < n; j++)
    {
      if (j != i)
        P *= (_x - x1[j]) / (x1[i] - x1[j]);
    }
    result += P * y[i];
  }
  return result;
}
 
int main()
{
  setlocale(LC_ALL, "Russian");
 
  int t = 0, a = 0, b = 2;
  double h = 0.08;
  double* x= new double[25];
  double* F = new double[25];
  double* dF = new double[25];
  double lagrange(double *x1, double* y, short n, double _x);
 
  cout << "Значения х" << "\t" << "Экспонента" << endl;
  while (t<25)  // подправлю сразу и первый пункт
  {
    x[t] = a + h * t;
    F[t] = sin(pow(x[t], 3));
    t++;
  }
// вот здесь разделила. Сначала посчитаем ВСЕ x[t] и F[t]
  t=0;
  while (t<25)
 {
    dF[t] = fabs(F[t] - lagrange(x, F, 25, t));
    cout << fixed << showpoint << setprecision(6);
    cout << x[t] << "\t" << F[t] << "\t" << lagrange(x, F, 25, t) << "\t" << endl;
    t++;
  }
 
  delete[] x;
  delete[] F;
  delete[] dF;
}

#include <math.h>
#include <stdio.h>
#include <process.h>
#include <iostream>
#include <iomanip>
#include <cmath>
#include <fstream>
 
using namespace std;
 
double lagrange(double *x1, double *y, const short n, const double _x)
{
   double result = 0;
   for (short i = 0; i < n; i++)
  {
     double P = 1;
     for (short j = 0; j < n; j++)
    {
      if (j != i)
        P *= (_x - x1[j]) / (x1[i] - x1[j]);
    }
    result += P * y[i];
  }
  return result;
}
 
int main()
{
  setlocale(LC_ALL, "Russian");
 
  const short pointsNumber=25; // количество точек
  const short polynomialDegree=3; // степень полинома
 
  int t = 0;
  const int a = 0;  // вы же не собираетесь менять значение этой переменной
                    // переменную b = 2; вы не используете
  const double h = 0.08; // и значение этой не меняете, так что const
  double* x= new double[pointsNumber];
  double* F = new double[pointsNumber];
  double* dF = new double[pointsNumber];
  double lagrange(double *x1, double* y, short n, double _x);
 
 
// первый раз проходим и считаем значения x[t] и F[t] в каждой точке
  cout << "Значения х" << "\t" << "Экспонента" << endl;
  while (t<pointsNumber)  // подправлю сразу и первый пункт
  {
    x[t] = a + h * t;
    F[t] = sin(pow(x[t], 3));
    t++;
  }
// второй раз идем и считаем лагранжа
  t=0;
  while (t<pointsNumber)
 {
    dF[t] = fabs(F[t] - lagrange(x, F, polynomialDegree, x[t]));
    cout << fixed << showpoint << setprecision(6);
    cout << x[t] << "\t" << F[t] << "\t" << lagrange(x, F, polynomialDegree, x[t]) << "\t" <<dF[t]<<"\t" <<endl;
    t++;
  }
 
  delete[] x;
  delete[] F;
  delete[] dF;
}