学習環境
- Surface
- Windows 10 Pro (OS)
- Nebo(Windows アプリ)
- iPad
- MyScript Nebo - MyScript(iPad アプリ(iOS))
- 参考書籍
ラング線形代数学(上) (ちくま学現文庫)(S.ラング (著)、芹沢 正三 (翻訳)、筑摩書房)の3章(行列)、1(1次方程式)、練習問題4の解答を求めてみる。
2020年1月31日金曜日
コード
#!/usr/bin/env python3
from unittest import TestCase, main
from sympy import symbols, solve, Rational, I
print('4.')
x, y, z = symbols('x, y, z', imag=True)
class MyTestCase(TestCase):
def test_a(self):
s = solve((I * x - 2 * y - 1, x + I * y - 2))
self.assertEqual(s[x].expand(), 4 + I)
self.assertEqual(s[y].expand(), -1 + 2 * I)
def test_b(self):
self.assertEqual(solve((2 * x + I * y - (1 + I) * z - 1,
x - 2 * y + I * z,
-I * x + y - (2 - I) * z - 1)),
{x: (4 - 7 * I) / 10, y: (1 - I) / 2, z: (-3 - 6 * I) / 10})
def test_c(self):
s = solve(((1 + I) * x - y,
I * x + y - (3 - I)))
self.assertEqual(s[x].expand(), ((1 - 7 * I) / 5).expand())
self.assertEqual(s[y].expand(), ((8 - 6 * I) / 5).expand())
def test_d(self):
s = solve((I * x - (2 + I) * y - 1,
x + (2 - I) * y - (1 + I)))
self.assertEqual(s[x].expand(), ((5 - I) / 6).expand())
self.assertEqual(s[y].expand(), ((-1 + 3 * I) / 6).expand())
if __name__ == '__main__':
main()
入出力結果(Zsh、PowerShell、Terminal、Jupyter(IPython))
% ./sample4.py -v
4.
test_a (__main__.MyTestCase) ... ok
test_b (__main__.MyTestCase) ... ok
test_c (__main__.MyTestCase) ... ok
test_d (__main__.MyTestCase) ... ok
----------------------------------------------------------------------
Ran 4 tests in 0.729s
OK
%
0 コメント:
コメントを投稿