% 
set(0,                           ...
   'defaultaxesfontsize', 20,   ...
   'defaultaxeslinewidth', 1.0, ...
   'defaultlinelinewidth', 1.2, ...
   'defaultpatchlinewidth', 0.7);
m = [0:.01:1];

a = -m*pi.*cos(m*pi)./sin(m*pi);

figure(1)
plot(a,m,'linewidth',2)
ylabel('\mu ','fontsize',20)
xlabel('a','fontsize',20)
axis([0 20 0 1])