 Solution

Split the microscopic sources into free and bound parts and defineSubstitution of $\rho ={\rho }_{\mathrm{f}}-\nabla \cdot P$ into Gauss's law gives $\nabla \cdot D={\rho }_{\mathrm{f}}$. Similarly, using $J=J_{\mathrm{f}}+{\partial }_{t}P+\nabla \times M$ in the Ampère-Maxwell equation gives $\nabla \times H=J_{\mathrm{f}}+{\partial }_{t}D$. The homogeneous Maxwell equations are unchanged. Thus