2019年2月5日火曜日

学習環境

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


  1. 0 2 π π r 2 2 π d θ = 1 2 0 2 π 2 + sin 2 θ 2 d θ = 1 2 0 2 π 4 + 4 sin 2 θ + sin 2 2 θ d θ = 1 2 4 θ - 2 cos 2 θ 0 2 π + 1 2 0 2 π sin 2 2 θ d θ = 1 2 · 8 π + 1 2 0 2 π sin 2 2 θ d θ t = 2 θ dt d θ = 2 sin 2 t · 1 2 dt = 1 2 - sin t cos t + 1 2 1 dt = 1 2 - sin t cos t + 1 2 t θ = 0 , t = 0 θ = 2 π , t = 4 π 0 4 π sin 2 t · 1 2 dt = 1 2 - sin t cos t 0 4 π + 1 2 t 0 4 π = π 4 π + 1 2 · π = 9 2 π

コード

Python 3

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

theta = symbols('θ')
r = 2 + sin(2 * theta)
x = r * cos(theta)
y = r * sin(theta)

I = Integral(r ** 2 / 2, (theta, 0, 2 * pi))

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

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


colors = ['red', 'green', 'blue', 'brown']
for i, s in enumerate(p):
    s.line_color = colors[i]
p.save('sample4.png')

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

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

π

$

落花生っぽい(?)形になった。

0 コメント:

コメントを投稿