Recursive legendre interpolatory product quadrature is working.
Playing around with setting up a PyPI package.