Interpolacion de hermite

// Funcion para calcular h(t)
float h(float t, int i, int j)
{
  float di, dj;
  if (i == 0)
    di = 0;
  else
    di = 1;
  if (j == 0)
    dj = 0;
  else
    dj = 1;
  if (i == 0 && j == 0)
    return 2 * pow(t, 3) - 3 * pow(t, 2) + 1;
  else
    return pow(t, 3) - 2 * pow(t, 2) + t;
}