clear all
N=20; % number of subdivisions

for n=1:N
    a(n)=1+(1.2-1)*n/N;
    A=[0 0 0 0 127 4 80; 0.6747 0.737 0 0 0 0 0; 0 0.0486 0.6610 0 0 0 0; 0 0 0.0147 0.6907*a(n) 0 0 0; 0 0 0 0.0518*a(n) 0 0 0; 0 0 0 0 0.8091*a(n) 0 0; 0 0 0 0 0 0.8091*a(n) 0.8089*a(n)]
    lambda(n)=eigs(A,1);
end

plot(a,lambda,a,ones(size(a)),'-')
xlabel('a')
ylabel('lambda')
