Advertisement
Matthen

Sierpinski

Jun 4th, 2011
392
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.73 KB | None | 0 0
  1. pinski[0] = { { {0, 0}, {1, 0}, {1/2, Sqrt[3]/2} } };
  2. f[m_] := 1/2 (-1 + 3^(1 + m));
  3. subtriangs[triang_] := {
  4. {triang[[1]], (triang[[2]] + triang[[1]])/
  5. 2, (triang[[3]] + triang[[1]])/2},
  6. {(triang[[1]] + triang[[2]])/2,
  7. triang[[2]], (triang[[3]] + triang[[2]])/2},
  8. {(triang[[1]] + triang[[3]])/2, (triang[[2]] + triang[[3]])/2,
  9. triang[[3]]}
  10. };
  11. pinski[n_] := pinski[n] = Module[{prev = pinski[n - 1], new},
  12. new =
  13. Join[prev,
  14. Flatten[Map[subtriangs, prev[[f[n - 2] + 1 ;; f[n - 1]]]], 1]]
  15. ];
  16. Manipulate[Graphics[
  17. Table[Join[{ColorData["AlpineColors"][(i + 1)/(n + 1)]},
  18. Map[Polygon, pinski[n][[f[i] + 1 ;; f[i + 1]]]]]
  19. , {i, -1, n - 1}]
  20. ]
  21. , {n, 0, 5, 1}]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement