Let $$ M = \begin{bmatrix} d_{1} - \lambda & 0 & \cdots & 0 & \zeta_{1} \\ 0& d_{2} - \lambda & \cdots & 0 & \zeta_{2} \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & d_{n-1} - \lambda & \zeta_{n-1} \\ \zeta_{1} & \zeta_{2} & \cdots & \zeta_{n-1} & \alpha - \lambda \end{bmatrix}. $$ Compute determinant of $M$.
So, my guess was to solve this using Laplace expansion, but I got stuck. Is there any faster and easier way to solve this, since $M$ is symmetric? Any hint helps!