Is it possible to derive Maxwell's laws ... only from that formula ... I think it is interesting if we can derive all classic electromagnetic field by a single formula ...
Hehl & Obulov provide a derivation of Maxwells equations starting only from relativistic charge conservation over any spacetime, this being a Lorentzian manifold. This principle is expressed in a 'single formula'. The usually equations are found after imposing a time field. This is detailed in their book, Foundations of Electromagnetism.
In a short paper on the Arxiv, they detail a similar derivation but over space. Here, they require both electric charge and magnetic flux conservation. This can be deduced from the preceding, as relativstic charge conservation implies electric charge and magnetic flux conservation. But here, they consider only space as opposed to time and hence they require two independent principles.
Their equations do not require a metric, but in this case they require four fields. When a spatial metric is imposed, the fields are reduced to two and the usual Maxwell equations are recovered.
The latter derivation is sketched out as follows. It uses the technology of differential forms but I use the term 'cofield' as a synonym for differential form as it's more closely aligned with the terminology in physics, aka fields. Thus a k-cofield is a differential k-form. In particular, a scalar cofield is a $0$-cofield and a coscalar cofield is a $0'$-cofield where $k':= m-k$.
Now, let $M$ be a spatial manifold, by this we mean a Riemannian manifold describing space. We let $m$ be it's dimension.
Let $\rho$ be the charge codensity. We call this a codensity because we describe it by by a coscalar cofield rather than a scalar cofield. This is because we want to avoid the use of the metric; and the Hodge star, which involves the metric, is required to transform a charge density into a charge codensity and it is the codensity rather than the density we integrate over a body to find its total charge (in the presence of a metric, the density is enough).
Now choose any closed ball $B$ of space $M$, then the total charge contained by $B$ is:
$Q:= \int_B \rho$
Now, since $\rho$ is a top cofield, it's exterior differential vanishes and hence it is a closed cofield. By Poincare's lemma, it will have a local potential $D$ which we call the electrical excitation. Thus:
$dD = \rho$, locally
This is Gauss's electric law. Notice it is only true locally. Now, global charge conservation means that we have current $j$ which is a cofield of degree $m-1$ and it satisfies:
$\partial_t \int_B *\rho + \int_{\partial B} j = 0$
We can rewrite this using Stoke's theorem as:
$\partial_t \int_B *\rho + \int_B dj = 0$
Since this is true for any closed ball $B$, we obtain the local form of charge conservation:
$\partial_t *\rho + dj = 0$
Locally, this is:
$\partial_t dD + dj = 0$
And this is equal to:
$d(\partial_t D + j)= 0$
The term in the brackets is again a closed form, so we can again introduce locally a potential for this term and which we call the magnetic excitation $H$. So:
$dH = \partial_t D + j$
This gives the Ampere-Maxwell law. Thus we have derived the inhomogeneous Maxwell's equations purely from electric charge conservation and without the use of any metric and all this over any smooth orientated manifold. However, to write it in traditional terms through curl and div, we would require a metric. First we introduce their cofield variants of div and curl:
$cocurl = *d$ and $codiv = *d*$
These are related to the traditional versions through the raising and lowering operators:
$curl = \sharp \circ cocurl \circ \flat$ and $div = \sharp \circ codiv \circ \flat$
And also the electric and magnetic cofields $E$ and $B$ are defined by what Hehl & Obukov call the constitutative equations:
$D := \epsilon_0. *E$ and $\mu_0.H:= *B$
Recall, Gauss's electric law is $dD = \rho$, locally. So replacing the electric excitation $D$ with the electric field we obtain, $d*E=\rho/\epsilon_0$ and then applying the Hodge star, we further get $(*d*)E =*\rho/\epsilon_0$. This is just $codiv(E) = *\rho/\epsilon_0$. Then introducing the raising and lowering operators, we get $\sharp \circ codiv \circ \flat \circ \sharp(E) = \sharp(*\rho)/\epsilon_0$. And this reduces to:
$div (E^{\sharp}) = (*\rho)^{\sharp}/\epsilon_0$
This is the traditional electric Gauss equation once we recognise as we introduced $E$ as an electric cofield, then $E^{\sharp}$ is the electric field. And as we introduced $\rho$ as the electric charge codensity which is a coscalar field, then $*\rho$ is the electric charge density which is a scalar field and as the raising operator, $\sharp$, does nothing on scalar cofields, ie they are equivalently scalar fields. Thus renaming $E^{\sharp} \rightarrow E$ and $(*\rho)^{\sharp} \rightarrow \rho$, we get the traditional equation:
$div(E) = \rho/\epsilon_0$, local electric Gauss law.
But unlike the traditional equation this is valid on any Riemannian manifold considered as the spatial manifold. Similarly, for the Ampere-Maxwell law.
Hehl & Obukhov derive the homogeneous equations from the conservation of magnetic flux. Again, this does not require a metric. Instead, the electric and magnetic cofields are introduced independently and the conservation of magnetic flux is used to derive Gauss's magnetic law and Faraday's law. To find the traditional forms again we apply the constitutive laws as above.