Approximating a continuous distribution with discrete points I want to estimate the pdf distribution of a polynomial transformation of a known continuous random variable $X$:
$Y = f(X) = a_nX^n + \cdots + a_1X + a_0$
Normally, the only numerical approach I know of is by using a simulation.
What I tried to do:
evaluate discrete points the pdf of $Y$ by solving a root finding problem:
$P(a_nX^n + \cdots + a_1X + a_0 = c) = P(X = \alpha)$,
where $\alpha$ is the root of  $a_nX^n + \cdots + a_1X + a_0 - c = 0$.
I do this for varying $c$.  I then evaluate these $\alpha$'s in the continuous distribution of the known $X$, and use the found points to fit a curve as approximation to the pdf of $Y$.
When I try this, the shape of the curve is right, however, it is not normalized.
I believe I am making an error by evaluating discrete points in a continuous distribution (which should be 0 by definition), but why is it that the shape of the curve is still right?
 A: Your problem is not really related to discretization, but with a bad understanding of changes of variables. It's true that the probability function of  discrete variables gives you directly the probability of the events $P(X=x)$, and hence it's enough to find the inverse mapping; but for continuous variables things are a little more subtle. A continuous density (say, for the sake of simplicity a continuous random variable X which has besides a continuous density) gives you (informally) the probability that the variable falls inside a neighbourhood of $x$, multiplied by the width of that interval:
$$f_X(x) \; \Delta x \approx P(x-\frac{\Delta x}{2} \leq X < x+\frac{\Delta x}{2}) \; $$
(the approximation is to be understood as a limit).
Now, if you have a transformation $Y = g(X)$, to get the corresponding $f_Y(y)$, then, it's not enough to find out the inverse mapping, $x = g^{-1}(y)$ (analytically or numerically... and leaving aside temporarily the issue of multiple solutions) and then write $f_Y(y) = f_X(g^{-1}(y))$. That's basically what you've apparently done, and that is wrong, because the events "X equals $x$" cannot be assigned a probability. The correct way would be to put in correspondence the events "X falls in this neighborhood of $x$" and  "Y falls in this neighborhood of $y$" and take the limit:
$$ f_Y(y) \; \Delta y   \approx f_X(x) \; \Delta x $$
$$ f_Y(y) \approx \frac{f_X(g^{-1}(y))}{\frac{\Delta y}{\Delta x}}   $$
If $g(X)$ is monotonous incresing, the denominator tends to $g'(x)=g'(g^{-1}(y))$ , and you get the standard formula for change of variables. If not, you have a little more work to do (as is your case) - but at least you get the idea.
