Let $(x_n, y_n, z_n) = (3, 1, 0)$ for $n=0$
For $n \ge 1$,
$$\begin{align} x_n &= x_{n-1} +3 z_{n-1}\\ y_n &= x_{n-1} +2 z_{n-1}\\ z_n &= 5 y_{n-1} \end{align}$$
Please let me know the formula to find $x_n,y_n,z_n$ values of any integer $n$.
The following would be the series:
$$\begin{array}{c:l} n & x_n,y_n,z_n\\ \hline 0 & 3,1,0\\ 1 & 3,3,5\\ 2 & 18,13,15\\ 3 & 63,48,65\\ \end{array}$$
By finding the characteristic polynomial of the matrix $$ M = \left(\begin{matrix}1&0&3\\1&0&2\\0&5&0\end{matrix}\right) $$ that is $p(t)=t^3-t^2-10t-5$, and exploiting the Cayley-Hamilton theorem, we have that: $$ x_{n+3} = x_{n+2}+10 x_{n+1} + 5 x_n $$ and the same relation holds by replacing $x_{*}$ with $y_{*}$ or $z_{*}$. This gives an iterative method for finding $x_n,y_n,z_n$, and a closed expression too, in terms of linear combinations of the $n$-th powers of the roots of $p(t)$.