Solution (source code)

= Solution

Direct expansion gives
$$
[D_\mu,D_\nu]\phi
=\{\partial_\mu A_\nu-\partial_\nu A_\mu+[A_\mu,A_\nu]\}\phi
=F_{\mu\nu}\phi.
$$
Since $D'_\mu(g\phi)=gD_\mu\phi$,
$$
F'_{\mu\nu}(g\phi)
=[D'_\mu,D'_\nu](g\phi)
=g[D_\mu,D_\nu]\phi
=gF_{\mu\nu}\phi.
$$
Therefore the non-Abelian <gauge field strength> transforms covariantly:
$$
F'_{\mu\nu}=gF_{\mu\nu}g^{-1}.
$$