%SMMOMENU 
%
function mom=smmomEnu(dig,N,nu,a)
disp('running smmomEnu')
digits(dig);
mom=vpa(zeros(2*N,1));
abj=sr_jacobi01(dig,N,0,nu-1);
abl=sr_laguerre(dig,N);
abla=sr_laguerre(dig,2*N,a);
xwj=sgauss(dig,N,abj);
xwl=sgauss(dig,N,abl);
for k=1:2*N
  sig=vpa(0);
  for r=1:N
    for s=1:N
      x=xwj(s,1)*xwl(r,1);
      w=xwj(s,2)*xwl(r,2);
      if k==1
        p=vpa(1);
      else
        p0=vpa(0); p=vpa(1);
        for j=1:k-1
          pm1=p0; p0=p;
          p=(x-abla(j,1))*p0-abla(j,2)*pm1;
        end
      end
      sig=sig+w*p;
    end
  end
  mom(k)=sig;
end
