- /
-
Paleo Climatology - Tree rings
on 18 Oct 2021
- 40
- 9
- 0
- 0
- 272
N=50;
x=2+rand(1,N)-1;
y=2+rand(1,N)-1;
P=[x;y];
P=P-mean(P,2);
t=atan2(P(2,:),P(1,:));
[t,idx]=sort(t);
P=P(:,idx);
P=[P P(:,1)];
figure('color','k')
hold on
for k=1:0.001:1.5
S1=[k 0;0 1];
S2=[1 0;0 k];
P=S1*(S2*P);
plot(P(1,:),P(2,:),'color',[0.8 0.4 0],'LineWidth',2);
end
axis([-2 2 -2 2])
axis off