Let $A\_n\in\mathbb{R}^{n\times n}$ be defined as
$$
A\_n=\begin{bmatrix} a & b & 0 & \cdots & \cdots & 0 & 0\\ b & a & b & \cdots & \cdots & 0 & 0\\ 0 & b & a & \cdots & \cdots & 0 & 0\\ \vdots & \vdots & \vdots &\ddots & \ddots & \vdots & \vdots \\ \vdots & \vdots & \vdots &\ddots & \ddots & \vdots & \vdots \\ 0 & 0 & 0 &\cdots & \cdots & a & b \\ 0 & 0 & 0 &\cdots & \cdots & b & a\end{bmatrix},
$$
where $a,b\in\mathbb{R}$. It is well-known that the eigenvalues of $A\_n$ are
$$
\text{eig}(A\_n) =\left\{ a+2b\cos\left(\frac{\pi}{n+1}k\right), \ k=1,2,\dots,n \right\}.
$$



> 
> **My question.** Does there exist a closed-form expression for the eigenvalue density of the sequence $\{A\_n\}$ as $n\to \infty$?
> 
> 
> 



Eigenvalue density of a symmetric tridiagonal matrix
For large $n$ we may treat $x\equiv k/n+1$ as a continuous variable with a uniform density in the interval $0 <x< 1$. The corresponding eigenvalue $\lambda(x)=a+2b\cos\pi x$ ranges from $a-2|b|$ to $a+2|b|$. Since
$$|d\lambda/dx|=\pi\sqrt{4b^2-(\lambda-a)^2}.$$
The eigenvalue density follows from
$$\rho(\lambda)d\lambda=ndx\Rightarrow \rho(\lambda)=\frac{n}{\pi}\frac{1}{\sqrt{4b^2-(\lambda-a)^2}},\;\;a-2|b|<x<a+2|b|.$$
As a check $\int \rho(\lambda)d\lambda=n$.


