% QUADCL A quadrature routine used in Figure2_2.m.
%
function xw=quadcl(N,i)
global c0
if i==1
  abc=r_jacobi(N,-1/2);
  xw=gauss(N,abc);
else
  abl=r_jacobi(N);
  xw=gauss(N,abl);
  xw(:,2)=c0*xw(:,2);
end
