Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- >>> from numpy import matrix
- >>> from numpy import linalg
- >>> Q = matrix([[ 1. , -0.25 , 0. , 0. ,
- ... 0. , 0. , 0. , 0. , 0. ,
- ... -0.25 ], [-0.25 , 1. , -0.25 , 0. ,
- ... 0. , 0. , 0. , 0. , 0. ,
- ... 0. ], [ 0. , -0.25 , 1. , -0.25 ,
- ... 0. , 0. , 0. , 0. , 0. ,
- ... 0. ], [ 0. , 0. , -0.36111111, 1. ,
- ... -0.13888889, 0. , 0. , 0. , 0. ,
- ... 0. ], [ 0. , 0. , 0. , -0.25 ,
- ... 1. , -0.25 , 0. , 0. , 0. ,
- ... 0. ], [ 0. , 0. , 0. , 0. ,
- ... -0.13888889, 1. , -0.36111111, 0. , 0. ,
- ... 0. ], [ 0. , 0. , 0. , 0. ,
- ... 0. , -0.25 , 1. , -0.25 , 0. ,
- ... 0. ], [ 0. , 0. , 0. , 0. ,
- ... 0. , 0. , -0.36111111, 1. , -0.13888889,
- ... 0. ], [ 0. , 0. , 0. , 0. ,
- ... 0. , 0. , 0. , -0.13888889, 1. ,
- ... -0.36111111], [-0.13235294, 0. , 0. , 0. ,
- ... 0. , 0. , 0. , 0. , -0.36764706,
- ... 1. ]])
- >>>
- >>> a = [0.5, 0.5, 0.5117647058823529, 0.6555555555555556, 0.5,
- ... 0.6555555555555556, 0.525, 0.6555555555555556, 0.6555555555555556,
- ... 0.5117647058823529]
- >>> p=a*linalg.inv(Q)
- >>> p*Q-a
- matrix([[ 2.22044605e-16, -1.66533454e-16, 0.00000000e+00,
- 0.00000000e+00, -5.55111512e-17, 2.22044605e-16,
- -2.22044605e-16, 0.00000000e+00, 0.00000000e+00,
- -1.11022302e-16]])
- >>> p=linalg.solve(Q,a)
- >>> p*Q-a
- matrix([[ 0.12919174, 0. , -0.13265237, -0.00413775, 0.26823418,
- 0.00397455, -0.27302125, 0.12658865, -0.00717732, -0.11324618]])
- >>>
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement