По долгу службы передо мной стоит задача взять третью производную от множество точек {х,y} n ~ 2000. Ничтоже сумняшися я решил для каждой точки построить интерполяционный многочлен Лагранжа (L(x), четвертой степени)), а потом его продифференцировать аналитически.
Ни смотря на то, что алгоритм получения L(x), как и дифференцирование многочлена работают верно по-отдельности, вместе они не дают правильного результата.
Существуют ли какие-либо извесные проблемы в этой области?