Given U-Unif(0,1) and U=u and X-Pois(u) compute E[X].
I could think of a way to compute this but it seems fairly convoluted.
$P(X=k|U=u)=u^k e^{-u}/ k!$ and using the law of total probability (? I guess not sure what exactly this is using the pdf of U seems weird but I saw it in https://math.stackexchange.com/a/567234)
$P(X=k)=\int_0^1 P(X=k|U=u)f_U(u) du=\int_0^1 u^k e^{-u} / k! du $
And lastly sum over this but this integral seems ugly to compute out it looks like an incomplete gamma function but the problem I was solving expects a numerical answer and im not sure how to arrive at that through this approach.
The MSE post mentions using function transformation (MGF, CF etc) but again not clear where that would enter the picture here.
Would appreciate any general advice for problems like these

. They now have 845