Как построить рекурсивную функцию сплайна в C ++

Сейчас я работаю над методом решения дифференциальных уравнений, который называется базис-сплайн-коллокацией. У меня возникли проблемы с созданием метода построения сплайна произвольного порядка с соотношением с начальным условием B 1 i (x)={ 1Вк + 1я( х ) = х - хяИкск + я- хяВКя+ хk + i + 1- хИксk...