Natural interpolate splines of arbitrary order

The natural critérium was proposed for extra coefficients of interpolative splines of arbitrary power: this is minimization of the highest derivative breaks. The simple algorithm was constructed for calculation of these splines coefficients. Errors of splines of different powers were investigated as theoretically so on numerical tests. The weight function was found which improved accuracy essentially near boundaries. Some recommendations were given for applied calculations.