Electrostatic interaction energies from a generalized Gaussian quadrature