Variational methods in a Harmonic Oscillator basis