https://www.bragitoff.com/2018/02/determining-roots-legendre-polynomials-weights-gaussian-quadrature-c-program/gauss-legendre-quadrature-output-c-program-n12/