Sunflower 0.0.2 (with comments)

on 15 Oct 2021

hold on;axis off;

% Petals

t=linspace(0,2*pi,1e3);% theta vector

for b=0:pi/2:pi% angle of rotation

x=275*cos(6*t+b).*cos(t+b);% amplitude = 250

y=275*cos(6*t+b).*sin(t+b);% number of petals = 6 (odd)

c=sqrt(x.^2+y.^2)>127.38;% max radius florets = sqrt(118 * 137.5) ~ 127.38

patch(x(c),y(c),[1,.8157,.2510])% color_petals = [1,.8157,.2510]

end

% Florets

z=linspace(1,118,1200)*137.508;% golden angle = (3-sqrt(5)) * 180 ~ 137.5 [deg]

v=sqrt(z).*cos(z);

w=sqrt(z).*sin(z);

plot(v,w,'ok','MarkerFaceColor',[.4431,.2588,.0667])% color_florets = [.4431,.2588,.0667]