% FIGURE2_2 Computations to generate Figure 2.2.
%
clear mc mp iq idelta irout
global mc mp iq idelta irout
global c0
N=5; Mmax=6; 
mc=2; mp=0; iq=1; idelta=2; irout=1; eps0=100*eps;
hold on
abl=r_jacobi(5);
b1=zeros(1,131); b2=zeros(1,131); b3=zeros(1,131);
b4=zeros(1,131);
for i=1:131
  if i<=51
    c0=(i-1)/5;
    c(i)=c0;
  elseif i<=71
    c0=10+(i-51)/2;
    c(i)=c0;
  else
    c0=i-51;
    c(i)=c0;
  end
  [abcl,Mcap,kount]=mcdis(N,eps0,@quadcl,Mmax);
  b1(i)=abcl(2,2);
  b2(i)=abcl(3,2);
  b3(i)=abcl(4,2);
  b4(i)=abcl(5,2);
end
plot(c,b1)
plot(c,b2)
plot(c,b3)
plot(c,b4)
axis([0 95 .2 .55])
xlabel('c')
text(83,.33,'\beta_1','FontSize',12)
text(83,.265,'\beta_2','Fontsize',12)
text(87,.255,'\beta_3','FontSize',12)
text(83,.24,'\beta_4','FontSize',12)
hold off
