import numpy as np
import sympy as sym
from IPython.display import display,Math
display(Math('x^2 + 1 = 0'))
display(Math('x^2 = -1'))
display(Math('x = \\pm \\sqrt 1'))
print(1j)
print(np.sqrt(-1))
print(np.sqrt(-1,dtype='complex'))
print((-1)**(1/2))
print(sym.I)
real_part = 4
imag_part = -5
cn1 = np.complex(real_part,imag_part)
cn2 = real_part + 1j*imag_part
print(type(cn1))
print(cn1)
print(type(cn2))
print(cn2)
z1 = np.complex(4,5)
z2 = np.complex(3,2)
z1+z2
z1-z2
np.real(z1)
np.imag(z1)
w = np.complex(2,4)
z = np.complex(5,6)
realsum = np.real(w) + np.real(z)
imagsum = np.imag(w) + np.imag(z)
print(realsum, imagsum)
wPlusZ_separate = np.complex(realsum,imagsum)
wPlusZ_direct = w+z
print(wPlusZ_separate)
print(wPlusZ_direct)