#範例3-3：方程式f(x)=0，求解x=? （solve，roots）
from sympy import *
x = symbols('x')
f1 = Eq(x**2-7*x+10, 0)
#方程式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])

#畫出f的圖形
#若要畫出圖形，必須是函數f(x)，才能畫圖
f2 = Function('f2')(x)
f2 = x**2-7*x+10
plot(f2,(x,-3,8))

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