function res = jprod( k,t,p )

s = 0;

for j = 1:(k+1)
    for jj = 1:(k+1)
        s = s +momexp(t,j,p)*conj(momexp(t,jj,p));
    end;
end;

res = s;
end

