Deje $\beta>0$ ser dado. Para cada una de las $n\geq 2$, vamos a $\Delta_n=\det M_n$ el valor del determinante de la siguiente matriz: \begin{align} M_n = \begin{pmatrix} 2+\epsilon^2 & -1 & 0 & 0 & 0 & -1 \\ -1 & 2+\epsilon^2 & -1 & 0 & 0 & \ddots \\ 0 & -1 & 2+\epsilon^2 & -1 & 0 & \ddots \\ 0 & 0 & -1 & 2+\epsilon^2 & -1 & \ddots \\ 0 & 0 & 0 & -1 & 2+\epsilon^2 & \ddots \\ -1 & \ddots & \ddots & \ddots & \ddots & \ddots \\ \end{pmatrix}, \qquad \epsilon = \frac{\beta}{n} \end{align} ¿Cómo se puede evaluar el límite \begin{align} \Delta = \lim_{n\to\infty}\Delta_n ? \end{align} Este problema surgió en el contexto de la evaluación de un determinado camino integral en un problema de física. Sé cómo determinar los valores propios de cada una de las $M_n$, pero calcular el determinante por tomar el producto de los autovalores conduce a un producto que no se puede evaluar. También probé la que se derive una relación de recursividad de la $\Delta_n$ y demostrando que en el límite de $n\to\infty$, la recursividad relación puede ser considerado como una ecuación diferencial cuya solución sujeto a ciertos datos iniciales determina $\Delta$, pero no así.
Cualquier conocimiento se agradece.
Adenda.
Cuando me calcular los autovalores $\lambda_k$$M_n$, puedo obtener \begin{align} \lambda_k = 4\sin^2\left(\frac{\pi k}{n}\right)+\left(\frac{\beta}{n}\right)^2, \qquad k=0,1,\dots, n-1 \end{align} a partir de la cual se deduce que \begin{align} \Delta_n = \prod_{k=0}^{n-1}\left[4\sin^2\left(\frac{\pi k}{n}\right)+\left(\frac{\beta}{n}\right)^2\right] \end{align} que no tengo la más remota idea de cómo evaluar. Yo, sin embargo, tienen una conjetura para la respuesta que viene del hecho de que este límite de determinantes vino de una ruta integral que puede ser calculado en otras formas. Mi conjetura (que he comprobado numéricamente usando mathematica hasta cierto punto) es \begin{align} \Delta = 4\sinh^2\frac{\beta}{2} \end{align}