In [1]:
from sympy import *
from IPython.display import display,Math
r1,r2,t1,t2 = symbols("rho_1 rho_2 theta_1 theta_2",constant=true,real=true)
z1 = r1*exp(I*t1)
z2 = r2*exp(I*t2)
In [2]:
produtoZ = z1*z2
produtoZ = produtoZ.simplify()
display(Math("z_1 z_2 = "+latex(produtoZ)))
In [3]:
divisaoZ = z1/z2
divisaoZ = divisaoZ.simplify()
display(Math("\\frac{z_1}{z_2} = "+ latex(divisaoZ)))
In [4]:
n = symbols("n",real=true,constant=true)
expZ = z1**n
expZ = expZ.simplify()
display(Math("z_1^n = "+ latex(expZ)))
Avaliando a raiz:
$$w = z^{\frac{1}{n}} = \sqrt[n]{z}$$Supondo $w$ tal que
$$w = \rho (\cos(\phi) + \sin(\phi))$$A partir do resultado B3 elevando $w^n$ temos
$$w^n = \rho^n (\cos(n \phi) + \sin(n \phi))$$Porem $w^n =z$ e portanto por igualdade de números complexos temos que $\rho^n = p$ e $n \theta = \phi$. Isolando $\rho$ e $\phi$:
$$ \left\{\begin{array}{l} \rho = r^{\frac{1}{n}}\\ \phi = \frac{\theta + 2k\pi}{n} \end{array}\right. \ ,\ n \neq 0 $$E portanto temos para todas as n raizes enésimas de $z$
$$w = \sqrt[n]{r} \left(\cos\left(\frac{\theta + 2k\pi}{n}\right) + \sin\left(\frac{\theta + 2k\pi}{n}\right)\right)$$com $0\leq k \leq n-1$.