global mc mp iq idelta irout AB c alpha beta eps0
alpha=-1/2; beta=alpha;
eps0=.5e-12;
mc=1; mp=0; iq=0; idelta=1;
irout=1; AB=[-1 1]; Mmax=3000;
for omega=pi*[1/2 9/10 99/100]
  c=sin(omega/2);
  for n=[5:5:20 30:10:100 200 300]
    [ab,Mcap,kount]=mcdis(n,eps0,@quadgp,Mmax);
    ab(1,2)=2*c*ab(1,2);
    abv=trigauss(n,omega,pi/2,pi);
    maxerr=max(abs(ab(:,2)-abv(:,2)));
    fprintf('c=%14.12f   n=%3.0f   maxerr=%11.4e\n',c,n,maxerr)
  end
  fprintf('\n')
end
