Processing math: 100%

2019年7月31日水曜日

学習環境

解析入門(上) (松坂和夫 数学入門シリーズ 4) (松坂 和夫(著)、岩波書店)の第5章(各種の初等関数)、5.1(対数関数・指数関数)、問題11の解答を求めてみる。



    1. ddxe-x2=-2xe-x2=(-1)(2x)e-x2d2dx2e-x2=(-1)(2e-x2+2x(-2xe-x2))=(-1)(2+2x(-2x))e-x2=(-1)2(22x2-2)e-x2d3dx3e-x2=(-1)2(23xe-x2+(22x2-2)(-2xe-x2))=(-1)2(23x+(22x2-2)(-2x))e-x2=(-1)3(-23x+24x3-22x)e-x2=(-1)3(24x3-(23+22)x)e-x2

      また、

      dndxne-x2=ddxdn-1dxn-1(e-x2)=ddx((-1)n-1Hn-1(x)e-x2]=(-1)n-1((ddxHn-1(x))e-x2+Hn-1(x)(-2xe-x2))=(-1)n(2xHn-1(x)+ddxHn-1(x))e-x2

      よって、帰納法により、

      Hn(x)

      は n 次の多項式で、 n が奇数ならば、奇関数、nが偶数ならば偶関数である。

コード

Python 3

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

print('11.')

x = symbols('x')
f = exp(-x ** 2)
fns = [Derivative(f, x, n).doit() for n in range(6)]

for n, fn in enumerate(fns):
    print(f'n = {n}')
    pprint(fn)
    print()

p = plot(*fns,
         (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('sample11.png')

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

C:\Users\...>py sample10.py
11.
n = 0
   2
 -x 
ℯ   

n = 1
        2
      -x 
-2⋅x⋅ℯ   

n = 2
                2
  ⎛   2    ⎞  -x 
2⋅⎝2⋅x  - 1⎠⋅ℯ   

n = 3
                  2
    ⎛       2⎞  -x 
4⋅x⋅⎝3 - 2⋅x ⎠⋅ℯ   

n = 4
                        2
  ⎛   4       2    ⎞  -x 
4⋅⎝4⋅x  - 12⋅x  + 3⎠⋅ℯ   

n = 5
                             2
    ⎛     4       2     ⎞  -x 
8⋅x⋅⎝- 4⋅x  + 20⋅x  - 15⎠⋅ℯ   


C:\Users\...>

0 コメント:

コメントを投稿