In the general model of a Pólya urn, we consider the balanced urn model for which $b = c = 0$, so $a = d$. Each time we draw a ball, we add $a$ balls of its color to the urn. The initial composition is $a_{0}$ white balls and $b_{0}$ black balls. The function $G$ is defined on $U = \{(x,u,v) \in \mathbb{R} \times \mathbb{R}_{+}^{*} \times \mathbb{R}_{+}^{*} ; axu^{a} < 1, axv^{a} < 1\}$ by $$G(x,u,v) = u^{a_{0}} v^{b_{0}} (1 - axu^{a})^{-a_{0}/a} (1 - axv^{a})^{-b_{0}/a}$$ We use the notation $D_{\rho} = ]-\rho, \rho[ \times ]0,2[^{2}$.
Using the preliminary results, prove that there exists $\rho > 0$ such that $D_{\rho} \subset U$ and, for all $(x,u,v) \in D_{\rho}$, $$G(x,u,v) = \sum_{n=0}^{+\infty} Q_{n}(u,v) \frac{x^{n}}{n!}$$ where $Q_{n}$ is a polynomial function of two variables to be specified.
In the general model of a Pólya urn, we consider the balanced urn model for which $b = c = 0$, so $a = d$. Each time we draw a ball, we add $a$ balls of its color to the urn. The initial composition is $a_{0}$ white balls and $b_{0}$ black balls. The function $G$ is defined on $U = \{(x,u,v) \in \mathbb{R} \times \mathbb{R}_{+}^{*} \times \mathbb{R}_{+}^{*} ; axu^{a} < 1, axv^{a} < 1\}$ by
$$G(x,u,v) = u^{a_{0}} v^{b_{0}} (1 - axu^{a})^{-a_{0}/a} (1 - axv^{a})^{-b_{0}/a}$$
We use the notation $D_{\rho} = ]-\rho, \rho[ \times ]0,2[^{2}$.

Using the preliminary results, prove that there exists $\rho > 0$ such that $D_{\rho} \subset U$ and, for all $(x,u,v) \in D_{\rho}$,
$$G(x,u,v) = \sum_{n=0}^{+\infty} Q_{n}(u,v) \frac{x^{n}}{n!}$$
where $Q_{n}$ is a polynomial function of two variables to be specified.