mpspline_est1 function

Estimate spline parameters