clear all;
close all;

gamma = 0.5;
omega = 5.0;
Dt = 0.01;
l1 = -gamma+sqrt(gamma^2-omega^2);
l2 = -gamma-sqrt(gamma^2-omega^2);

ar = linspace(-6,6,300);
ai = linspace(-6,6,300);

figure(1)
handle = axes
set(handle,'FontSize',16)
Z = stab_theta(0.0,ar,ai);
[cs,handle] = contour(ar,ai,Z,[1 1],'-');
set(handle,'LineWidth',3);
axis('square'), axis([-0.6 0.6 -0.6 0.6]), grid on
hold on

Z = stab_theta(0.25,ar,ai);
[cs,handle] = contour(ar,ai,Z,[1 1],'-.');
set(handle,'LineWidth',3);

handle = plot(real(l1*0.1),imag(l1*0.1),'o')
set(handle,'MarkerSize',12);
handle = plot(real(l2*0.1),imag(l2*0.1),'s')
set(handle,'MarkerSize',12);

handle = plot(real(l1*0.01),imag(l1*0.01),'x')
set(handle,'MarkerSize',12);
handle = plot(real(l2*0.01),imag(l2*0.01),'+')
set(handle,'MarkerSize',12);

legend('\theta=0.0','\theta=0.25','0.1\lambda_1','0.1\lambda_2','0.01\lambda_1','0.01\lambda_2',3)

hold off

figure(2)
handle = axes
set(handle,'FontSize',16)
Z = stab_theta(0.0,ar,ai);
[cs,handle] = contour(ar,ai,Z,[1 1],'-');
set(handle,'LineWidth',3);
axis('square'), axis([-0.1 0.1 -0.1 0.1]), grid on
hold on

Z = stab_theta(0.25,ar,ai);
[cs,handle] = contour(ar,ai,Z,[1 1],'-.');
set(handle,'LineWidth',3);

handle = plot(real(l1*0.01),imag(l1*0.01),'x')
set(handle,'MarkerSize',12);
handle = plot(real(l2*0.01),imag(l2*0.01),'+')
set(handle,'MarkerSize',12);

legend('\theta=0.0','\theta=0.25','0.01\lambda_1','0.01\lambda_2',3)

hold off