Maxwell's equations

$$\triangledown \times \mathbf{E} = -\frac{\partial \mathbf{H}}{\partial t}; \, \,\, \triangledown \times \mathbf{H} = \frac{\partial \mathbf{E}}{\partial t};$$

in 2D with starting condition $\mathbf{H} = (0, 0, H_z)$ have the following form:

$$\frac{\partial E_x}{\partial t} = \frac{\partial H_z}{\partial y}; \,\,\, \frac{\partial E_y}{\partial t} = -\frac{\partial H_z}{\partial x}; \,\,\, \frac{\partial H_z}{\partial t} = \frac{\partial E_x}{\partial y} - \frac{\partial E_y}{\partial x};$$

Example of numerical simulations with periodic boundary conditions: