Response power spectrum of multi-degree-of-freedom nonlinear systems by a Galerkin technique