Multi-Gaussian closure method for randomly excited non-linear systems