% ex4

L1=[];
L2=[];
L3=[];
for p = 0:0.001:0.2
   L=rloc_ex4(p);
   L1=[L1;L(:,1)];
   L2=[L2;L(:,2)];
   L3=[L3;L(:,3)];  
end

plot([0,0],[-10,10],'k',real(L1),imag(L1),'.',real(L2),imag(L2),'x',real(L3),imag(L3),'.');
hold on
L=rloc_ex4(0.06361); % val. caract. duplo
plot(real(L(:,1)),imag(L(:,1)),'ko');
L=rloc_ex4(0.0889); % ponto limite
plot(real(L(:,1)),imag(L(:,1)),'mo',real(L(:,2)),imag(L(:,2)),'mo');
L=rloc_ex4(0.10574); % ponto limite
L0=rloc_ex4(0.1057555);
plot(real(L0(:,1)),imag(L0(:,1)),'go',real(L(:,2)),imag(L(:,2)),'go');
L=rloc_ex4(0.1309); % hopf
plot(real(L(:,1)),imag(L(:,1)),'ro');

xlabel('Re(\lambda)'); ylabel('Im(\lambda)');