close all;
h = 340;
% for a = [sqrt(0.5) 1/4 0.50280655082539]
for a = [sqrt(0.5) 1/4]
h = h+1;
h=figure(h); 
hold on;
b = a;
L = 5;
t = [-L:1/200:-eps;eps:1/200:L];
j = a*cosh(t)-i*b*sinh(t);
c = -i*0.25; j2 = a*cosh(t+c)-i*b*sinh(t+c); j3 = a*cosh(t-c)-i*b*sinh(t-c);
plot(j','k-','LineWidth',2);
if h==341
    plot(j2','r:','LineWidth',2);
    plot(j3','r:','LineWidth',2);
end
axis equal;
axis([-5 5 -4 4]);ax = axis;
set(gca,'FontSize',18,'FontWeight','normal','FontName','Times',...
    'XTick',[-5:5],'YTick',[-4:4],'Box','on');
k = [(1-2*j+sqrt(8*j-7))/2;(1-2*j-sqrt(8*j-7))/2];
plot(k','b-.','LineWidth',2);
k = [(3-2*j+sqrt(1-8*j))/2;(3-2*j-sqrt(1-8*j))/2];
plot(k','b-.','LineWidth',2);
plot([0 0],ax(3:4),'k-')
plot(ax(1:2),[0,0],'k-')
plot([1:5],0*[1:5],'.','MarkerSize',18);
xlabel('Re({\itz})','FontWeight','normal','FontSize',24);
ylabel('Im({\itz})','FontWeight','normal','FontSize',24);

% print(gcf,'-deps2c',strcat(strcat('loci',num2str(h-340)),'.eps'))

end