function maxeigen(n);
e=ones(n,1);
A=spdiags([-e 2.01*e -e],[-1 0 1],n,n);
D=diag(2.01*e,0);
I=eye(n);

fprintf('maximun eigenvalue =%.8f \n',eigs(I-D\A,1));