I want to compute $\mathbb{E}\left[e^{\lambda X_t} X_s^2 \right]$ with $X_t= \mu t + W_t$.
Also $0 0$ and $W_t$ is a standard Brownian Motion.
To do this I first computed \begin{align*} \mathbb{E}\left[e^{\lambda X_t}X_s^2\right] &= e^{\lambda \mu (t-s)} \mathbb{E}\left[e^{\lambda (W_t-W_s)} \right] \mathbb{E}\left[e^{\lambda X_s} X_s^2\right] \\ &=e^{\lambda \mu (t-s)} \mathbb{E}\left[e^{\lambda (W_t-W_s)} \right] \mathbb{E}\left[e^{\lambda X_s}\right] e^{\lambda \mu s} \mathbb{E}\left[ e^{\lambda W_s} \left(W_s^2 +2 \mu s W_s + \mu^2 s^2 \right) \right] \end{align*} Now since $W_t-W_s \sim N(0,t-s)$ we can substitute this in $\mathbb{E}\left[e^{\lambda (W_t-W_s)} \right]$ to compute this using integral manipulation.
Can I also do this for $\mathbb{E}\left[ e^{\lambda W_s} \left(W_s^2 +2 \mu s W_s + \mu^2 s^2 \right) \right]$? (I know I can split this into 3 parts by the linearity of the expectation)
Or is there another easier way of evaluating this expectation?