function fx = xlnx(x)
%
% f(x) = x .* log(x), assuming non-negative arguments
%
fx       = NaN(size(x));

mask     = find(x>0);
fx(mask) = x(mask).*log(x(mask));

mask     = find(x==0);
fx(mask) = 0;