Ocean Science Hack: Sageで方程式を解く

2010年9月17日金曜日

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)

$\left[x = \frac{-b + \sqrt{-4 \, a c + b^{2}}}{2 \, a}, x = \frac{-b - \sqrt{-4 \, a c + b^{2}}}{2 \, a}\right]$

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

$x=2$
$x=3$

Latexの書式を出すときには

latex(_)

などとする。