use of polarization and magnetization , separately, to describe the response of a medium to electromagnetic field is useful for media which are not spatially dispersive or which have a particular form of spatial dispersion [@24]. There are in general media for which it is not clear how the separation of the total current induced in medium is to be made into electric and magnetic parts. The disturbances, caused by the electromagnetic field, in such media are described completely in terms of only a response equation relating the total current induced in medium to the electric field [@24]. The quantization method outlined in the previous sections can cover such cases applying the Lagrangian $$\begin{aligned}
\label{c58}
 L(t)&=& \int _0^\infty d\omega \int d^3r \left[ \frac{1}{2}
 \dot{\vec{X}}_\omega \cdot  \dot{\vec{X}}_\omega
 -\frac{1}{2}\omega^2
 \vec{X}_\omega \cdot  \vec{X}_\omega\right]\nonumber\\
&+&\int d^3r \left[ \frac{1}{2}\varepsilon _0
\vec{E}^2-\frac{\vec{B}^2}{2\mu_0}\right]\nonumber\\
&+& \int _0^\infty d \omega \int d^3 r \int d^3 r' f_{ij}( \omega ,
\vec{r} , \vec{r'}) E^i(\vec{r},t) \vec{X}_\omega ^j(\vec{r'} ,t)\end{aligned}$$ which is the same as the Lagrangian (\[c1\])-(\[c4\]) with the exceptional that, now the medium is modeled by a single set of three dimensional harmonic oscillators $ \vec{X}_\omega$. In this case the Euler- Lagrange equation (\[c11\]) takes the form $$\begin{aligned}
\label{c59}
&& \frac{d}{d 