var('x')
solve(x^2 + 2 == 0, x)
solve(e^x == 5, x)
var('y')
solve([x+y==6, x-y==4], x, y)
var('b c')
solve([x^2 + b*x + c == 0],x)
f = x - (x/2 + x/3 + x/5) + (x/6 + x/10 + x/15) - x/30
type(f)
show(f)
f = (1/2)*log(2*x) + (1/2)*log(1/x)
show(f)
show(f.canonicalize_radical())