legendre-gauss Compute Legendre quadrature in C++ Uses the Legendre polynomial recursion relation and Newton's method to iteratively compute the Legendre Gauss nodes. Based on my old Matlab code below http://www.mathworks.com/matlabcentral/fileexchange/4540-legendre-gauss-quadrature-weights-and-nodes/content/lgwt.m