Формулировка / код Эрмита Интерполяция через четыре точки

Поиск C / C ++ реализации Hermite Interpolation. Пожалуйста, предложите простую / читаемую / академическую программу.

Входные данные: массивы x и y двойников с 4 значениями каждый. Представляя 4 2d баллов. Обратите внимание, что у меня нет касательных со мной.

Вывод: теперь, учитывая любое новое значение x (в пределах диапазона заданных значений x), нужно найти новое y. Обратите внимание, что у меня нет параметра ‘t’ со мной.

Пожалуйста, предложите формулировку и покажите реализацию C / C ++.