2022-06-12, 16:27 | #1 |
Jun 2021
3×17 Posts |
Ulam spirale, envolved to
Hi.
Plot code was taken from here https://community.wolfram.com/groups/-/m/t/1543056 Code:
mgonLayer[q_, m_] := Ceiling[1/2 (1 + Sqrt[-8 + m + 8 q]/Sqrt[m])] - 1 maround[q_, m_] := q - (1 + m/2 ((mgonLayer[q, m] - 1) (mgonLayer[q, m]))) mpoint[q_, m_] := (mgonLayer[q, m] AngleVector[(2 \[Pi])/m (#1 - 1)] + #2 2 Sin[\[Pi]/m] AngleVector[2 \[Pi]/m (#1 - 1) + \[Pi]/2 + \[Pi]/m]) & @@ Quiet[QuotientRemainder[maround[q, m], mgonLayer[q, m]]] ListPlot[mpoint[#, M] & /@ (Prime[#] & /@ Range[10000]), Axes -> False, AspectRatio -> Automatic] Lets make a little changes in the plot code: Code:
ListPlot[mpoint[#, M] & /@ (Prime[#*#] & /@ Range[20000]), Axes -> False, AspectRatio -> Automatic] Sharingan https://www.mersenneforum.org/attach...1&d=1655050840 M = 503 Rinnegan https://www.mersenneforum.org/attach...1&d=1655051029 Last fiddled with by RomanM on 2022-06-12 at 16:33 Reason: Range[20000] for obtaine the same pictures |
2022-06-18, 19:57 | #2 |
Jun 2021
3×17 Posts |
Hello.
Code:
ListPlot[mpoint[#, 4] & /@ ((Prime[2*#] - 2*Prime[#]) & /@ Range[15000]), Axes -> True, AspectRatio -> Automatic] Code:
Prime[2*#] - 2*Prime[#] |
2022-06-23, 12:20 | #3 |
Jun 2021
33_{16} Posts |
some symmery
Code:
ListPlot[mpoint[#, 4] & /@ ((Prime[#]/EulerPhi[#]) & /@ Range[15000]), Axes -> True, AspectRatio -> Automatic] |
Thread Tools | |
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
RH and Ulam's spiral | Damian | Lounge | 0 | 2009-09-26 03:17 |