The integral bottleneck in evaluating molecular energies arises from the two-electron contributions. These are difficult and time-consuming to evaluate, especially over exponential type orbitals, used here to ensure the correct behaviour of atomic orbitals. In this work, it is shown that the two-centre Coulomb integrals involved can be expressed as one-electron kinetic-energy-like integrals. This is accomplished using the fact that the Coulomb operator is a Green's function of the Laplacian. The ensuing integrals may be further simplified by defining Coulomb forms for the one-electron potential satisfying Poisson's equation therein. A sum of overlap integrals with the atomic orbital energy eigenvalue as a factor is then obtained to give the Coulomb energy. The remaining questions of translating orbitals involved in three and four centre integrals and the evaluation of exchange energy are also briefly discussed. The summation coefficients in Coulomb forms are evaluated using the LU decomposition. This algorithm is highly parallel. The Poisson method may be used to calculate Coulomb energy integrals efficiently. For a single processor, gains of CPU time for a given chemical accuracy exceed a factor of 40. This method lends itself to evaluation on a parallel computer. †Dedicated to Professor Nicholas Handy.