• Remix
  • Share
  • New Entry

on 13 Oct 2021
  • 10
  • 17
  • 0
  • 0
  • 259
t = 0:5e-4:10;
[~,y]=ode15s(@f,t,[.1 0 0]);
%% plot
patch(y(:,1),y(:,2),-y(:,3),t,'EdgeC','w','linew',3);
set([gcf,gca],'Color','b');
colormap(flipud(prism));
view(110,33);
axis off
%% ordinary differential equation
function w2 = f(t,X)
x = X(1);
y = X(2);
z = X(3);
w2 = [200*(y-x)+x*z/2; 50*y-x*z;.833*z+x*y-.65*x*x];
end
Remix Tree