数学 - Python - 代数学 - 連立方程式と高次方程式 - 連立2元1次方程式の解、不定、係数

2020年5月4日月曜日

数学 - Python - 代数学 - 連立方程式と高次方程式 - 連立2元1次方程式の解、不定、係数

学習環境

代数への出発 (新装版数学入門シリーズ) (松坂和夫(著)、岩波書店)の第5章(連立方程式と高次方程式)、問2の解答を求めてみる。

1. $\begin{array}{l} {\begin{cases} - 2 a x - 6 y = - 10 \\ 4 x - 6 y = b \end{cases} \\ - 2 a = 4 \\ b = - 10 \\ a = - 2 \end{array}$
2. $\begin{array}{l} {\begin{cases} 3 x - 6 a y - 42 = 0 \\ 3 x + (5 - a) y - b = 0 \end{cases} \\ - 6 a = 5 - a \\ a = - 1 \\ b = 42 \end{array}$

コード

#!/usr/bin/env python3
from unittest import TestCase, main
from sympy import symbols, solve

print('2.')

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


class TestCube(TestCase):
    def test1(self):
        self.assertEqual(
            len(solve([eq.subs({a: -2, b: -10})
                       for eq in [a * x + 3 * y - 5,
                                  4 * x - 6 * y - b]])),
            1)

    def test2(self):
        self.assertEqual(
            len(solve([eq.subs({a: -1, b: 42})
                       for eq in [x / 2 - a * y - 7,
                                  3 * x + (5 - a) * y - b]])),
            1)


if __name__ == "__main__":
    main()

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

% ./sample2.py -v
2.
test1 (__main__.TestCube) ... ok
test2 (__main__.TestCube) ... ok

----------------------------------------------------------------------
Ran 2 tests in 0.305s

OK
%

0 コメント:

コメントを投稿

コメントの投稿(Atom)