illcond.m

sareli<matlab科学计算>配套程序

function errore=illcond(alpha,n)

theta = pi/(n+1);
lambda = alpha - 2*cos([1:n]*theta);

lmax = max(lambda);
lmin = min(lambda);
A=alpha*eye(n)-diag(ones(n-1,1),-1)-diag(ones(n-1,1),1);
cond2=cond(A);
errore = abs(cond2-lmax/lmin)/lmax*lmin;

return