zpherm.m

various matlab programs to slove various problems

function [x,w]=zpherm(n) 
if (n <= 1), disp(' n must be > 1 '); return; end
[a,b]=coefherm(n); 
JacM=diag(a)+diag(sqrt(b(2:n)),1)+diag(sqrt(b(2:n)),-1);
[w,x]=eig(JacM); x=diag(x); scal=sqrt(pi); w=w(1,:)'.^2*scal; 
[x,ind]=sort(x); w=w(ind);