Eigenproblem for Jacobi matrices: hypergeometric series solution

We study the perturbative power series expansions of the eigenvalues and eigenvectors of a general tridiagonal (Jacobi) matrix of dimension d. The (small) expansion parameters are the entries of the two diagonals of length d-1 sandwiching the principal diagonal that gives the unperturbed spectrum. The solution is found explicitly in terms of multivariable (Horn-type) hypergeometric series in 3d-5 variables in the generic case. To derive the result, we first rewrite the spectral problem for the Jacobi matrix as an equivalent system of algebraic equations, which are then solved by the application of the multivariable Lagrange inversion formula. The corresponding Jacobi determinant is calculated explicitly. Explicit formulae are also found for any monomial composed of eigenvector's components.