zplagu.m

various matlab programs to slove various problems

function [x,w]=zplagu(n)
if (n <= 1), disp(' n must be > 1 '); return; end
[a,b]=coeflagu(n); 
JacM=diag(a)+diag(sqrt(b(2:n)),1)+diag(sqrt(b(2:n)),-1);
[w,x]=eig(JacM); x=diag(x); w=w(1,:)'.^2;