m = 1000;
x = 1:m;
y = mylogfactorial(m);

plot(x, y)

hold on

% educated guess
plot(x, x .* (log(x) - 1))

hold off
