Sageで方程式を解く

Sageで方程式を解いてみる。

x,a,b,c=var('x a b c')
solve(x^2+1==0,x)

f=a*x^2+b*x+c==0
solve(f,x)

Q=solve(f,x)
for q in Q:
    q.substitute({a:1,b:-5,c:6})

Latexの書式を出すときには

latex(_)

などとする。