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