clear all;
close all;

xr = solve('abs((1+(1-th)*xr)/(1-th*xr))=1.0')

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


th = [0.0 0.25 0.5 0.53 0.75 1.0];

filename(1).s = '00.eps';
filename(2).s = '025.eps';
filename(3).s = '05.eps';
filename(4).s = '053.eps';
filename(5).s = '075.eps';
filename(6).s = '10.eps';

linestyle = char('r','g','b','m','c','y');
%linestyle = char('k','k','k','k','k','k');

i=1;
figure(i)
Z = stab_theta(th(i),ar,ai);
[cs,handle] = contourf(ar,ai,-Z,[-1 -1],linestyle(i));
%set(handle,'LineWidth',3);
axis('square'), axis([-6 6 -6 6]), grid on, colormap([0.8 0.8 0.8])
%title('theta')
hold on
%arrow([-5.5 0],[5.5 0],'Width',1), text(5.0, 0.5,'Re')
%arrow([0 -5.5],[0 5.5],'Width',1), text(0.5, 5.0,'Im')
arrow([-5.5 0],[5.5 0]), text(5.0, 0.5,'Re')
arrow([0 -5.5],[0 5.5]), text(0.5, 5.0,'Im')
hold off
print('-deps2',filename(i).s)
for i=2:6
  figure(i)
  Z = stab_theta(th(i),ar,ai);
  [cs,handle] = contourf(ar,ai,-Z,[-1 -1],linestyle(i));
  %set(handle,'LineWidth',3);
  axis('square'), axis([-6 6 -6 6]), grid on, colormap([0.8 0.8 0.8])
  %title('theta')
  hold on
  %arrow([-5.5 0],[5.5 0],'Width',1), text(5.0, 0.5,'Re')
  %arrow([0 -5.5],[0 5.5],'Width',1), text(0.5, 5.0,'Im')
  arrow([-5.5 0],[5.5 0]), text(5.0, 0.5,'Re')
  arrow([0 -5.5],[0 5.5]), text(0.5, 5.0,'Im')
  hold off
  print('-deps2',filename(i).s)
end