Use a better orthogonalization procedure for basis finding in quadrature finder

The procedure in [this paper](https://arxiv.org/pdf/2307.14494) could replace this code:

https://github.com/inducer/modepy/blob/8f06f1c6e3573ea03dd22f7006d98ccc30ec7873/modepy/quadrature/construction.py#L171-L188

(not exactly, since the initial basis is fixed, but morally)

Read the procedure in modepy/quadrature/construction.py at lines 171-188 alongside the linked paper. Determine how the paper's orthogonalization procedure applies when the initial basis is fixed. Done means the existing basis-finding procedure is replaced consistently and the quadrature construction still works.

Written by the indexing model from the issue text.