In [1]:
# Example by Simone Conradi in Python/NumPy: https://molab.marimo.io/notebooks/nb_jjSUUNEbg1u3k9obMzWgqb
import matplotlib.pyplot as plt
In [2]:
import apple
In [3]:
v=apple.jit("λN.{k ⟜ 2*𝜋%ℝN; ([cos.(k*x)]'(frange 0 (ℝN-1) N), [sin.(k*x)]'(frange 0 (ℝN-1) N))}")
In [4]:
step=apple.jit("λxs. { φ ← (1+√5)%2; ♭(xs [φ^2*x-φ*y]⊗ xs) }")
In [5]:
N=3
v_x,v_y=v(5)
for _ in range(N): v_x = step(v_x);v_y = step(v_y)
In [6]:
fig, ax=plt.subplots(figsize=(6,6))
fig.set_facecolor("k")
ax.set_facecolor("k")
ax.scatter(v_x,v_y,s=1,c="w")
plt.show()
