Schauen Sie sich diese Kamillenblüte an:

Schön, nicht wahr? Was ist, wenn ich dir sage, dass dies nicht wirklich eine Blume ist?

Viele Blumen (einschließlich Sonnenblumen, Kamille, Gänseblümchen und andere) bestehen tatsächlich aus vielen sehr kleinen Blüten (die schwarzen Punkte auf Sonnenblumen) auf einem Blütenkopf. Diese Miniatur - Blumen genannt Röschen , und sie sind in einer ganz besonderen Art und Weise angeordnet sind .

Grundsätzlich ist die Position des n-ten Blütchens auf einem Blütenkopf (in Polarkoordinaten):

Dabei ist c = 1 (Beachten Sie, dass 137,508 Grad = goldener Winkel. Sie müssen nicht genau diese Genauigkeit verwenden.)

Dies führt dazu, dass die Blütchen in einer Spirale namens Fermat-Spirale geformt werden. Die Positionierung der Blütchen hängt auch mit Fibonnaci-Zahlen zusammen, aber das ist eine Geschichte für ein anderes Mal.

Also, hier ist die Herausforderung. Berechnen Sie mit einer Ganzzahl n als Eingabe die Positionen der ersten n Florets und zeichnen Sie sie auf . Dies ist eine grafische Ausgabe , daher möchte ich, dass Sie die Punkte entweder in einem Fenster anzeigen oder als Daten in einem gängigen Bildformat an STDOUT oder in einer Datei ausgeben. Davon abgesehen sollte diese Herausforderung ziemlich einfach sein. Es ist Code-Golf , also gewinnt der kürzeste Code. VIEL GLÜCK UND VIEL SPASS!

Hier ist ein Beispielbild, wie eine Ausgabe aussehen könnte:

TI-BASIC, 34 Bytes

Für die Taschenrechner der Serie TI-83 + / 84 +.

