#範例3-3：多項式f(x)=0，求解x=? （solve，roots）
from sympy import *
x = symbols('x')
f1 = 'x**2-7*x+10'
#多項式f1(x)=0，求解x=?
#方法1：roots(f)求解
print('f1(x)=0，求解x=?',roots(f1))
#方法2：solve(f)求解
print('f1(x)=0，求解x=?',solve(f1))
print('f1(x)=0，求解x1=',solve(f1)[0])
print('f1(x)=0，求解x2=',solve(f1)[1])

#畫出f1的圖形
plot(f1,(x,-3,8))

#多項式f2(x)=x**2-9*x+10=0，求解x=?
f2 = 'x**2-9*x+10'
print('f2(x)=0，求解x=?',solve(f2))
print('f2(x)=0，求解x1=',solve(f2)[0])
print('f2(x)=0，求解x2=',solve(f2)[1])