2019年6月11日火曜日

読書環境

ということで、PythonのSymPyを利用して微分してみた。

コード

Python 3

#!/usr/bin/env python3
from sympy import pprint, symbols, Derivative

常滑市 = symbols('常滑市')
a, b, c = symbols('a, b, c')
f = a * 常滑市 ** 2 + b * 常滑市 + c
fs = [Derivative(f, 常滑市, n) for n in range(4)]

for i, g in enumerate(fs):
    print(f'{i}階微分')
    for o in [g, g.doit()]:
        pprint(o)
        print()
    print()

入出力結果(Bash、cmd.exe(コマンドプロンプト)、Terminal、Jupyter(IPython))

C:\Users\...>py sample.py
0階微分
     2            
a⋅常滑市  + b⋅常滑市 + c

     2            
a⋅常滑市  + b⋅常滑市 + c


1階微分
 ∂  ⎛     2            ⎞
────⎝a⋅常滑市  + b⋅常滑市 + c⎠
∂常滑市                    

2⋅a⋅常滑市 + b


2階微分
   2                     
  ∂  ⎛     2            ⎞
─────⎝a⋅常滑市  + b⋅常滑市 + c⎠
    2                    
∂常滑市                     

2⋅a


3階微分
   3                     
  ∂  ⎛     2            ⎞
─────⎝a⋅常滑市  + b⋅常滑市 + c⎠
    3                    
∂常滑市                     

0


C:\Users\...>

pprint、日本語を含む式を綺麗に表示するのは得意ではないみたい。

0 コメント:

コメントを投稿

関連コンテンツ