2020年1月17日金曜日


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

学習環境

解析入門 原書第3版 (S.ラング(著)、松坂 和夫(翻訳)、片山 孝次(翻訳)、岩波書店)の第Ⅵ部(多変数の関数)、第17章(ベクトル)、4(ベクトルのノルム)の練習問題15の解答を求めてみる。



    • ⟨ sin 3 x , sin 3 x ⟩ = ∫ - π π sin 2 3 x dx = ∫ - π π 1 - cos 6 x 2 dx = 1 2 x - π π = π

    • ⟨ cos x , cos x ⟩ = ∫ - π π cos 2 x dx = ∫ - π π cos 2 x + 1 2 dx = π

コード

#!/usr/bin/env python3
from unittest import TestCase, main
from sympy import symbols, Integral, plot, sqrt, pi, sin, cos, plot

print('15.')

x = symbols('x')


def dot(f, g):
    return Integral(f * g, (x, -pi, pi)).doit()


def norm(f):
    return sqrt(dot(f, f))


class MyTestCase(TestCase):
    def test_sine(self):
        self.assertEqual(norm(sin(3 * x)), sqrt(pi))

    def test_cosine(self):
        self.assertEqual(norm(cos(x)), sqrt(pi))


p = plot(sin(3 * x), cos(x),
         (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(f'sample15.png')
if __name__ == '__main__':
    main()

入出力結果(Zsh、PowerShell、Terminal、Jupyter(IPython))

% ./sample15.py -v
15.
test_cosine (__main__.MyTestCase) ... ok
test_sine (__main__.MyTestCase) ... ok

----------------------------------------------------------------------
Ran 2 tests in 0.148s

OK
%
時刻: 12:00 ラベル: , , , ,

0 コメント:

コメントを投稿

コメントの投稿(Atom)