function mom=mmom_Enufin(N,nu,c)
%   MOM=MMOM_ENUFIN(N,NU,C) computes in Matlab double
%   precision the first 2N modified moments, relative
%   to monic shifted Jacobi01 polynomials with Jacobi
%   parameters alpha=0, beta=nu-1, transformed from the 
%   interval [0,1] to [0,C] of the subrange exponential 
%   integral weight function w(x)=E_nu(x) on [0,C].     
 
mom=zeros(2*N,1); can=zeros(N); canmin=zeros(2*N,1);
abl=r_laguerre(N); xwl=gauss(N,abl);
abj=r_jacobi01(2*N-1,0,nu-1); xwj=gauss(N,abj);
abjc(:,1)=c*abj(:,1); abjc(1,2)=(c^nu)*abj(1,2);
abjc(2:2*N-1,2)=c^2*abj(2:2*N-1,2); 
abstar=loadvpa('ab_powexp64_1_5.txt',9,64,100,2); 
abstar=double(abstar); xwstar=gauss(N,abstar);
tic
for k=1:2*N
  u1=0;
  for i=1:N
    u2=0;
    for j=1:N
      t1=xwj(j,2)*pk(xwl(i,1)*xwj(j,1),k-1,abjc);
      t2=xwstar(j,2)*pk(c+xwl(i,1)*xwstar(j,1),k-1,abjc);
      u2=u2+t1-t2;
    end
    u1=u1+xwl(i,2)*u2;
  end
  mom(k)=u1;
end 
toc
