I am interested in finding an expression (closed form or recursive) for the matrix exponential of this banded matrix:
$$ \begin{pmatrix} 0 & 1 & 0 & 0 & \cdots & 0 & 0 \\ 0 & a_1 & 1 & 0 & \cdots & 0 & 0 \\ 0 & 0 & a_2 & 1 & \cdots & 0 & 0 \\ 0 & 0 & 0 & a_3 & \cdots & 0 & 0 \\ \vdots & \vdots & \vdots & \vdots &\ddots&\vdots&\vdots \\ 0 & 0 & 0 & 0 & \cdots & a_{n-1} & 1 \\ 0 & 0 & 0 & 0 & \cdots & 0 & a_n \end{pmatrix} $$
For simplicity, assume that $a_k>0~\forall k$. I am quite certain one must exist having played around with it for a while, and looking at the solution for small values of $n$.
Has anyone seen this structure before? Does it have a name? Do you know if there is a solution published somewhere?
If you go ahead and compute the answer for small values of $n$, you get:
$$ \exp\left( \begin{array}{cc} 0 & 1 \\ 0 & a_1 \\ \end{array} \right) = \left( \begin{array}{cc} 1 & \frac{-1+e^{a_1}}{a_1} \\ 0 & e^{a_1} \\ \end{array} \right) $$
$$ \exp \left( \begin{array}{ccc} 0 & 1 & 0 \\ 0 & a_1 & 1 \\ 0 & 0 & a_2 \\ \end{array} \right) = \left( \begin{array}{ccc} 1 & \frac{-1+e^{a_1}}{a_1} & \frac{-e^{a_2} a_1+a_1+e^{a_1} a_2-a_2}{a_1 \left(a_1-a_2\right) a_2} \\ 0 & e^{a_1} & \frac{e^{a_1}-e^{a_2}}{a_1-a_2} \\ 0 & 0 & e^{a_2} \\ \end{array} \right) $$
Unfortunately, $n=3$ is to large to print here, but a pattern remains.