  We describe a procedure to systematically improve direct diagonalization
results for few-particle systems trapped in one-dimensional harmonic potentials
interacting by contact interactions. We start from the two-body problem to
define a renormalization method for the interparticle interactions. The
procedure is benchmarked with state-of-the-art numerical results for three and
four symmetric fermions.
