Spline interpolation library in cpp
So you have a series of points and, at each point, a supplied derivative?
Is a piecewise cubic sufficient or does it need higher derivative continuity?  If the former is ok, then Cubic Hermite Splines will do the job. (If you need, they can be trivially mapped into cubic Beziers)