2020年5月10日日曜日

学習環境

解析入門 原書第3版 (S.ラング(著)、松坂 和夫(翻訳)、片山 孝次(翻訳)、岩波書店)の第Ⅵ部(多変数の関数)、第20章(合成微分律と勾配ベクトル)、2(接平面)の練習問題5の解答を求めてみる。



    1. gradx2+y2-z=2x,2y,-12·3,2·4,-1=6,8,-1

      よって、求める接平面の方程式は、

      6x+8y-z=6·3+8·4-256x+8y-z=18+32-256x+8y-z=25

    2. gradxx2+y212-z=x2+y212-x·12·1x2+y212·2xx2+y2,-xyx2+y212x2+y2,-1

      点 P における勾配ベクトル。

      5-9525,125·25,-1=165·25,125·25,-1

      接平面の方程式。

      165·25x+125·25y-z=16·35·25+12·-45·25-3516x+12y-125z=48-48-7516x+12y-125z=-75

    3. gradsinxy-z=ycosxy,xcosxy,-1πcosπ,cosπ,-1=-π,-1,-1-πx-y-z=-π-ππx+y+z=2π

コード

#!/usr/bin/env python3
from sympy import symbols, sin, sqrt, pi
from sympy.plotting import plot3d

print('5.')

x, y = symbols('x, y')

fs = [x ** 2 + y ** 2,
      x / sqrt(x ** 2 + y ** 2),
      sin(x * y)]
gs = [6 * x + 8 * y - 25,
      (16 * x + 12 * y + 75) / 125,
      (2 * pi - pi * x - y) / 2]

for i, (f, g) in enumerate(zip(fs, gs)):
    c = chr(ord("a") + i)
    print(f'({c})')
    p = plot3d(f, g,
               (x, -5, 5),
               (y, -5, 5),
               show=False)
    p.save(f'sample5_{c}.png')

p.show()

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

% ./sample5.py
5.
(a)
(b)
(c)
%

0 コメント:

コメントを投稿