- /
- 
        Total solar eclipse
        on 13 Oct 2021
        
        
 
    - 27
- 118
- 1
- 0
- 279
% Remix of original work by Eric Ogier
M=120;
N=90;
A=0:4.0449:360;
X=NaN(3,N);
Y=X;
C=copper(M);
figure(Color='k')
hold on
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,:),EdgeC='n');
end
axis equal off
camva(4)
Image
 
            

 

 
             
             
 
             
             
              