Solution (source code)

= Solution

The same fitted <factorial> approximation gives
$$
\boxed{\ln(1000!)\approx L(1000)\simeq5912.1174.}
$$