2020年6月30日火曜日

学習環境

続 解析入門 (原書第2版) (S.ラング(著)、松坂 和夫(翻訳)、片山 孝次(翻訳)、岩波書店)の第1章(ベクトル)、6(平面)の練習問題14の解答を求めてみる。



    1. {2a-b+c=03a+b+c=0
      5a+2c=0c=-52a
      b=2a-52a=-a2
      (2,-1,-5)

      2平面の 交わりの直線上の1点を求める。

      {2x-y+z=13x+y+z=2
      5x+2z=3z=-5x+32
      y=2x+-5x+32-1=-x+12
      (1,0,-1)

      よって、 求めるパラメーター方程式は

      (x,y,z)=(1,0,-1)+t(2,-1,-5)

    2. {2a+b+5c=03a-2b+c=0
      7a+11c=0c=-711a
      b=-2a+3511a=1311a
      (11,13,-7)
      {2x+y+5z=23x-2y+z=3
      7x+11z=7z=-7x+711
      y=-2x+35x-3511+2=13x-1311
      (x,y,z)=(1,0,0)+t(11,13,-7)

コード

#!/usr/bin/env python3
from unittest import TestCase, main
from sympy.plotting import plot3d, plot3d_parametric_line
from sympy.abc import x, y, t

print('14.')

fgs = [(-2 * x + y + 1, -3 * x - y + 2),
       ((2 - (2 * x + y)) / 5, 3 - (3 * x - 2 * y))]
ls = [(1 + 2 * t, -t, -1 - 5 * t),
      (1 + 11 * t, 13 * t, -7 * t)]
ts = [(-10, 10),
      (-2, 2)]
for i, ((f, g), l, (t1, t2)) in enumerate(zip(fgs, ls, ts), 12):
    p = plot3d(f, g,
               show=False)
    p.append(
        plot3d_parametric_line(
            *l,
            (t, t1, t2),
            legend=True,
            show=False,
        )[0]
    )
    p.save(f'sample14_{i}.png')

p.show()

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

% ./sample14.py
14.
%

0 コメント:

コメントを投稿