2017年5月24日水曜日

開発環境

ふと、SymPy で三角関数の値が正確に求まるのか気になったから試してみた。(結果を求めるので過程は関係ない。)

コード(Emacs)

#!/usr/bin/env python3
# -*- coding: utf-8 -*-

from sympy import Symbol, solve, pprint, factor, cos, pi

z = Symbol('z')
expr = z ** 4 + z ** 3 + z ** 2 + z + 1

for s in solve(expr, z):
    pprint(s)

for s in solve(expr, z):
    pprint(factor(s))

a = cos(2 * pi / 5)
pprint(a)
pprint(factor(a))

入出力結果(Terminal)

$ ./sample3.py 
                 ________
  1   √5        ╱ √5   5 
- ─ + ── - ⅈ⋅  ╱  ── + ─ 
  4   4      ╲╱   8    8 
                 ________
  1   √5        ╱ √5   5 
- ─ + ── + ⅈ⋅  ╱  ── + ─ 
  4   4      ╲╱   8    8 
                 __________
  √5   1        ╱   √5   5 
- ── - ─ - ⅈ⋅  ╱  - ── + ─ 
  4    4     ╲╱     8    8 
                 __________
  √5   1        ╱   √5   5 
- ── - ─ + ⅈ⋅  ╱  - ── + ─ 
  4    4     ╲╱     8    8 
 ⎛                 ________⎞ 
-⎝-√5 + 1 + √2⋅ⅈ⋅╲╱ √5 + 5 ⎠ 
─────────────────────────────
              4              
                 ________
-1 + √5 + √2⋅ⅈ⋅╲╱ √5 + 5 
─────────────────────────
            4            
 ⎛                _________⎞ 
-⎝1 + √5 + √2⋅ⅈ⋅╲╱ -√5 + 5 ⎠ 
─────────────────────────────
              4              
                 _________
-√5 - 1 + √2⋅ⅈ⋅╲╱ -√5 + 5 
──────────────────────────
            4             
  1   √5
- ─ + ──
  4   4 
-1 + √5
───────
   4
$   

0 コメント:

コメントを投稿

Comments on Google+: