Nuclear attraction integrals are implemented (one nucleus at a time).
Clarified the source of chebyshev-based Boys
Chebyshev-based Boys can be initialized with negative mmax, this will avoid initialization of the table, etc. Aimed at generic code.
Simple API now can handle nuclear attraction integrals.
TwoBodyEngine works for Coulomb kernel