  We have developed a McMurchie-Davidson-like recursion formula for efficient
evaluation of the Coulomb attraction and interaction matrix elements between
two-dimensional primitive Cartesian Gaussian type orbitals. We also present
recurrence schemes for combined position and differential operator integrals,
and three-center Gaussian integrals. The Cartesian Gaussian orbitals are
isotropic in the exponent, but with arbitrary centers and angular momentum.
