数学 - Python - 解析学 - 積分 - 積分の応用 - 曲線の長さ(三角関数(余弦)、半角、加法定理、極座標表示)

解析入門 原書第3版
原著
楽天ブックス(Kobo) Yahoo!

学習環境

• Surface Go、タイプ カバー、ペン(端末)
• Windows 10 Pro (OS)
• Nebo(Windows アプリ)
• iPad Pro + Apple Pencil
• MyScript Nebo - MyScript(iPad アプリ(iOS))
• 参考書籍
  • Pythonからはじめる数学入門 (楽天ブックス(Kobo) Yahoo!)
  • JavaScriptリファレンス 第6版 (楽天ブックス(Kobo) Yahoo!)
  • インタラクティブ・データビジュアライゼーション (楽天ブックス(Kobo) Yahoo!)

解析入門 原書第3版 (S.ラング(著)、松坂 和夫(翻訳)、片山 孝次(翻訳)、岩波書店)の第3部(積分)、第13章(積分の応用)、補充問題、曲線の長さの練習問題19の解答を求めてみる。

19. 
    
    $\begin{array}{}\underset{0}{\overset{\mathrm{\pi \; <!--\; greek\; small\; letter\; pi\; -->}}{\int \; <!--\; integral\; -->}}\sqrt{{\left(1-\cos \mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}\right)}^2+{\left(\frac{d}{d\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}}\left(1-\cos \mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}\right)\right)}^2}d\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}\\ =\underset{0}{\overset{\mathrm{\pi \; <!--\; greek\; small\; letter\; pi\; -->}}{\int \; <!--\; integral\; -->}}\sqrt{1-2\cos \mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}+{\cos }^2\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}+{\sin }^2\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}}d\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}\\ =\underset{0}{\overset{\mathrm{\pi \; <!--\; greek\; small\; letter\; pi\; -->}}{\int \; <!--\; integral\; -->}}\sqrt{2-2\cos \mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}}d\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}\\ =\sqrt{2}\underset{0}{\overset{\mathrm{\pi \; <!--\; greek\; small\; letter\; pi\; -->}}{\int \; <!--\; integral\; -->}}\sqrt{1-\cos \left(\frac{\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}}{2}+\frac{\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}}{2}\right)}d\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}\\ =\sqrt{2}\underset{0}{\overset{\mathrm{\pi \; <!--\; greek\; small\; letter\; pi\; -->}}{\int \; <!--\; integral\; -->}}\sqrt{1-\left({\cos }^2\left(\frac{\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}}{2}\right)-{\sin }^2\left(\frac{\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}}{2}\right)\right)}d\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}\\ =\sqrt{2}\underset{0}{\overset{\mathrm{\pi \; <!--\; greek\; small\; letter\; pi\; -->}}{\int \; <!--\; integral\; -->}}\sqrt{2{\sin }^2\left(\frac{\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}}{2}\right)}d\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}\\ =2\underset{0}{\overset{\mathrm{\pi \; <!--\; greek\; small\; letter\; pi\; -->}}{\int \; <!--\; integral\; -->}}\sin \left(\frac{\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}}{2}\right)d\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}\\ =2{\left[-2\cos \left(\frac{\mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}}{2}\right)\right]}_0^{\mathrm{\pi \; <!--\; greek\; small\; letter\; pi\; -->}}\\ =-4\left(\cos \frac{\mathrm{\pi \; <!--\; greek\; small\; letter\; pi\; -->}}{2}-\cos \mathrm{\theta \; <!--\; greek\; small\; letter\; theta\; -->}\right)\\ =-4\left(0-1\right)\\ =4\end{array}$

コード

Python 3

#!/usr/bin/env python3
from sympy import pprint, symbols, Integral, Derivative, sqrt, cos, sin, pi
from sympy.plotting import plot_parametric

theta = symbols('θ')
r = 1 - cos(theta)
x = r * cos(theta)
y = r * sin(theta)

I = Integral(sqrt(r ** 2 + Derivative(r, theta, 1) ** 2),
             (theta, 0, pi))

for o in [I, I.doit()]:
    pprint(o.simplify())
    print()

I = 2 * Integral(sin(theta / 2), (theta, 0, pi))
for o in [I, I.doit()]:
    pprint(o.simplify())
    print()

p = plot_parametric((x, y, (theta, 0, pi)),
                    (x, y, (theta, pi, 2 * pi)),
                    show=False)

colors = ['red', 'green']
for i, color in enumerate(colors):
    p[i].line_color = color
p.save('sample19.png')

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

$ python3 sample19.py
π                                             
⌠                                             
⎮      ____________________________________   
⎮     ╱                                  2    
⎮    ╱              2   ⎛d              ⎞     
⎮   ╱   (cos(θ) - 1)  + ⎜──(-cos(θ) + 1)⎟   dθ
⎮ ╲╱                    ⎝dθ             ⎠     
⌡                                             
0                                             

π                     
⌠                     
⎮   _______________   
⎮ ╲╱ -2⋅cos(θ) + 2  dθ
⌡                     
0                     

  π          
  ⌠          
  ⎮    ⎛θ⎞   
2⋅⎮ sin⎜─⎟ dθ
  ⎮    ⎝2⎠   
  ⌡          
  0          

4

$