• Remix
  • Share
  • New Entry

on 23 Oct 2021
  • 1
  • 9
  • 0
  • 0
  • 260
X=[];
for m=2:-0.03:0.01
f=@(t,y)([y(2);m*(1-y(1)^2)*y(2)-y(1)]);
[t,y]=ode23(f,0:0.01:19.99,[2; 0]);
X=[X;[y,zeros(size(y,1),1)-m+10]];
hold on
end
r=@(x)(reshape(x(:,1),2000,[]));
surf(r(X(:,1)),r(X(:,2)),r(X(:,3)),EdgeColor='none')
axis off
daspect([1 1 1])
colormap winter
Remix Tree