Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- b = 10;
- n = 30;
- l = First@RealDigits[Pi, b, n];
- mask = RegionPlot[x^2 + y^2 >= 1, {x, -1, 1}, {y, -1, 1},
- PlotStyle -> White, BoundaryStyle -> None];
- frame[z_] := Show[
- Graphics[
- {
- Table[
- If[Exp[z] Sqrt[(b^-i) ] < 1,
- a = l[[i]];
- {Darker[Hue[a/(b - 1)], 0.4],
- Disk[{0, 0}, Exp[z] Sqrt[(b^(-i + 1)) ], {0, -a 2 Pi/(b - 1)}],
- White, Disk[{0, 0}, Exp[z] Sqrt[(b^-i) ]], White,
- Text[Style[a, Medium, Bold,
- FontSize -> 50 Exp[z] Sqrt[(b^(-i + 1)) ]],
- Exp[z] (Sqrt[(b^(-i + 1)) ] + Sqrt[(
- b^-i) ]) {Cos[Pi/(b - 1)], -Sin[Pi/(b - 1)]}/2]
- }, {}]
- , {i, Length@l}],
- Thick,
- RGBColor[0.8, 0.15, 0.2],
- Table[
- Line[{{0, 0}, {Cos[\[Theta]], Sin[\[Theta]]}}]
- , {\[Theta], 0, 2 Pi, 2 Pi/(b - 1)}]
- }, PlotRange -> 1], mask
- ];
- Manipulate[frame[z],{z,0,34}]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement