Suppose that X is distributed Poisson with a known rate and Y is a normal distributed with a know mean and variance. My goal is to approximate the distribution Z where P(Z) = P(X) * P(Y), where Z is a non-negative integer. I could get a good approximation by sampling, but I'd really like to have a fast solution, ideally closed-form.