If $\displaystyle f^{(n)}=a_{n-1}f^{(n-1)}+\cdots+a_1f'+a_0$ then $$\frac{d}{dt} \begin{pmatrix}f \\ f' \\ \vdots \\ f^{(n-2)} \\ f^{(n-1)}\end{pmatrix}=\begin{pmatrix} 0 & 1 & 0 & \cdots & 0 \\ 0 & 0 & 1 & \cdots & 0 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\
0 & 0 & 0 & \dots & 1 \\ a_0 & a_1 & a_2 & \cdots & a_{n-1} \end{pmatrix} \begin{pmatrix}f \\ f' \\ \vdots \\ f^{(n-2)} \\ f^{(n-1)}\end{pmatrix}.$$ Writing the column vector above as $\vec{f}$ and the constant matrix as $A$, we have $\vec{f}^{(k)}=A^k\vec{f}$, and we also have that $f^{(n+k)}$ is the last component of $\vec{f}^…
0 & 0 & 0 & \dots & 1 \\ a_0 & a_1 & a_2 & \cdots & a_{n-1} \end{pmatrix} \begin{pmatrix}f \\ f' \\ \vdots \\ f^{(n-2)} \\ f^{(n-1)}\end{pmatrix}.$$ Writing the column vector above as $\vec{f}$ and the constant matrix as $A$, we have $\vec{f}^{(k)}=A^k\vec{f}$, and we also have that $f^{(n+k)}$ is the last component of $\vec{f}^…