function r = legroots(N);%  The function r = legroots(N) computes the roots of the %  Legendre polynomial of degree N.%  J.A.C. Weideman, S.C. Reddy 1998.n = [1:N-1];                   %  Indicesd = n./sqrt(4*n.^2-1);         %  Create subdiagonalsJ = diag(d,1)+diag(d,-1);      %  Create Jacobi matrixr = sort(eig(sparse(J)));      %  Compute eigenvalues