Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- evenper[x_] := Select[Permutations[x], Signature[#] == 1 &]
- evenper[x_] := Select[Permutations[x], Signature[#] == 1 &]
- manual=evenper[{a, b, c, d}]
- phi = 1.61803398875;
- p1 = manual /. {a -> phi, b -> 1, c -> 1/phi, d -> 0};
- p2 = manual /. {a -> phi, b -> 1, c -> -1/phi, d -> 0};
- p3 = manual /. {a -> phi, b -> -1, c -> 1/phi, d -> 0};
- p4 = manual /. {a -> phi, b -> -1, c -> -1/phi, d -> 0};
- p5 = manual /. {a -> -phi, b -> 1, c -> 1/phi, d -> 0};
- p6 = manual /. {a -> -phi, b -> 1, c -> -1/phi, d -> 0};
- p7 = manual /. {a -> -phi, b -> -1, c -> 1/phi, d -> 0};
- p8 = manual /. {a -> -phi, b -> -1, c -> -1/phi, d -> 0};
- list1=Union[p1, p2, p3, p4, p5, p6, p7, p8]//Sort;
- evenper[x_] := Select[Permutations[x], Signature[#] == 1 &]
- phi = 1.61803398875;
- list20 = #*{phi, 1, 1/phi, 0} & /@ (Tuples[{{1, -1}, {1, -1}, {1, -1}, {1}}]);
- list2 = Flatten[evenper[#] & /@ case30, 1] // Sort;
- Position[list1, {-1.618033988749895`, 0, 0.6180339887498948`, -1}]
- Position[list2, {-1.618033988749895`, 0, 0.6180339887498948`, -1}]
- (*{}*)
- (*{{6}}*)
- Union[{2., (Sqrt[5] + 1)/2}]
- % // N
- Union[{2., (Sqrt[5] + 1.0)/2}]
- SortBy[{2., (Sqrt[5] + 1)/2}, N]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement