Advertisement
Guest User

Untitled

a guest
Mar 4th, 2015
187
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.77 KB | None | 0 0
  1. using namespace arma;
  2.  
  3. vec eigval;
  4. mat eigvec;
  5.  
  6. // C++11 initialization
  7. mat A = { 1, -1, 0, -1, 2, -1, 0, -1, 1};
  8.  
  9. eig_sym(eigval, eigvec, A);
  10. eigvec.print("Eigen Vectors");
  11.  
  12. Eigen Vectors
  13. -5.7735e-01 -7.071068e-01 0.4082483
  14. -5.7735e-01 9.714451e-e17 -0.8164966
  15. -5.7735e-01 7.017068e-01 0.4082483
  16.  
  17. import numpy as np
  18. w,v = np.linalg.eig(np.array([[1,-1,0],[-1,2,-1],[0,-1,1]]))
  19. v
  20.  
  21. array([[ -4.08248290e-01, -7.07106781e-01, 5.77350269e-01],
  22. [ 8.16496581e-01, 2.61214948e-16, 5.77350269e-01],
  23. [ -4.08248290e-01, 7.07106781e-01, 5.77350269e-01]])
  24.  
  25. eigen(matrix(c(1,-1,0,-1,2,-1,0,-1,1), 3, byrow=TRUE)$vectors
  26.  
  27. -4.082483e-01 -7.071068e-01 5.773503e-01
  28. 8.164966e-01 9.420555e-16 5.773503e-01
  29. -4.082483e-01 7.071068e-01 5.773503e-01
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement