- /
-
RGB Limit Cycle
on 15 Oct 2021
- 7
- 20
- 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
Image
