def fit(lcs, knotstep=20.0, n=None, knots=None, stab=True,
	stabext=300.0, stabgap=20.0, stabstep=5.0, stabmagerr=-2.0, stabrampsize=0, stabrampfact=1.0,
	bokit=1, bokeps=2.0, boktests=5, bokwindow=None, k=3, verbose=True):
	"""
	The highlevel function to make a spline fit.