Let $F(x,y) = M(x,y)i + N(x,y)j$ be a vector field on the real plane.
Suppose that Div F = $\frac{\partial M}{\partial x} + \frac{\partial N}{\partial y} = 0$
It follows that there exists a function $g: \mathbb{R}^2 \rightarrow \mathbb{R}$ (oft called a stream function) such that:
$M = \frac{\partial g}{\partial y}$ and $N = \frac{-\partial g}{\partial x}$
Now, suppose also that Curl F = $\frac{\partial M}{\partial y} - \frac{\partial N}{\partial x} = 0$
It follows that there exists a function $f: \mathbb{R}^2 \rightarrow \mathbb{R}$ (oft called a potential function) such that:
$M = \frac{\partial f}{\partial x}$ and $N = \frac{\partial f}{\partial y}$
Thus, if $F(x,y) = M(x,y)i + N(x,y)j$ has vanishing curl and divergence, we can see that:
$M = \frac{\partial f}{\partial x}=\frac{\partial g}{\partial y}$
$N = \frac{\partial f}{\partial y}=\frac{-\partial g}{\partial x}$
That is it to say, the components of $F$ are equal to the Cauchy-Riemann Equations.
This is amazing to me. I really would like to understand what's going on more. I understand that vector fields with vanishing curl and divergence are very important in mathematics. Why should we expect that such vector fields have such a direct connection with complex functions? I have a hard time knowing the exact question I want to ask because I don't know what I don't know.
I also recommend Ian's answer in this post: Is there any intuition or meaning regarding Cauchy-Riemann equations?
Where he explains the Cauchy-Riemann equations as the condition for the derivative of a complex function to act as multiplication of a complex number (and thus a linearization of the complex function).