We consider the weak closure of the set of all feasible pairs (solution, flow) of the family of potential elliptic systems $$\begin{array}{c}\text{div}\left(\sum _{s=1}^{s_0}{\sigma }_s\left(x\right)F_s^{\text{'}}(\nabla u\left(x\right)+g\left(x\right))-f\left(x\right)\right)=0\text{in}\Omega ,u=(u_1,\cdots ,u_m)\in H_0^1(\Omega ;{𝐑}^m),\sigma =({\sigma }_1,\cdots ,{\sigma }_{s_0})\in S,\end{array}$$ where Ω ⊂ is a bounded Lipschitz domain, are strictly convex smooth functions with quadratic growth and $S=\{\sigma measurable\mid {\sigma }_s\left(x\right)=0\text{or}1,s=1,\cdots ,s_0,{\sigma }_1\left(x\right)+\cdots +{\sigma }_{s_0}\left(x\right)=1\}$. We show that is the zero level set for an integral functional with the integrand $Qℱ$ being the -quasiconvex envelope for a certain function $ℱ$ and the operator = (curl,div). If the functions ...