2018年1月30日火曜日

学習環境

解析入門〈3〉(松坂 和夫(著)、岩波書店)の第14章(多変数の関数)、14.2(高次偏導関数、テイラーの定理)、問題2-(c).を取り組んでみる。


    1. f x = 1 1 + y x 2 - y x 2 = - y x 2 + y 2
      2 f x 2 = y · 2 x x 2 + y 2 2
      f y = 1 1 + y x 2 · 1 x = x x 2 + y 2
      2 f y 2 = - 2 x y x 2 + y 2 2

      よって、

      Δ f = 0

      ゆえに問題の関数は調和関数である。

コード(Emacs)

Python 3

#!/usr/bin/env python3
from sympy import pprint, symbols, atan, Limit, Derivative

x, y, t = symbols('x, y, t')
f = atan(y / x)
Ds = [Derivative(f, t, n) for n in range(1, 3)]

for D in Ds:
    for t0 in [x, y]:
        D = D.subs({t: t0})
        for s in [D, D.doit()]:
            pprint(s)
            print()
        print()
    print()

pprint(sum([Ds[1].subs({t: t0}).doit().factor() for t0 in [x, y]]))

入出力結果(Terminal, Jupyter(IPython))

$ ./sample2.py
∂ ⎛    ⎛y⎞⎞
──⎜atan⎜─⎟⎟
∂x⎝    ⎝x⎠⎠

    -y     
───────────
   ⎛     2⎞
 2 ⎜    y ⎟
x ⋅⎜1 + ──⎟
   ⎜     2⎟
   ⎝    x ⎠


∂ ⎛    ⎛y⎞⎞
──⎜atan⎜─⎟⎟
∂x⎝    ⎝x⎠⎠

    -y     
───────────
   ⎛     2⎞
 2 ⎜    y ⎟
x ⋅⎜1 + ──⎟
   ⎜     2⎟
   ⎝    x ⎠



  2         
 ∂ ⎛    ⎛y⎞⎞
───⎜atan⎜─⎟⎟
  2⎝    ⎝x⎠⎠
∂x          

    ⎛          2    ⎞
    ⎜         y     ⎟
2⋅y⋅⎜1 - ───────────⎟
    ⎜       ⎛     2⎞⎟
    ⎜     2 ⎜    y ⎟⎟
    ⎜    x ⋅⎜1 + ──⎟⎟
    ⎜       ⎜     2⎟⎟
    ⎝       ⎝    x ⎠⎠
─────────────────────
        ⎛     2⎞     
      3 ⎜    y ⎟     
     x ⋅⎜1 + ──⎟     
        ⎜     2⎟     
        ⎝    x ⎠     


  2         
 ∂ ⎛    ⎛y⎞⎞
───⎜atan⎜─⎟⎟
  2⎝    ⎝x⎠⎠
∂x          

    ⎛          2    ⎞
    ⎜         y     ⎟
2⋅y⋅⎜1 - ───────────⎟
    ⎜       ⎛     2⎞⎟
    ⎜     2 ⎜    y ⎟⎟
    ⎜    x ⋅⎜1 + ──⎟⎟
    ⎜       ⎜     2⎟⎟
    ⎝       ⎝    x ⎠⎠
─────────────────────
        ⎛     2⎞     
      3 ⎜    y ⎟     
     x ⋅⎜1 + ──⎟     
        ⎜     2⎟     
        ⎝    x ⎠     



0
$

0 コメント:

コメントを投稿