Input N
2πe^(-2sinh⁻¹(.5→θstep
AnsN→θmax
"√(θ→r₁
Polar                      ;Displays polar mode graphs to graph screen
Dot                        ;Prevents lines from connecting points
DispGraph                  ;Displays graph screen

Dies betrachtet den Punkt am Ursprung als den 0. Punkt.

Dank des Ein-Byte- sinh⁻¹(Tokens 2πe^(-2sinh⁻¹(.5ist dies ein kurzer Weg, um den goldenen Winkel im Bogenmaß zu erhalten. Dies ergibt sich aus der Tatsache, dass e^(sinh⁻¹(.5es sich um den goldenen Schnitt handelt.

Hier sind Screenshots für N = 50.

(Ja, das ist ein 96x64-Monochrom-Display auf einem TI-84 +. Die neueren Farbrechner haben ein Auflösungs-Upgrade, aber immer noch nur 3,7% der Pixel eines iPhones.)

Drücken Sie TRACE, um die einzelnen Punkte zu durchlaufen.

Python 2, 85 82 81 Bytes

from pylab import*
for i in arange(0,input(),2.39996):polar(i,sqrt(i),'o')
show()

Verkürzt um ein Byte durch Marinus.

Verwenden Sie den goldenen Winkel im Bogenmaß. Die Bytelänge ist die gleiche, wenn ich stattdessen 137.508 verwendet habe, sieht aber irgendwie nicht so gut aus. Erzeugt eine Polarkurve mit pylab. Unten ist angegeben, wann 300 (für die ältere Version) die Eingabe ist und 7000 (für die neuere Version) die Eingabe ist. Kann den Winkel auf 2,4 runden, um die Anzahl der Bytes auf 77 zu verringern.

Hier ist eine längere Version, die durch Entfernen des Gitters und der Achse ein klareres Erscheinungsbild ergibt:

from pylab import *
def florets(n):
    for i in arange(0, n, 2.39996):polar(i, sqrt(i), 'o')
    grid(0)#turn off grid
    xticks([])#turn off angle axis
    yticks([])#turn off radius axis
    show()

Der Grund für die unterschiedlichen Farben liegt darin, dass jeder Punkt separat gezeichnet und als eigener Datensatz behandelt wird. Wenn die Winkel und Radien als Listen übergeben würden, würden sie als eine Menge behandelt und von einer Farbe sein.

Blitz 2D / 3D , 102 Bytes

(Die allererste Blitz 2D / 3D Antwort auf dieser Seite!)

Graphics 180,180
Origin 90,90
n=Input()
For i=1To n
t#=i*137.508
r#=Sqr(t)
Plot r*Cos(t),r*Sin(t)
Next

Eine Eingabe von 50füllt das Fenster. (Ja, ich könnte zwei Bytes abschneiden Graphics 99,99, aber das ist optisch nicht so interessant oder nützlich.)

Schönere Version (und schöner beendet):

Graphics 400,400
Origin 200,200

n=Input("How many florets? ")

For i = 1 To n
    t# = i * 137.508
    r# = Sqr(t)

    Oval r*Cos(t)-3,r*Sin(t)-3,7,7,1
Next

WaitKey
End

Mathematica, 43 42 Bytes

ListPolarPlot@Array[(2.39996#)^{1,.5}&,#]&

Dies ist eine unbenannte Funktion, die ein ganzzahliges Argument verwendet, z

^{Der Screenshot verwendet eine ältere Version, aber die Ausgabe sieht gleich aus.}

Mathematica verfügt über eine integrierte Funktion GoldenAnglefür noch genauere Ergebnisse, die jedoch länger dauert als 2.39996.

MATLAB, 42 Bytes

t=2.39996*(1:input(''));polar(t,t.^.5,'.')

Ruft die eingegebene Nummer ab und erstellt einen Bereich von 1 bis zu dieser Nummer.

Multipliziert den Bereich mit dem goldenen Winkel im Bogenmaß (der verwendete Wert liegt näher am wahren Wert als 137,508 Grad bis 6 sf).

Zeichnen Sie dann einfach Theta vs. R in einem Polarkoordinatendiagramm mit Punkten. Hier gezeigt mit 2000 Punkten

Ein etwas hübscher aussehendes Diagramm (keine Gitterlinien) wäre dieser Code:

t=2.39996*(1:input(''));[x,y]=pol2cart(t,t.^.5);plot(x,y,'.');axis equal

Dies geht jedoch zu Lasten von 31 Bytes. Auch hier ist es mit 2000 Punkten gezeigt

R, 58 55 54 Bytes

x=2.39996*1:scan();plotrix::radial.plot(x^.5,x,rp="s")

Dazu muss das plotrixPaket installiert werden, es muss jedoch nicht importiert werden, da wir explizit auf den Namespace verweisen.

Ungolfed:

# Read a number of STDIN
n <- scan()

x <- 2.39996*(1:n)

# The rp.type = "s" option specifies that we want to plot points rather
# than lines (the default)
plotrix::radial.plot(lengths = sqrt(x), radial.pos = x, rp.type = "s")

Beispielausgabe für n = 1500:

3 Bytes gespart dank Plannapus!

R, 55 54 Bytes

t=1:scan()*2.39996;r=t^.5;plot(r*cos(t),r*sin(t),as=1)

Hier ist das Ergebnis für n = 1000:

Bearbeiten: Dank @AlexA 1 Byte mit teilweiser Übereinstimmung der Argumente ( asanstelle von asp) gespeichert.!

R, 48 47 Bytes

Ich denke, das ist ausreichend anders als die anderen R-Lösungen bisher. Dieser verwendet komplexe Vektoren, um die Koordinaten zu konstruieren. Die Quadrate von t und t werden in die Modul- und Argumentparameter eingefügt, und die XY-Werte werden aus dem Real und dem Imaginären abgeleitet. Vielen Dank an @AlexA. für das Byte.

plot(complex(,,,t^.5,t<-1:scan()*2.39996),as=1)

HTML + JavaScript 179

<canvas id=C></canvas><script>n=1500,C.width=C.height=400,T=C.getContext('2d');for(n=prompt(M=Math);n--;)r=M.sqrt(t=n*2.4)*9,T.fillRect(M.cos(t)*r+200,M.sin(t)*r+200,2,2)</script>

Jolf, 25 Bytes

yG@@KyoΜzXDyOUH*Hm°yT'.}

(Ausgabe für n = 5000)

Probieren Sie es online aus. (Beachten Sie, dass die resultierende Spirale klein ist)

Nicht konkurrierend, da Jolf nach dieser Herausforderung geschaffen wurde. Dies sind 25 Byte, wenn es mit ISO-8859-7 codiert ist, und es enthält eines, das nicht gedruckt werden kann (hier ist ein Hexdump):

0000000: 7947 4096 404b 796f cc7a 5844 794f 5548  yG@.@Kyo.zXDyOUH
0000010: 2a48 6db0 7954 272e 7d                   *Hm.yT'.}

Erläuterung

yG@@KyoΜzXDyOUH*Hm°yT'.}
yG@@K                      goto (150,75) (center of the canvas)
     yo                    set current location as the origin
       MzX                 map over range 1...input
          D                start of function
           yO              goto polar coordinates ....
             UH            radius: square root of argument
               *Hm°        angle: argument times golden angle
                   yT'.    draw a dot there
                       }

Python 2, 74 Bytes

from pylab import*
i=arange(1,input(),2.39996)
polar(i,sqrt(i),'o')
show()

MATL , 20 Bytes (nicht konkurrierend)

Als nicht konkurrierend markiert, da die Sprache die Herausforderung nachträglich angibt

:2.4*tX^wJ*Ze*'.'&XG

Probieren Sie es bei MATL Online!

Der goldene Winkel 137.708deg = pi*(3-sqrt(5))rad = 2.39996...rad wird als 2.4rad angenähert .

In der folgenden Version ( 25 Byte ) wird der exakte Wert bis zur doubleGleitkommagenauigkeit verwendet:

:YPI5X^-**tX^wJ*Ze*'.'&XG

Probieren Sie es bei MATL Online!

Tcl / Tk, 114

grid [canvas .c]
proc P n {time {incr i
.c cr o [lmap h {cos sin cos sin} {expr sqrt($i*2.4)*$h\($i*2.4)+99}]} $n}

Anwendungsbeispiel:

P 1024

und gibt das Fenster aus