% QUADFERMI A quadrature routine used in Example 2.41.
%
function xw=quadfermi(N,i)
global r
ab=r_laguerre(N); xw=gauss(N,ab);
xw(:,1)=xw(:,1)/r;
xw(:,2)=xw(:,2)./(r*(1+exp(-xw(:,1))).^r);
