+ vector<double> smoothSpline(vector<double>, vector<double>, int);
+ vector<double> calc_qvalues(vector<double>&);
+ vector<double> sknot1(vector<double>);
+ int nkn(int);
+ int OrderPValues(int, int, vector<double>&, vector<int>&);
+ int swapElements(int, int, vector<double>&, vector<int>&);
+ vector<int> getSequence(int, int, int);
+
+ int spline(vector<double>&, vector<double>&, int, int, vector<double>&);
+ int splint(vector<double>&, vector<double>&, double&, double&, vector<double>&);
+
+