NB: 0th, 1st, and 2nd moments generated using: + 2,000 pt. Gauss-Legendre grid/interval + relative error tolerance 10**-8 + 0.02 *omega_p(0) intervals + Im omega = 10**-10 * omega_p(0) 3rd moment sum rule data generated with: + relative error tolerance 10**-7 for both the 3rd frequency moment, and the 2D integral convergence + min_cutoff=2 + cubic spline interpolation of spacing 0.005 k_F + same parameters driving the sum rule calculation, besides relative error tolerance The parameters for the eps_c integration were unchanged from the routines (i.e., the defaults set in get_ec.py)