## 2019年6月23日日曜日

### 数学 - Python - 解析学 - 級数 - テイラーの公式 - 三角関数(正弦)、累乗(べき乗、平方)、差、商、極限

1. $\begin{array}{l}\frac{d}{\mathrm{dx}}\left(\mathrm{sin}x-x\right)\\ =\mathrm{cos}x-1\\ \frac{{d}^{2}}{{\mathrm{dx}}^{2}}\left(\mathrm{sin}x-x\right)\\ =-\mathrm{sin}x\\ \frac{{d}^{3}}{{\mathrm{dx}}^{3}}\left(\mathrm{sin}x-x\right)\\ =-\mathrm{cos}x\\ \frac{{d}^{4}}{{\mathrm{dx}}^{4}}\left(\mathrm{sin}x-x\right)\\ =\mathrm{sin}x\\ \frac{{d}^{5}}{{\mathrm{dx}}^{5}}\left(\mathrm{sin}x-x\right)\\ =\mathrm{cos}x\\ \frac{{d}^{6}}{{\mathrm{dx}}^{6}}\left(\mathrm{sin}x-x\right)\\ =-\mathrm{sin}x\\ \mathrm{sin}x-x=-\frac{1}{3!}{x}^{3}+\frac{1}{5!}{x}^{5}-\dots \\ \underset{x\to 0}{\mathrm{lim}}\frac{\mathrm{sin}x-x}{{x}^{2}}=0\end{array}$

コード

Python 3

#!/usr/bin/env python3
from sympy import pprint, symbols, plot, Limit, sin

print('26.')

x = symbols('x')
f = (sin(x) - x) / x ** 2

for dir in ['+', '-']:
l = Limit(f, x, 0, dir=dir)
for o in [l, l.doit()]:
pprint(o)
print()

p = plot(sin(x), -x, sin(x) - x, x ** 2, f,
(x, -5, 5),
ylim=(-5, 5),
legend=True,
show=False)
colors = ['red', 'green', 'blue', 'brown', 'orange',
'purple', 'pink', 'gray', 'skyblue', 'yellow']

for o, color in zip(p, colors):
o.line_color = color

p.show()
p.save('sample26.png')


C:\Users\...>py sample26.py
26.
⎛-x + sin(x)⎞
lim ⎜───────────⎟
x─→0⁺⎜      2    ⎟
⎝     x     ⎠

0

⎛-x + sin(x)⎞
lim ⎜───────────⎟
x─→0⁻⎜      2    ⎟
⎝     x     ⎠

0

C:\Users\...>