Advertisement
Guest User

Untitled

a guest
Oct 29th, 2010
135
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 2.05 KB | None | 0 0
  1. >>> from numpy import matrix
  2. >>> from numpy import linalg
  3. >>> Q = matrix([[ 1.        , -0.25      ,  0.        ,  0.        ,
  4. ... 0.        , 0.        ,  0.        ,  0.        ,  0.        ,
  5. ... -0.25      ], [-0.25      ,  1.        , -0.25      ,  0.        ,
  6. ... 0.        , 0.        ,  0.        ,  0.        ,  0.        ,
  7. ... 0.        ], [ 0.        , -0.25      ,  1.        , -0.25      ,
  8. ... 0.        , 0.        ,  0.        ,  0.        ,  0.        ,
  9. ... 0.        ], [ 0.        ,  0.        , -0.36111111,  1.        ,
  10. ... -0.13888889, 0.        ,  0.        ,  0.        ,  0.        ,
  11. ... 0.        ], [ 0.        ,  0.        ,  0.        , -0.25      ,
  12. ... 1.        , -0.25      ,  0.        ,  0.        ,  0.        ,
  13. ... 0.        ], [ 0.        ,  0.        ,  0.        ,  0.        ,
  14. ... -0.13888889, 1.        , -0.36111111,  0.        ,  0.        ,
  15. ... 0.        ], [ 0.        ,  0.        ,  0.        ,  0.        ,
  16. ... 0.        , -0.25      ,  1.        , -0.25      ,  0.        ,
  17. ... 0.        ], [ 0.        ,  0.        ,  0.        ,  0.        ,
  18. ... 0.        , 0.        , -0.36111111,  1.        , -0.13888889,
  19. ... 0.        ], [ 0.        ,  0.        ,  0.        ,  0.        ,
  20. ... 0.        , 0.        ,  0.        , -0.13888889,  1.        ,
  21. ... -0.36111111], [-0.13235294,  0.        ,  0.        ,  0.        ,
  22. ... 0.        , 0.        ,  0.        ,  0.        , -0.36764706,
  23. ...           1.        ]])
  24. >>>
  25. >>> a = [0.5, 0.5, 0.5117647058823529, 0.6555555555555556, 0.5,
  26. ... 0.6555555555555556, 0.525, 0.6555555555555556, 0.6555555555555556,
  27. ... 0.5117647058823529]
  28. >>> p=a*linalg.inv(Q)
  29. >>> p*Q-a
  30. matrix([[  2.22044605e-16,  -1.66533454e-16,   0.00000000e+00,
  31.            0.00000000e+00,  -5.55111512e-17,   2.22044605e-16,
  32.           -2.22044605e-16,   0.00000000e+00,   0.00000000e+00,
  33.           -1.11022302e-16]])
  34. >>> p=linalg.solve(Q,a)
  35. >>> p*Q-a
  36. matrix([[ 0.12919174,  0.        , -0.13265237, -0.00413775,  0.26823418,
  37.           0.00397455, -0.27302125,  0.12658865, -0.00717732, -0.11324618]])
  38. >>>
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement