A Galerkin method for the steady state analysis of harmonically excited non-linear systems