Matlab functions implementing non-parametric local linear smoothing (LLS) for one-dimensional curve fitting: yi = f(xi) + ei for i = 1, 2, ..., n. The toolbox contains: (i) the usual LLS estimator, see e.g. [1] (ii) a jump-preserving (two-sided) LLS