 
d=30;
%--------------------------------FIRST AREA
r=110;
cx=d;
 plot( -cx,0,'+');
 hold on
 plot( cx,0,'+');
 cy=0;
 t=0.5*pi:pi/100:3*0.5*pi;
 plot( cx+cos(t) * r,cy + sin(t) * r);
 t=-0.5*pi:pi/100:0.5*pi;
 cx=-d;
 plot( cx+cos(t) * r,cy + sin(t) * r);
 
 %----------------------------AXIS
 x=-100:0.1:100;
 y=-80:0.1:80;
 plot( x,0,'k');
 hold on
plot( 0,y,'k');
 
 
 %------------------------------------------2nd area
     beta=deg2rad(50);
     r=d*csc(beta);
     cx=0;
     cy=d*cot(beta);
     plot(0,cy)
     t=0:pi/100:2*pi;
     hold on
     plot( cx+cos(t) * r,cy + sin(t) * r,'r');
     cy=-d*cot(beta);
     plot( cx+cos(t) * r,cy + sin(t) * r,'r');
     plot( 0,-cy,'+r');
     plot(0,cy,'+r')
