function r = herroots(N);
% The function r = herroots(N) computes the roots of the
% Hermite polynomial of degree N.
% J.A.C. Weideman, S.C. Reddy 1998.
J = diag(sqrt([1:N-1]),1)+diag(sqrt([1:N-1]),-1); % Jacobi matrix
r = sort(eig(sparse(J)))/sqrt(2); % Compute eigenvalues
