function tournesol(nb)

n=1:nb;						% Nombre de points
nbdor=(1+sqrt(5))/2;		% Nombre d'or
alpha=2*pi/(nbdor+1);		% alpha=2*pi/(nbdor^2)=2*pi/(nbdor+1)=2*pi-2*pi/nbdor=137.5078°

% Courbes en coordonnées polaires
r=sqrt(n);					% rayon
t=alpha.*n;					% angle polaire

x=r.*cos(t);				% conversion polaire->cartésien
y=r.*sin(t);

% représentation des points
plot(x,y,'.','markersize',12)
axis equal
axis off
