Suppose $u\sim N(0,I_p)$ and $Y|U\sim N(x(t),\sigma_e^2I_m)$, and the marginal distribution of $y$ is $f(y)=\int f(y|u)f(u)du$.
$x(t)$ is composite function of $u$. The problem is I need to generate random variable from $f(y)$, I am reading Monte Carlo method textbook, but there is no detail about it. The textbook is giving an example of normal mixture and the example is given Generating random variables from a mixture of Normal distributions here too, but it's not similar what I am dealing with. The fact is how do I assume the distribution(marginal distribution) of $y$. Or do I need to go for indirect method (eg. accept-reject method)?