add Gauss-Laguerre Quadrature to HF code