Molecular Integral for the Double Electron-Repulsion Potential: A Closed Form Underived for Forty Years

A closed form of the title integral formula over the Gaussian-type orbitals is derived for the first time. The present closed form is analytical as the multiple hyper-geometric function of five variables.