Concise, open‐ended implementation of Rys polynomial evaluation of two‐Electron integrals

A quadrature‐point‐driven implementation of the standard Rys polynomial method for computing two‐electron repulsion integrals of gaussian basis functions has been found to be both concise and openended with respect to the angular momentum of the gaussian functions (i.e., s,p,d,f,g,…). These are important features in certain applications, such as molecular properties and property gradients.