zplagu.m

来自「various matlab programs to slove various」· M 代码 · 共 6 行

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;