Loading [MathJax]/jax/output/HTML-CSS/jax.js

2018年12月6日木曜日

学習環境

数学読本〈1〉数・式の計算/方程式/不等式 (松坂 和夫(著)、岩波書店)の第3章(数学の威力を発揮する - 方程式)、3.3(高次方程式)、高次方程式の解法の問25.を取り組んでみる。


  1. x3-a3=0(x-a)(x2+ax+a2)=0x2+ax+a2=0x=-a±a2-4a22=-a+-3ai2a=-1±3i2a

    よって、問題の方程式の解は、

    x=a,-1±3i2a=a,aω,-1-3i2a

    また、

    ω2=(-1+3i2)2=1-3-23i22=-2-23i22=-1-3i2

    よって、 問題の方程式の解は、

    x=a,ωa,ω2a

    である。

コード(Emacs)

Python 3

#!/usr/bin/env python3
from sympy import pprint, symbols, solve, sqrt, I

print('25.')

x = symbols('x')
a = symbols('a', real=True)
eq = x ** 3 - a ** 3
xs = solve(eq, x)
pprint(xs)
w = (-1 + sqrt(3) * I) / 2

for l, r in zip(xs[1:], [w ** 2 * a, w * a]):
    for o in [l, r, l.expand() == r.expand()]:
        pprint(o)
        print()
    print()

入出力結果(Terminal, cmd(コマンドプロンプト), Jupyter(IPython))

$ ./sample25.py
25.
⎡     ⎛  1   √3⋅ⅈ⎞    ⎛  1   √3⋅ⅈ⎞⎤
⎢a, a⋅⎜- ─ - ────⎟, a⋅⎜- ─ + ────⎟⎥
⎣     ⎝  2    2  ⎠    ⎝  2    2  ⎠⎦
  ⎛  1   √3⋅ⅈ⎞
a⋅⎜- ─ - ────⎟
  ⎝  2    2  ⎠

              2
  ⎛  1   √3⋅ⅈ⎞ 
a⋅⎜- ─ + ────⎟ 
  ⎝  2    2  ⎠ 

True


  ⎛  1   √3⋅ⅈ⎞
a⋅⎜- ─ + ────⎟
  ⎝  2    2  ⎠

  ⎛  1   √3⋅ⅈ⎞
a⋅⎜- ─ + ────⎟
  ⎝  2    2  ⎠

True


$

0 コメント:

コメントを投稿