- /
-
Gordion
on 12 Oct 2021
- 10
- 57
- 0
- 0
- 240
n=40;
t=2*pi*(1:n)/n;
t=t(randperm(length(t)));
p=cos(t)+1i*sin(t);
delta=0.02;
npl=200;
pall=zeros(length(p),npl);
pall(:,1)=p;
for k=2:size(pall,2)
p=pall(:,k-1);
p2=circshift(p,1);
d=p2-p;
d=delta*d./abs(d);
pall(:,k)=p+d;
end
plot(pall',LineWidth=3)
axis square off