Coulomb matrix elements in 2D harmonic oscillator basis