回答編集履歴
1
追記
test
CHANGED
|
@@ -1,4 +1,6 @@
|
|
|
1
1
|
numpyは使う必要はありません。
|
|
2
|
+
|
|
3
|
+
基本的にはsympyで代数的な処理を行った後、最終的な結果をnumpyに渡して処理するという形になります。
|
|
2
4
|
|
|
3
5
|
```python
|
|
4
6
|
|
|
@@ -25,3 +27,29 @@
|
|
|
25
27
|
print(C)
|
|
26
28
|
|
|
27
29
|
```
|
|
30
|
+
|
|
31
|
+
ついでに、kはいらないですね。
|
|
32
|
+
|
|
33
|
+
```python
|
|
34
|
+
|
|
35
|
+
import sympy as sym
|
|
36
|
+
|
|
37
|
+
|
|
38
|
+
|
|
39
|
+
φ = sym.Symbol('φ')
|
|
40
|
+
|
|
41
|
+
|
|
42
|
+
|
|
43
|
+
C = sym.Matrix([[φ, 0],
|
|
44
|
+
|
|
45
|
+
[0, 0]])
|
|
46
|
+
|
|
47
|
+
for i in range(2):
|
|
48
|
+
|
|
49
|
+
for j in range(2):
|
|
50
|
+
|
|
51
|
+
C[i, j] += 8*φ
|
|
52
|
+
|
|
53
|
+
print(C)
|
|
54
|
+
|
|
55
|
+
```
|