Advertisement
Guest User

Untitled

a guest
Sep 19th, 2014
241
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.14 KB | None | 0 0
  1. evenper[x_] := Select[Permutations[x], Signature[#] == 1 &]
  2.  
  3. evenper[x_] := Select[Permutations[x], Signature[#] == 1 &]
  4. manual=evenper[{a, b, c, d}]
  5. phi = 1.61803398875;
  6. p1 = manual /. {a -> phi, b -> 1, c -> 1/phi, d -> 0};
  7. p2 = manual /. {a -> phi, b -> 1, c -> -1/phi, d -> 0};
  8. p3 = manual /. {a -> phi, b -> -1, c -> 1/phi, d -> 0};
  9. p4 = manual /. {a -> phi, b -> -1, c -> -1/phi, d -> 0};
  10. p5 = manual /. {a -> -phi, b -> 1, c -> 1/phi, d -> 0};
  11. p6 = manual /. {a -> -phi, b -> 1, c -> -1/phi, d -> 0};
  12. p7 = manual /. {a -> -phi, b -> -1, c -> 1/phi, d -> 0};
  13. p8 = manual /. {a -> -phi, b -> -1, c -> -1/phi, d -> 0};
  14. list1=Union[p1, p2, p3, p4, p5, p6, p7, p8]//Sort;
  15.  
  16. evenper[x_] := Select[Permutations[x], Signature[#] == 1 &]
  17.  
  18. phi = 1.61803398875;
  19. list20 = #*{phi, 1, 1/phi, 0} & /@ (Tuples[{{1, -1}, {1, -1}, {1, -1}, {1}}]);
  20. list2 = Flatten[evenper[#] & /@ case30, 1] // Sort;
  21.  
  22. Position[list1, {-1.618033988749895`, 0, 0.6180339887498948`, -1}]
  23. Position[list2, {-1.618033988749895`, 0, 0.6180339887498948`, -1}]
  24. (*{}*)
  25. (*{{6}}*)
  26.  
  27. Union[{2., (Sqrt[5] + 1)/2}]
  28.  
  29. % // N
  30.  
  31. Union[{2., (Sqrt[5] + 1.0)/2}]
  32.  
  33. SortBy[{2., (Sqrt[5] + 1)/2}, N]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement