• Remix
  • Share
  • New Entry

on 15 Oct 2021
  • 7
  • 16
  • 1
  • 0
  • 267
figure('color','white');
hold on
l=2;
r=-l:.04:l;
axis equal
axis off
q=[-1.5 1.5];
xlim(q)
ylim(q)
D=@(t,x)[x(2)-x(1)*(x(1)^2+x(2)^2-1)^2;-x(1)-x(2)*(x(1)^2+x(2)^2-1)^2];
for x=r
for y=r
[~,o]=ode23(D,[0 50],[x;y]);
plot(o(:,1),o(:,2),'color',[x x*y/l y]./(l*2)+.5,'linewidth',2)
end
end
Remix Tree