Exploiting the fact that the Fock matrix is symmetric when building it.
When building Hermite Coeffs E: iterating up to i=iMax and i=jMax and not i=m_angMom and j=m_angMom.
Optimizing the calculation of pow(-Dx,k) in function tabulated() of class Boys.
Fixed small bug: in function tabulated() of class Boys: initialize Dx_pow to 1.