Wrong polynomial order for Gauss Lobatto Quadratures

[inversecrime]

Hi, thanks for creating this package. It is written that the order of the Gauss Lobatto Quadrature is 2N-3 for N+1 points, which is not true. It's 2N-1. For example, with 2 points, it is of order 1 (and not -1).
This needs to be changed in the documentation and also in the exact_to property.
Thanks and best wishes!

(Obviously this depends on what you define as "order". I use it synonymously to "exact up to degree")

[inducer]

Thanks for the report! You're right, obviously. It doesn't make sense that Gauss-Lobatto would lose four exactness degrees compared to Gauss. Hard to believe this was overlooked all these years, but here we are. 🤦

#163 should fix this.