Solution (source code)

= Solution

The <Lorenz gauge> is $\partial_\mu A^\mu=0$. Substituting the <electromagnetic four-potential> into the covariant <Maxwell equations> then removes the mixed derivative and gives
$$
\boxed{\Box A^\mu=\mu_0J^\mu},
\qquad
\Box=\frac1{c^2}\partial_t^2-\nabla^2.
$$
Equivalently,
$$
\left(\nabla^2-\frac1{c^2}\partial_t^2\right)A^\mu=-\mu_0J^\mu.
$$