Numerical Integration in Density Functional Methods with Linear Combination of Atomic Orbitals