Sageで方程式

Sageで200万が5年後に202万になる利息を計算する。

x=var(x)
q=solve(200*(1+x)**5==202,x)
for n in range(len(q)):
    print N(q[n].rhs())

結果

-0.690367429042489 + 0.952951066209182*I
-1.81062859479078 + 0.588956148532726*I
-1.81062859479078 - 0.588956148532726*I
-0.690367429042489 - 0.952951066209181*I
0.00199204766653335

したがって0.00199204766653335=約0.199%

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(_)

などとする。

Sageで数学 素数

Sageで素数を遊んでみる。

sage: A=Primes()
sage: A.first()
2
sage: A.next(2)
3
sage: A.next(3)
5
sage: A.next(5)
7
sage: A.next(4)
5
sage: A.unrank(0)
2
sage: A.unrank(4)　#Returns the 4-th prime number.
11

sage: A.unrank(100)

547
sage: 17 in A
True
sage: A.cardinality() # Counts the elements of the enumerated set.
+Infinity




sage: prime_pi(548) # Return the number of primes <=549
101
sage: plot(prime_pi,1,600)