Loading [MathJax]/jax/output/HTML-CSS/jax.js

2018年1月30日火曜日

学習環境

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


    1. fx=11+(yx)2-yx2=-yx2+y2
      2fx2=y·2x(x2+y2)2
      fy=11+(yx)2·1x=xx2+y2
      2fy2=-2xy(x2+y2)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 コメント:

コメントを投稿