- /
-
Sunflower
on 13 Oct 2021
- 5
- 55
- 2
- 0
- 279
M = 120;
N = 90;
A = 0:4.0449:360;
X = NaN(3,N);
Y = NaN(3,N);
C = hsv(M);
for m = 1:M
k = mod(m,2);
s = sign(0.5-k);
a = 2.0225*k;
r = 0.97^m;
X = [X(2:3,:); r*cosd(A+a)];
Y = [Y(2:3,:); r*sind(A+a)];
P = @(X)[X(1,:); X(2,:); X(3,:); circshift(X(2,:),s,2)];
patch(P(X),P(Y),C(m,:));
hold on;
end
axis equal off