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