In this paper, we consider the following nonlinear Maxwell system M\( \left \{ \textstyle\begin{array}{l@{\quad}l} \nabla \times \nabla \times E_{1}-\lambda E_{1}= \frac{\partial H}{\partial E_{2}}(x,E_{1},E_{2}) & \text{ $x\in \Omega $} \\ \nabla \times \nabla \times E_{2}-\lambda E_{2}= \frac{\partial H}{\partial E_{1}}(x,E_{1},E_{2}) & \text{ $x\in \Omega $} \\ \nu \times E_{1} =0,~~~ \nu \times E_{2} =0& \text{ $x\in \partial \Omega $} \end{array}\displaystyle \right . \) on a bounded domain $\Omega \subset \mathbb{R}^{3}$ with exterior normal $\nu :\partial \Omega \rightarrow \mathbb{R}^{3}$, $E_{1}$, $E_{2}:\Omega \rightarrow \mathbb{R}^{3}$ are vectors, $\lambda \in \mathbb{R}$ is a constant, ∇× denote the curl operator in $\mathbb{R}^{3}$, $H\in C_{0}^{1}(\Omega \times \mathbb{R}^{3}\times \mathbb{R}^{3}, \mathbb{R})$ is a real-valued subquadratic function. By utilizing a generalized Clark’s theorem, we establish the existence of infinitely many solutions for the sublinear Maxwell system (M).
