Processing math: 75%

2017年12月13日水曜日

学習環境

数学読本〈5〉微分法の応用/積分法/積分法の応用/行列と行列式(松坂 和夫(著)、岩波書店)の第21章(もう1つの数学の基盤 - 行列と行列式)、21.2(行列式)、3次の行列式の諸性質、問29.を取り組んでみる。


    1. |0101-a2a1-a1-a1a-1|=-((1-a2)(a-1)-(1-a)2)=(1-a2)(1-a)+(1-a)2=(1-a)2(1+a+1)=(1-a)2(a+2)=(a-1)2(a+2)

    2. |1ab-aa10b10|=(b-a)(a-b)=-(a-b)2

    3. a(bc-a2)-b(b2-ac)+c(ab-c2)=abc-a3-b3+abc+abc-c3=3abc-a3-b3-c3

    4. | 1 0 0 a b - a c - a a 2 b 2 - a 2 c 2 - a 2 | = b - a c 2 - a 2 - c - a b 2 - a 2 = c - a b - a c - a - b + a = c - a b - a c - b = a - b b - c c - a

コード(Emacs)

Python 3

#!/usr/bin/env python3
from sympy import pprint, symbols, Matrix

a, b, c = symbols('a, b, c')
XS = [
    [a, 1, 1, 1, a, 1, 1, 1, a],
    [1, a, b, a, 1, 1, b, 1, 1],
    [a, b, c, b, c, a, c, a, b],
    [1, 1, 1, a, b, c, a ** 2, b ** 2, c ** 2]
]


for i, t in enumerate(XS, 1):
    print(f'({i})')
    X = Matrix(t).reshape(3, 3)
    d = X.det()
    for s in [X, d, d.factor()]:
        pprint(s)
        print()
    print()

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

$ ./sample29.py
(1)
⎡a  1  1⎤
⎢       ⎥
⎢1  a  1⎥
⎢       ⎥
⎣1  1  a⎦

 3          
a  - 3⋅a + 2

       2        
(a - 1) ⋅(a + 2)


(2)
⎡1  a  b⎤
⎢       ⎥
⎢a  1  1⎥
⎢       ⎥
⎣b  1  1⎦

   2            2
- a  + 2⋅a⋅b - b 

        2
-(a - b) 


(3)
⎡a  b  c⎤
⎢       ⎥
⎢b  c  a⎥
⎢       ⎥
⎣c  a  b⎦

   3              3    3
- a  + 3⋅a⋅b⋅c - b  - c 

             ⎛ 2                2          2⎞
-(a + b + c)⋅⎝a  - a⋅b - a⋅c + b  - b⋅c + c ⎠


(4)
⎡1   1   1 ⎤
⎢          ⎥
⎢a   b   c ⎥
⎢          ⎥
⎢ 2   2   2⎥
⎣a   b   c ⎦

   2      2        2      2    2        2
- a ⋅b + a ⋅c + a⋅b  - a⋅c  - b ⋅c + b⋅c 

-(a - b)⋅(a - c)⋅(b - c)


$

0 コメント:

コメントを投稿