Gauss-Jacobi implementation
draft stroud quadrature