Solution P3.3.3
The following program is hopefully self-explanatory.
import numpy as np
import matplotlib.pyplot as plt
# Number seeds in the sunflower seedhead
n = 200
# The golden ratio, 1.61803398874989...
phi = (1 + np.sqrt(5))/2
# Calculate and plot the seed positions, (r, theta) for each seed, s
s = np.linspace(0, n-1, n)
r = np.sqrt(s)
theta = 2 * np.pi * s / phi
plt.polar(theta, r, marker='o', linestyle='', markersize=12, color='r')
plt.show()
