I have learned that one can derive the Cauchy-Riemann equations by using the "typical" definition of differentiation (i.e. $f'(z)=\lim_{h \rightarrow 0} \frac{f(z+h)-f(z)}{h}$) by letting $h=h_1+i0$ and calculating the limit and then letting $h=0+ih_2$ and calculating the limit.
I know that besides the (what I called typical) definition of differentiability there also is a equivalent formulation for it:
A function $f:(a,b) \rightarrow \mathbb{R}$ is differentiable in the point $p \in (a,b)$ if and only if $f(p+h)-f(p)=f'(p)h+r(h)$ with $\lim_{h \rightarrow 0} \frac{r(h)}{h}=0$.
I would like to use this definition of differentiability to derive the Cauchy-Riemann equations.
Here is my attempt:
Let $f$ be a complex-valued function defined on some open set $U \subseteq \mathbb{C}$.
If $f$ complex differentiable at some point $z \in U$, then we can write
$f(z+h)-f(z)=f'(z)h+r(h)$ with $\lim_{h \rightarrow 0} \frac{r(h)}{h}=0$. (with $h=h_1+ih_2,r(h):\mathbb{C} \rightarrow \mathbb{C}$)
Let $f(x+iy)=u(x,y)+iv(x,y)$, $f'(z)=a+ib$ and $r(h)=r_1(h_1,h_2)+ir_2(h_1,h_2)$.
Now let $g(x,y):=(u(x,y),v(x,y))$, where we view g as a function $\mathbb{R}^2 \rightarrow \mathbb{R}^2$.
For a function $g:\mathbb{R}^2 \rightarrow \mathbb{R}^2$ to be differentiable in the point $(x,y)$ is equivlanet that $g(x+h_1,y+h_2)-g(x,y)=A \cdot (h_1,h_2)+(r_1(h_1,h_2),r_2(h_1,h_2))$ with $\lim_{h \rightarrow 0} \frac{(r_1(h_1,h_2),r_2(h_1,h_2))}{\|(h_1,h_2)\|}=0$, where $A$ is the Jacobi Matrix of $g$.
Then I did the following calculations:
$f'(z)h=(a+ib)(h_1+ih_2)=h_1 a-h_2 b+i(h_1b+h_2a)$
$A \cdot (h_1,h_2)=(h_1 u_x+h_2 u_y,h_1 v_x+h_2v_y)$.
Doing now the comparison we get the CR equations $u_x=v_y$ and $u_y=-v_x$.
Now I am stuck on showing that if $\lim_{h \rightarrow 0} \frac{r(h)}{h}=0$ then so does $\lim_{h \rightarrow 0} \frac{(r_1(h_1,h_2),r_2(h_1,h_2))}{\|(h_1,h_2)\|}=0$ (in $\mathbb{R^2}$).
How can I argue here? (And also, is my approach this far even correct?)