close all
a=2;b=1.2;
c=sqrt(a^2+b^2);
C='black';
y=-3:.1:3;
xa=sqrt(1+y.^2/b^2)*a;
xd=-sqrt(1+y.^2/b^2)*a;
figure;
plot(xa,y,'linewidth',2,'color',C);
hold on;
plot(xd,y,'linewidth',2,'color',C);
ylim([-4,4]);
plot([0,0],[-4,4],'linewidth',2,'color',C);
fill([-0.15,0,0.15],[3.5,4,3.5],'k');
plot([-6,6],[0,0],'linewidth',2,'color',C);
fill([5.5,6,5.5],[0.12,0,-0.12],'k')
xf=a^2/c;
plot([xf,xf],[-3.2,3.2],'linewidth',2,...
'linestyle',':','color',C);
plot(c,0,'.','color',C,'markersize',16);
plot(-c,0,'.','color',C,'markersize',16);
text(5.6,-0.5,'X',&#