complex geometry - Newlander-Nirenberg for surfaces

Quite a long ago, I tried to work out explicitly the content of the Newlander-Nirenberg theorem. My aim was trying to understand wether a direct proof could work in the simplest possible case, namely that of surfaces. The result is that the most explicit statement I could get is a PDE I was not able to solve.

Assume a quasi-complex structure $J$ is given on the surface $S$; we want to prove that this is induced by a complex one (in this case there are no compatibility conditions). This can be easily transformed in the problem of local existence for a second order PDE, as follows.

We look for local charts on $S$ which are holomorphic (with respect to the quasi-complex structure on $S$). Two such charts are then automatically compatible. So the problem is local.

Fix a small open set $U subset S$ and identify it with a neighboorhood of $0 in mathbb{R}^2$ via a differentiable chart. Locally we can write $J = left(begin{matrix}[a | b] \ [c | d ]end{matrix}right)$ for some functions $a, cdots, d$ (pretend it is a two by two matrix).

A chart is given by a complex valued function $f = u + iv$. The condition that the differential is $mathbb{C}$-linear can be verified on a basis of the tangent space; moreover if it is true for a vector v, it remains true for Jv, which is linearly independent. Here we have used that $J^2 = -1$.

So we need only to check it for the vector $partial_x$. Since $J partial_x = a partial_x + c partial_y$, the condition says

$-v_x = a u_x + c u_y$

$u_x = a v_x + c v_y$

Hence we need to solve this system, with $f = u + i v$ non singular ($f$ will be then locally invertible). Since $a$ and $c$ do not vanish simultaneously, we can assume $c(0) neq 0$, hence $c neq 0$ on $U$ provided $U$ is small.

We can then solve for $u_y$ and get the equivalent system

$u_x = a v_x + c v_y$

$-u_y = frac{1 + a^2}{c}v_x +a v_y$

Moreover the Jacobian $J_f = u_x v_y + u_y v_x = frac{1}{c}(v_x^2 + (a v_x + c v_y)^2)$, so $f$ is nonsigular if $v$ is. By Poincaré's lemma, the system admits a local solution if and only if

$frac{partial}{partial_y} left( a v_x + c v_y right) - frac{partial}{partial_x} left( frac{1 + a^2}{c}v_x +a v_y right) = 0$.

Hence we are looking for a local solution of the last equation with $(v_x(0), v_y(0)) neq (0, 0)$.

So my question is:

Is there a simple way to prove local existence for a nonsingular solution of the last displayed equation?

I should make clear that I'm not looking for a proof of Newlander-Nirenberg; of this there are plenty. I am more interested in seeing what Newlander-Nirenberg becomes in terms of PDE in the simplest possible case, and then see that the PDE thus obtained is solvable. According to the answer of Andy, the equation which comes out is the Beltrami equation, so I will have a look at it. Still, I'm curious if any standard PDE technique can solve the equation I derived in the most stupid way above.