Coulomb interaction matrix elements in 2D harmonic oscillator basis