Sabemos que el valor principal de la función argumento de un número complejo $z$ se define como
$$\text{Arg} (z)=\begin{cases} \arctan\left(\frac{\text{Im}(z)}{\text{Re}(z)}\right)&,\text{Re}(z)>0\\\\ \arctan\left(\frac{\text{Im}(z)}{\text{Re}(z)}\right)+\pi&,\text{Re}(z)<0,\,\text{Im}(z)\ge 0\\\\ \arctan\left(\frac{\text{Im}(z)}{\text{Re}(z)}\right)-\pi&,\text{Re}(z)<0,\,\text{Im}(z)< 0\\\\ \pi/2&,\text{Re}(z)=0,\,\text{Im}(z)> 0\\\\ -\pi/2&,\text{Re}(z)=0,\,\text{Im}(z)< 0\\\\ \end{cases}$$
El argumento de $z$ es
$$\arg (z)=\text{Arg}(z)+2\ell \pi$$
para valores enteros de $\ell$ .
Ahora, para $z=x+iy$ con valores fijos de $x$ y $y$ y $n>-x$ tenemos
$$\begin{align} \arg\left(\left(1+\frac zn\right)^n\right)&=\arg\left(\left(1+\frac xn+i\frac yn\right)^n\right) \tag 1\\\\ &=n\arctan\left(\frac{y/n}{1+x/n}\right)+2\ell \pi \tag 2\\\\ &\to y+2\ell \pi \,\,\text{as}\,\,n\to \infty \tag 3 \end{align}$$
Por lo tanto, ya que $\cos (y+2\ell \pi)=\cos y$ y $\sin (y+2\ell \pi)=\sin y$ tenemos
$$\lim_{n\to \infty}\left(1+\frac zn\right)^n=e^x\left(\cos y+i\sin y\right)$$
¡como se iba a mostrar!
NOTA $1$ :
Al pasar de $(1)$ a $(2)$ , hicimos uso de forma inductiva de la relación
$$\arg(z^2)=2\arg(z)+2\ell \pi \tag 4$$
Podemos verificar $(4)$ como sigue. En primer lugar, dejemos que $z=x+iy$ . Entonces, para $x^2-y^2>0$ y $x>0$ tenemos
$$\begin{align} \arg(z^2)&=\text{Arg}(z^2)+2\ell \pi\\\\ &=\arctan\left(\frac{2xy}{x^2-y^2}\right)+2\pi \ell\\\\ &=\arctan\left(\frac{2(y/x)}{1-(y/x)^2}\right)+2\pi \ell\\\\ &=2\arctan(y/x)+2\pi \ell\,\,\dots \,\text{where we used the double angle formula for the tangent}\\\\ &=2\arg(z)+2\ell \pi \tag 5 \end{align}$$
Ahora, dejemos que $z=\left(1+\frac xn\right)+i\left(\frac yn\right)$ en $(5)$ . Entonces, para $n$ suficientemente grande, $\left(1+\frac xn\right)^2-\left(\frac yn\right)^2>0$ y $1+\frac xn >0$ encontramos
$$\arg\left(z^2\right)=2\arg(z)+2\pi \ell$$
El procedimiento inductivo revela que
$$\arg(z^n)=n\arg (z)+2\ell \pi$$
como se iba a demostrar.
NOTA $2$ :
Al pasar de $(2)$ a $(3)$ , hicimos uso de ESTA RESPUESTA en el que establecí las desigualdades para la función arctangente
$$\frac{|x|}{\sqrt{1+x^2}} \le|\arctan(x)|\le |x|$$
Entonces, tenemos
$$\left|\frac{\frac{y}{1+x/n}}{\sqrt{1+\left(\frac{y/n}{1+x/n}\right)^2}} \right|\le \left| n\arctan\left(\frac{y/n}{1+x/n}\right)\right| \le \left|\frac{y}{1+x/n}\right|$$
con lo que aplicando el teorema de la compresión se obtiene el límite
$$\lim_{n\to \infty}n\arctan\left(\frac{y/n}{1+x/n}\right)=y$$
¡como se iba a mostrar!