spline_interpolation 样条插值代码

/***************************************************************************
 * Func: spline_interpolation                                              *
 *                                                                         *
 * Desc: takes 16 unsigned characters in an array and outputs the results  *
 *       of a B-spline interpolation.  The sample points of the            *
 *       interpolation function f(-1-d) and f(-d) are evaluated as f(1+d)  *
 *       and f(d).  The computed constants are a0 = f(1+d)                 *
 *                                             a1 = f(d)                   *
 *                                             a2 = f(1-d)                 *
 *                                             a3 = f(2-d)                 *
 *                                                                         *