% test
set(0,                           ...
   'defaultaxesfontsize', 20,   ...
   'defaultaxeslinewidth', 1.0, ...
   'defaultlinelinewidth', 1.2, ...
   'defaultpatchlinewidth', 0.7);
thp=[0.01:0.01:0.99]*pi;

f = (sin(thp)./thp-cos(thp))./sin(thp);
 
g =thp-atan(f);
 
figure(1)
plot(thp/pi,pi./g, 'linewidth',2) 
xlabel('\theta_p/\pi','fontsize',20)
ylabel('\alpha R \pi/V','fontsize',20)
axis([0 1 0 10])
