I'm looking for an upper bound to the following integral or equivalent when $n$ leads to $ +\infty $ to the following expression
$$I_n:=\left|\int_{0}^1 \int_{0}^1 \frac{p_n(x) p_n(y)}{(1-xy)} dx dy \right| $$
with
$$ p_n(t):=\frac{1}{n!}(t^n(1-t)^n)^{(n)}.$$
This integral is similar to Beukers integral; after integrating $n$ times to $y$ ,I obtain
$$I_n:=\left|\int_{0}^1 \int_{0}^1 \frac{p_n(x) x^n y^n(1-y)^n}{(1-xy)^{(n+1)}} dx dy \right| ,$$$$I_n:=\left|\int_{0}^1 \int_{0}^1 \frac{p_n(x) x^n y^n(1-y)^n}{(1-xy)^{n+1}} dx dy \right| ,$$
I don't know what to do after this since I can't have interesting expression inside the integral when I derive n times $x$ .