clear all
mu1 = [1:99]*pi/200;
mu2 = pi/2+mu1;
top = 5;
bottom = -5;
figure(1)
plot([mu1,mu2],[mu1,mu2]/0.5, 'linewidth',2)
hold on
plot([mu1,mu2],[mu1,mu2] ,[mu1,mu2],[mu1,mu2]/2.5,'linewidth',2)
 
plot(mu1,-tan(mu1),mu2,-tan(mu2),[0,pi],[0,0],'k--','linewidth',2)
 
plot([pi/2,pi/2],[top,bottom],'--')
axis([0 pi bottom top])
xlabel('\mu','fontsize',20)
text(1,3,'a=0.5','fontsize',18)
text(2.5,3,'a=1','fontsize',18)
text(2.5,1.5,'a=2.5','fontsize',18)

hold off

