The following equations are equations of electrostatics: $$\nabla \times \vec E=0$$ $$\nabla\cdot\vec E=\dfrac{\rho}{\epsilon_0}.$$
These are 4 independent equations, while $\vec E$ has only 3 independent components. Yet these equations do not completely specify the field, as after adding the gradient of a scalar $\nabla \lambda$ that satisfies Laplace equation ($\nabla^2 \lambda$=0) to $\vec E$ leaves the equations unchanged: $$\cases{\nabla \times \vec E=0\\\nabla\cdot\vec E=\dfrac{\rho}{\epsilon_0}}\xrightarrow[\nabla^2\lambda=0]{\vec E'=\vec E+\nabla \lambda}\cases{\nabla \times \vec {E'}=0\\\nabla\cdot\vec {E'}=\dfrac{\rho}{\epsilon_0}}$$ (note the primes on the RightHandSide $\vec E$s)
The system should be overdetermined (4 equations, 3 unknowns) but apparently it is underdetermined.
- Is the system overdetermined or underdetermined?
- How do we usually choose the arbitrary $\lambda$ in a problem with $\rho$ given and (say) Neumann boundary condition?
- Why the first equation (curl) is not enough?