Programmatically generate Gauss-Legendre quadrature table

[GregVernon]

Lets see if we can programmatically (using symbolic toolbox) generate the Gauss-Legendre quadrature table (see below)

Things we'll need:

• [x] Legendre basis polynomials
• [x] Solve roots symbolically via solve(_,'MaxDegree',p) (documentation)
• [ ] Solve roots numerically via vpasolve (documentation)

[GregVernon]

Good work so far @hopdahl

Note that to complete the GLQ table we need to compute the x_i and w_i values.

x_i is simply the roots of the Legendre polynomials w_i is then: