Analytic evaluation of Coulomb integrals for one, two and three-electron distance operators, $R_{C1}^{-n}R_{D1}^{-m}$, $R_{C1}^{-n}r_{12}^{-m}$ and $r_{12}^{-n}r_{13}^{-m}$ with $n, m=0,1,2$