Add Gauss-Lobatto quadrature rules up to 43rd-order

These rules, which include the "end points" of the region, are less efficient than the standard Gauss rules, but may be useful in certain applications, or for verification purposes.

See also issue #350.