Tenemos una función $F(\beta,\lambda)$ (polinomio de grado $n$) que es analítica cerca de $\beta_0$ y $\lambda_0$. Por lo tanto, podemos escribir $$F(\beta,\lambda)=\sum_{m=0}^n(\lambda-\lambda_0)^mf_m(\beta)$$ donde $f_0(\beta_0)=F(\beta_0,\lambda_0)=0$ y $f_1(\beta_0)=\frac{\partial F}{\partial\lambda}(\beta_0,\lambda_0)\neq0$. Queremos resolver la ecuación $F(\beta,\lambda)=0$, que es equivalente a \begin{equation} \lambda=\lambda_0-\frac{f_0(\beta)}{f_1(\beta)}-\sum_{m=2}^n(\lambda-\lambda_0)^m\frac{f_m(\beta)}{f_1(\beta)} \end{equation} Intentamos resolver esta ecuación con una solución en la forma $$\lambda(\beta)=\sum_{k=1}^\infty\alpha_k(\beta-\beta_0)^k$$ Mediante el uso de la sustitución recursiva en la última ecuación podemos calcular los $\alpha_k$'s $$\alpha_1=-\left[\frac{f_0(\beta)}{f_1(\beta)}\right]'\bigg|_{\beta=\beta_0}\ ,\quad \alpha_2=-\frac{1}{2}\left[\frac{f_0(\beta)}{f_1(\beta)}\right]''\bigg|_{\beta=\beta_0}-\alpha_1^2\frac{f_2(\beta_0)}{f_1(\beta_0)}$$
Realmente no veo cómo utilizamos la sustitución recursiva y de dónde provienen las fórmulas de $\alpha_1$ y $\alpha_2$. Cuando comparo las dos ecuaciones y expando en serie de Taylor $f_0(\beta)$ y $f_1(\beta)$ cerca de $\beta_0$, termino con todas las derivadas superiores de la expansión. Cualquier ayuda sería apreciada.