Comment construire une fonction spline récursive en C ++

En ce moment, je travaille sur une méthode de résolution d'équations différentielles appelée collocation base-spline. Ce avec quoi j'ai du mal à construire une méthode pour construire une spline d'ordre arbitraire, avec la relation avec la condition initiale B 1 i (x)={ 1Bk + 1je( x ) = x - xjeXk +...