読者です 読者をやめる 読者になる 読者になる

tak0kadaの何でもノート

発声練習、生存確認用。

医学関連は 医学ノート

sympyで行列計算

python 線形代数 sympy

物理数学の教科書を読んでいて文字入りの逆行列の計算が合わなかったのでsympyを使った。sympyがインストールされていなければSymPy Liveが少し遅いが便利。

from sympy import *
l, m, s = symbols("λ μ σ")
A = Matrix([[l + 2 * m, l, l], [l, l + 2 * m, l], [l, l, l + 2 * m]])
y = Matrix([s, 0, 0])
x = A.inv() * y
print( x.symplify() )
Matrix([
[σ*(λ + μ)/(μ*(3*λ + 2*μ))],
[   -λ*σ/(2*μ*(3*λ + 2*μ))],
[   -λ*σ/(2*μ*(3*λ + 2*μ))]])