I'm interested in continuous conformal maps $\mathbb{R}^2 \rightarrow \mathbb{R}^2$ that are ONLY angle preserving (not necessarily orientation preserving). I would like to write down a generic system of PDEs such that ALL conformal maps must obey this PDE (or maybe some other suitable functional equation).
It's a well known fact that all holomorphic functions are conformal and a holomorphic function can basically be characterized as a map $D: \mathbb{R}^2 \rightarrow \mathbb{R}^2$ with two components $u(x,y),v(x,y)$ such that the cauchy riemann equations hold.
$$ \frac{\partial u}{\partial x} = \frac{ \partial v}{\partial y} \\ \frac{\partial u}{\partial y} = - \frac{ \partial v}{\partial x} $$
But not ALL conformal maps are holomorphic... namely if we move that minus sign to the upper equation we would get an "anti holomorphic" cauchy riemann equation. (See here for an explanation).
So there are at least two different systems of two PDEs whose general solutions result in conformal maps.
Now I haven't seen any result to suggest that harmonic functions are conformal but they certainly seem related and moreover both holomorphic and antiholomorphic functions are harmonic.
So that leads us to the question, can we some characterize ALL conformal maps by a single system of PDEs?
An attempt at an approach:
We try to do this the ugly way... Suppose we have two parametrized curves $\gamma_1 = (\gamma_{1,x}(t), \gamma_{1,y}(t))$ and $\gamma_2 = (\gamma_{2,x}(t), \gamma_{2,y}(t))$
Now suppose these intersect at $t = t_0$ we can then discuss the angle between the curves as
$$ \theta_{\text{original}} = \left| \tan^{-1} \left( \frac{\frac{d\gamma_{2,y}}{dt}}{ \frac{d\gamma_{2,x}}{dt} } \right) - \tan^{-1} \left( \frac{\frac{d\gamma_{1,y}}{dt}}{ \frac{d\gamma_{1,x}}{dt} } \right) \right| @ t = t_0 $$
Now we consider some arbitrary map on $\mathbb{R^2} \rightarrow \mathbb{R^2}$ consisting of two components $u(x,y), v(x,y)$ where the $u$ is our x coordinate and $v$ is our y coordinate.
So we talk about the angle between the two curves after applying this map and set the two things equal to each other. So that angle is given as:
$$ \theta_{\text{new}} = \left| \tan^{-1} \left( \frac{\frac{d v\left( \gamma_{2,x}, \gamma_{2,y} \right) }{dt}}{ \frac{d u\left( \gamma_{2,x}, \gamma_{2,y} \right)}{dt} } \right) - \tan^{-1} \left( \frac{\frac{d v\left( \gamma_{1,x}, \gamma_{1,y} \right)}{dt}}{ \frac{d u\left( \gamma_{1,x}, \gamma_{1,y} \right) }{dt} } \right) \right| @ t = t_0 $$
which we can expand via the chain rule to find
$$ \left| \tan^{-1} \left( \frac{\frac{\partial v }{\partial x } \frac{d \gamma_{2,x}}{dt} + \frac{\partial v}{\partial y} \frac{d \gamma_{2,y}} {dt} }{ \frac{\partial u }{\partial x } \frac{d \gamma_{2,x}}{dt} + \frac{\partial u}{\partial y} \frac{ d\gamma_{2,y}}{dt}} \right) - \tan^{-1} \left( \frac{\frac{\partial v }{\partial x } \frac{d \gamma_{1,x}}{dt} + \frac{\partial v}{\partial y} \frac{d \gamma_{1,y}} {dt} }{ \frac{\partial u }{\partial x } \frac{d \gamma_{1,x}}{dt} + \frac{\partial u}{\partial y} \frac{ d\gamma_{1,y}}{dt}} \right) \right|. $$
So a characterization of all conformal curves maps is the following:
$$ \left| \tan^{-1} \left( \frac{\frac{\partial v }{\partial x } \frac{d \gamma_{2,x}}{dt} + \frac{\partial v}{\partial y} \frac{d \gamma_{2,y}} {dt} }{ \frac{\partial u }{\partial x } \frac{d \gamma_{2,x}}{dt} + \frac{\partial u}{\partial y} \frac{ d\gamma_{2,y}}{dt}} \right) - \tan^{-1} \left( \frac{\frac{\partial v }{\partial x } \frac{d \gamma_{1,x}}{dt} + \frac{\partial v}{\partial y} \frac{d \gamma_{1,y}} {dt} }{ \frac{\partial u }{\partial x } \frac{d \gamma_{1,x}}{dt} + \frac{\partial u}{\partial y} \frac{ d\gamma_{1,y}}{dt}} \right) \right| = \left| \tan^{-1} \left( \frac{\frac{d\gamma_{2,y}}{dt}}{ \frac{d\gamma_{2,x}}{dt} } \right) - \tan^{-1} \left( \frac{\frac{d\gamma_{1,y}}{dt}}{ \frac{d\gamma_{1,x}}{dt} } \right) \right| $$
Of course this definition sucks because of the dependence on the arbitrary $\gamma_1$ and $\gamma_2$ curves, which we would like to remove somehow. It's not even obvious if $u,v$ is holomorphic (meaning obyeing cauchy riemann) that this equation holds true (perhaps some trig identity is required there).