Maxwell's equations of classical electrodynamics are, in vacuo,
$$
\begin{array}{ll}
\nabla \cdot \mathbf{E}=\rho, & \nabla \times \mathbf{E}+\frac{\partial \mathbf{B}}{\partial t}=0 \\
\nabla \cdot \mathbf{B}=0, & \nabla \times \mathbf{B}-\frac{\partial \mathbf{E}}{\partial t}=\mathbf{j}
\end{array}
$$
(where we are using Heaviside-Lorentz rationalized units, see Appendix C of Aitchison and Hey). Show that these equations are equivalent to the following covariant equation for $A^{\mu}$ :
$$
\square^{2} A^{\mu}-\partial^{\mu}\left(\partial_{\nu} A^{\nu}\right)=j^{\mu},
$$
with $j^{\mu}=(\rho, \mathbf{j})$, and where $A^{\mu}=(\phi, \mathbf{A})$, the four-vector potential, is related to the electric and magnetic fields by
$$
\mathbf{E}=-\frac{\partial \mathbf{A}}{\partial t}-\nabla \phi, \quad \mathbf{B}=\nabla \times \mathbf{A}
$$