I was just playing around in desmos and I think I found something that approximates Lambert W for whole numbers.
$$f(x)=a(\int_{0}^{1}(\sum_{n=1}^{x}t^n-1)dt)+b$$
Where $a\approx0.765424$ And $b\approx0.944602$
$a$ and $b$ were achieved through linearization, plugging $xe^x$ in for x.
Please point me in the right direction for proofing or disproving that as $x \to \infty$ $f(x) \to W(x)$