A general way to solve recurrences is generating functions. Define:
$\begin{align*}
A(z)
&= \sum_{n \ge 0} a_n z^n
\end{align*}$
Shift your recurrence, multiply by $z^n$, sum over $n \ge 0$, recognize the resulting sums:
$\begin{align*}
\sum_{k \ge 0} a_{k + 1} z^k
&= \sum_{k \ge 0} a_k z^k
+ 3 \sum_{k \ge 0} (k + 1) z^k + \sum_{k \ge 0} z^k \\
\frac{A(z) - a_0}{z}
&= A(z) + 3 \frac{z}{(1 - z)^2} + \frac{4}{1 - z}
\end{align*}$
Here we used:
$\begin{align*}
\frac{1}{1 - z}
&= \sum_{k \ge 0} z^k \\
\sum_{k \ge 0} k z^k
&= z \frac{d}{d z} \sum_{k \ge 0} z^k \\
&= z \frac{d}{d z} \frac{1}{1 - z} \\
&= \frac{z}{(1 - z)^2}
\end{align*}$
Solve for $A(z)$, as partial fractions:
$\begin{align*}
A(z)
&= \frac{z (4 - z)}{(1 - z)^3} \\
&= \frac{3}{(1 - z)^3} - \frac{2}{(1 - z)^2} - \frac{1}{1 - z}
\end{align*}$
To extract coefficients, remember:
$\begin{align*}
(1 + u)^{-m}
&= \sum_{k \ge 0} \binom{-m}{k} u^k \\
&= \sum_{k \ge 0} (-1)^k \binom{k + m - 1}{m - 1} u^k
\end{align*}$
Thus:
$\begin{align*}
[z^k] A(z)
&= 3 (-1)^k \binom{k + 3 - 1}{3 - 1} (-1)^k
- 2 (-1)^k \binom{k + 2 - 1}{2 - 1} (-1)^k
- (-1)^k \binom{k + 1 - 1}{1 - 1} (-1)^k \\
&= 3 \binom{k + 2}{2} - 2 \binom{k + 1}{1} - \binom{k}{0} \\
&= 3 \frac{(k + 2) (k + 1)}{2} - 5 \frac{k + 1}{1} + 2 \\
&= \frac{k (3 k + 5)}{2}
\end{align*}$