Assuming the following equation $$ \sum_{m=0}^M b_m r^m \left[\sum_{i,j}i j\lambda_i \lambda_j r^{i+j-2}-\sum_j j(j+1) \lambda_j r^{j-2}+2E\right]-2\sum_{k=0}^K a_k r^k=0 \tag{1} $$ the authors in a paper says that
Equating the coefficients of $r^p, p=0, 1, 2, . . .$ in (1) to zero, we get a set of relations connecting $\lambda_n$s. For example, in a case where M=1 and K=3, these are:
$$ \lambda_1=Z; \\ \lambda_2=(2E+\lambda^2_1)/6; \\ \lambda_3=(2\lambda_1 \lambda_2-\alpha)/6;\\ \lambda_4=(4\lambda_2^2+6\lambda_1 \lambda_3-2\beta)/20;\\ \lambda_n=\frac 1 {n(n+1)}\sum_{i=1}^{n-1}i (n-i)\lambda_i \lambda_{n-1}, n>4;\\ \tag{2} $$ In fact they have obtained an expression for this special potential in the previous section as follows (before they write $\Psi$ as a sum, see this question) $$ n(n-1)\mu_{n-1}+2Z\mu_n+2E\mu_{n+1}-2\alpha\mu_{n+2}-2\beta\mu_{n+3}=0 \tag{3} $$ where $$ \mu_j=\int_0^\infty r^j \Psi(r)dr \tag{4} $$ this is while after writing $\Psi$ as a sum they define $$ \int_0^\infty r^p exp[-S(r)]dr=\mu_p, \,\, p=0,1,2,...,N \tag{5} $$
I can't understand how they have reached equation (2) by comparison equations (1) and (3)? or if they have reached it in another way! Any idea?
You can find more details in the original paper. Of course I have cropped the relevant parts. In this question I want to know how authors have reached from equation (11) to (12) in the paper.