2020年3月22日日曜日

学習環境

解析入門(上) (松坂和夫 数学入門シリーズ 4) (松坂 和夫(著)、岩波書店)の第10章(n次元空間)、10.1(ユークリッド空間)、問題4の解答を求めてみる。



    1. x-a=x-a+b2-a-a+b2=x-a+b2-a-b2x-b=x-a+b2-b-a+b2=x-a+b2+a-b2

      ここで、

      x-a+b2=y=y1,,yaa-b2=c=c1,,cn

      とおくと、

      x-a=y-cx-b=y+c

      よって、問題の方程式は、

      y-c=y+c=r

      この解を考えてみる。

      仮定の

      2r>d

      より

      2r>a-br>a-b2=c

      また、

      y-c=y+cy-c2=y+c2y-c·y-c=y+c·y+c-2y·c=2y·cy·c=0c1y1++cnyn=0y-c=ry-c2+y+c2=2r2y2+c2=r2y2=r2-c2y12++yn2=r2-c2>0

      問題の仮定より

      n3

      なので、 y、 すなわち x は無限に存在する。

      (証明終)


    2. 2r=d

      ならば、

      r=d2=cr2-c2=0

      よって、 方程式の解は

      y=0x=a+b2

      で、 ただ1つである。

      (証明終)


    3. 2r<dr<cr2-c2<0

      の場 合、

      y·y=r2-c2

      を満たす y、 すなわち x は存在しない。

      (証明終)

コード

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

print('4.')

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

for i, n in enumerate([Rational(1, 2), 1, 2]):
    p = plot3d(*[s2 * sqrt(-((x + s1 * n) ** 2 + y ** 2) + 1)
                 for s1 in [-1, 1]
                 for s2 in [-1, 1]],
               (x, -5, 5),
               (y, -5, 5),
               show=False)
    p.xlabel = x
    p.ylabel = y
    p.save(f'sample4_{i}.png')

p.show()

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

% ./sample3.py -v
3.
test1 (__main__.MyTestCase) ... ok
test2 (__main__.MyTestCase) ... ok

----------------------------------------------------------------------
Ran 2 tests in 0.343s

OK
%

0 コメント:

コメントを投稿