You should first probably not include $n$ in your polynomial. You just assume that $|p(x)-f(x)|<\epsilon$. You then use the fact that $|e^{-nx/(1-x)}|\le 1$ so (regardless of $n$):
$$\left|\int_0^1e^{-nx/(x-1)}p(x)dx\right| < \epsilon$$
Then you set $q(x-1) = p(x)$ integrate by parts:
$$\int_0^1e^{-nx/(1-x)}q(x-1)dx = \left[e^{-nx/(1-x)}Q(x-1)\right]_0^1 + \int_0^1{ne^{-nx/(1-x)}Q(x-1)\over (x-1)^2}dx$$
If we know $q(0) = 0$ (which means it's $x^0$ term is zero) and $Q$ is any antiderivate of $q$ so we can choose it so that the constant term is zero we see that $Q(x-1)/(x-1)^2$ is again a polynomial. And we can continue to integrate by parts. This first step we get:
$$\int_0^1e^{-nx/(1-x)}q(x-1)dx = -Q(-1) + \int_0^1{ne^{-nx/(1-x)}Q(x-1)\over (x-1)^2}dx$$
What we then want is the last integral to converge to zero which would prove that $|Q(-1)| = |P(0)| < \epsilon$. Repeating this results in an estimate of the terms of $p(x)$.
Now one would of course need to compute the result of the repeated integration by parts so that this results in an uniform estimate of $p(x)$ only in terms of $\epsilon$. We know it's less than the absolute sums of it's coefficient, so we need to make sure that the estimate don't depend on the degree of $p$.
Now we used that $q(0)=0$, to handle the case where this isn't we have to separate out it from the polynomial. We would solve this by choosing a polynomial on the form $C+q(x)$ instead where $q(0)=0$ and we would get an extra term in the integration by parts:
$$C\int_0^1 e^{-nx/(1-x)}$$
and we would like this to converge to zero as well. The integrand converges uniformly on any compact subset of the interval $[0,1]$ and it's otherwise bounded. That is:
$$C\int_0^1 e^{-nx/(1-x)} = C\int_0^\eta e^{-nx/(1-x)}dx + C\int_\eta^1 e^{-nx/(1-x)}dx $$
where the second term aproaches zero and the absolute of the first term is no larger than $C\eta$.