I have two functions $P(r)$ and $Q(r)$ that can be expressed in the following power series
$$P(r) = \sum^\infty_{n=0} a_n r^n$$
$$Q(r) = \sum^\infty_{n=0} b_n r^n$$
where $r \in \mathbb R_{\ge 0}$.
The coefficents $a_n$ and $b_n$ are defined by the following recurrence relation
$$a_n = C\, n b_n - a_{n-1}$$
$$b_n = C\, n a_n + b_{n-1}$$
$$a_0 = B \,b_0$$
where $a_n, b_n, C,B \in \mathbb R$.
This recurrence relation stems from inserting the power series into a set of coupled ordinary differential equations. I performed some numerical experiments already to determine the rate and radius of convergence. However, it would be nice to have some analytical expression here. Sadly, I have great difficulties applying common convergence criterion (e.g. root test, ratio test) here because of the interdependence of the coefficients.
Is there any way to get at least an approximation for the radius of convergence?