%% 
clear;
clc;
close;
k=1;
% 
% 
% J=[4-k 2;-26 -8-17*k];
% [det(J) trace(J)]

ev=zeros(2,0)

k=-0:0.01:10;
 for kv=k
     J=[4-kv 2;-26 -8-17*kv];
     ev=[ev real(eig(J))];
     eig(J)
     real(eig(J));
     
end
plot(k,ev','LineWidth',2)
ylim([-8 2]);

