201702GraduateEx0500 - dknife/Graduate2017_Autumn GitHub Wiki

import numpy as np

v = np.array([1,2,3])
A = np.array(
    [
        [1,2,0],
        [1,1,1],
        [2,2,4]
    ])

print(v)
print(A)
b = A.dot(v)
print(b)
Ainv = np.linalg.inv(A)
print(Ainv)
print(Ainv.dot(b))

[1 2 3]

[[1 2 0]
 [1 1 1]
 [2 2 4]]

[ 5  6 18]

[[-1.   4.  -1. ]
 [ 1.  -2.   0.5]
 [-0.  -1.   0.5]]

[ 1.  2.  3.]