Abstract: A spline fit method and program for multi-set data have been developed.Improvements have been made to have new functions: any order of spline as base.Knot optimization and accurate calculation for error of fit value. The programhas been used for practical evaluation of nuclear data.