I want to show that if $\phi$ transforms as the adjoint representation of SU(n), its covariant derivative is given by $\textbf{D}_\mu \phi = \partial_\mu \phi + i [\textbf{A}_\mu, \phi]$. (Exercise in chapter 6 of P. Ramond QFT book)
The covariant derivative is defined as
$$\textbf{D}_\mu= \partial_\mu + i\textbf{A}_\mu$$
If $\phi$ have the transformation properties assumed, we know that $\phi= \phi^a T^a$, where $T^a$ are the generators of SU(n) and in the adjoint representation $(T^b)_{cd}= -if^{bcd} $. Besides that, we have $\textbf{A}_\mu=A_\mu^a T^a$, etc.
The attempted solution:
$$\begin{split}(\textbf{D}_\mu \phi)_{ab} & = (\partial_\mu \phi)_{ab} + i (\textbf{A}_\mu \phi)_{ab} \\ & = (\partial_\mu \phi)_{ab} + i A_\mu^c \phi^e (T^c)_{ad}(T^e)_{db} \end{split}$$
And now I can use the relation $(T^b)_{cd}= -if^{bcd} $, but it doesn't get me to the right answer because there isn't one $f$ that contracts with both $A$ and $\phi$ as that should be to get the correct result. We can see this writing the answer with explicit components
$$(\textbf{D}_\mu \phi)_{ab} = (\partial_\mu \phi)_{ab} - A_\mu^d \phi^e f^{dec}(T^c)_{ab}.$$
How should I proceed? Any help will be appreciated.
$\textbf{Edit}$
Following the comments, I got to a point and now I am stuck:
$$\begin{split}(\textbf{D}_\mu \phi)_{be} & = (\partial_\mu \phi)_{be} + i A_\mu^a \phi^d (T^a)_{bc}(T^d)_{ce} \\ & = (\partial_\mu \phi)_{be} + i A_\mu^a \phi^d (-if^{abc})(-if^{dce}) \\ & = (\partial_\mu \phi)_{be} + i A_\mu^a \phi^d (-f^{adc}f^{ceb}-f^{aec}f^{cbd}) \\ & = (\partial_\mu \phi)_{be} - A_\mu^a \phi^d f^{adc}(T^c)_{be}-A_\mu^a \phi^d f^{aec}f^{cbd} \\ & =(\partial_\mu \phi)_{be} +i A_\mu^a \phi^d [T^a,T^d]_{be}-A_\mu^a \phi^d f^{aec}f^{cbd}\end{split} $$
Which shows that I got the right answer plus $-A_\mu^a \phi^d f^{aec}f^{cbd}$. Is this right so far?