“…In terms of the earlydeveloped cubic spline (Wang and Kahawita [12]), if uniform grids are taken, then the first derivative and second derivative have fourth-order and second-order accuracy, respectively, hence numerical solution of the spline has been widely applied. There is much recent research based on the spline theory, such as the high-precision C 1 -cubic spline collocation method [13], the quartic spline method [14], the residual error-correction method [15], and the parameter spline method [16,17].…”