from sympy import *
from IPython.display import display
x,y,z, t,n = symbols('x y z t n')
eqc = Eq(x * (x-2)**2 * (x-5),0)
display(eqc)
print('Soluções:',solve(eqc))
print('Raízes :',roots(eqc))
Soluções: [0, 2, 5]
Raízes : {5: 1, 2: 2, 0: 1}
eqc = Eq(x**2 + y**2,oo)
display(eqc)
eqc = Eq(x**2 + y**2,10)
display(eqc)
display(solveset(eqc,x,domain=Reals))
eqs1 = [Eq(x*2 + 5*y,x + 10),Eq(x-y,y)]
for i in eqs1: display(i)
print('Soluções:',solve(eqs1))
print('\n\n')
eqs2 = [Eq(x*y - 1,x),Eq(x-2,y), Eq(z*10+2*y,E)]
for i in eqs2: display(i)
a = solve(eqs2)
print('Soluções:',solve(eqs2))
Soluções: {x: 20/7, y: 10/7}
Soluções: [{x: 3/2 - sqrt(13)/2, y: -sqrt(13)/2 - 1/2, z: 1/10 + E/10 + sqrt(13)/10}, {x: 3/2 + sqrt(13)/2, y: -1/2 + sqrt(13)/2, z: -sqrt(13)/10 + 1/10 + E/10}]
f,g,h = symbols('f g h', cls=Function)
display(Derivative(f(x),x))
display(f(x).diff(x))
eq1 = Eq(f(x).diff(x,2) - f(x).diff(x) + f(x),E**x + 30*x)
display(eq1)
print()
display(dsolve(eq1))
eq2 = Eq(f(x).diff(x,2), f(x) * g(y))
display(eq2)
print()
soluc = dsolve(eq2)
display('Solução:',soluc)
print('\n\n')
display('x = 2',
soluc.subs('x',2))
print('\n\n')
display('x = 5, C1 = log(5), C2 = 0',
soluc.subs([('x',5),('C1',log(5)),('C2',0)])
)
'Solução:'
'x = 2'
'x = 5, C1 = log(5), C2 = 0'