print(log(pi).n())
print(e^(1/3).n())
print(tan(arccot(sqrt(3)/3)).n())
print(i^(i).n())

W = 18*x^4 + 123*x^3 + 89*x^2 - 8*x - 12
print(W.factor())
pretty_print(W(x=(1/2)))
W2 = W*(x+1)
print(W2.expand())
W3= W(x=(-sqrt(3)))
print(round(W3, 2)) #przybliza cale wyrazenie do liczb calkowitych 

a, b, c = -1, 3, 5
f = (-1)*x^2 + 3*x + 5

p = -b/(2*a)
q = f(x=p)

rys = plot(f, (x, -4, 4), color="#990666", thickness=2, ymin=-10, ymax=10, gridlines=True, title='Wykres funkcji $f(x)$')
krp = point((p,q), color="#444666", size=40)
tekst = text(f'W({p},{q})', (1.5, 8))

rys + krp + tekst

fun = plot((2*x+1)/(x-1), (x, -6, 8), ymin=-10, ymax=12, detect_poles=True, axes_labels=['$x$', '$y$'], legend_label='$y=(2*x+1)/(x-1)$')
as1 = line([(1, -12), (1, 14)], linestyle="--", color="#990666", legend_label='$as1: x=1$')
as2 = line([(-8, 2), (10, 2)], linestyle="--", color="#990666", legend_label='$as2: y=2$')
fun.set_legend_options(loc='lower right') # z racji tego ze w ostatniej lini sumuje wszystkie rysunki wystarczy ze zmienie polozenie legendy tylko dla jednego

fun + as1 + as2

solve(x^4 - 4*x^3 - 3*x^2 - 2*x - 1 == 0, x)
x in RR
find_root(x^4 - 4*x^3 - 3*x^2 - 2*x - 1 == 0, -1, 1)
#cdn