Return weighted sum of squared residuals of the spline approximation: sum((w[i] * (y[i]-spl(x[i])))**2, axis=0).
scipy.interpolate.LSQUnivariateSpline.get_knots
scipy.interpolate.LSQUnivariateSpline.integral