% QUADLAG  A special quadrature routine used in Table2_18.
%
function xw=quadlag(N,i)
ab=r_laguerre(N);
xw=gauss(N,ab);
xw(:,2)=xw(:,2)./((1+exp(-xw(:,1))).^2);
if i==1, xw(:,1)=-xw(:,1); end
