% Generates Figure 8.2 (left)

q = linspace(0,1,1000);
th = 0*q;
for j=1:length(q)
    th(j)=theta(q(j));
end
plot(q,th,'k-','LineWidth',2);
axis([-0.008 1.008 -0.02 pi/2+0.02])
set(gca,'XTickMode','manual','XTick',[0 1],...
        'YTickMode','manual','YTick',[0 pi/2],'YTickLabel',{'0';''});
set(gca,'FontSize',18,'FontWeight','normal','FontName','Times');
grid on;
xlabel('{\itq}','FontWeight','normal','FontSize',24);
ylabel('\theta({\itq})','FontWeight','normal','FontSize',24);
text(-0.015,pi/2,'\pi/2','FontWeight','normal','FontSize',18,...
         'HorizontalAlignment','right','FontName','Times')
% print -deps2c theta1.eps