• Remix
  • Share
  • New Entry

  • haz

  • /
  • Corn Lorenz Attractor

on 30 Oct 2021
  • 3
  • 13
  • 0
  • 0
  • 273
figure(Color=[0 0.7 0])
% Parameters
r=28;
s=10;
b=8/3;
% Timespan and initial position
% T = [0 1];
P0=ones(3,1);
% Jacobian matrix as anonymous function
J = @(T,P)[s*(P(2)-P(1)); r*P(1)-P(2)-P(1)*P(3); P(1)*P(2)-b*P(3)];
% Solving through ODE45
for i=1:100
T=[0 101-i];
[T,P]=ode45(J,T,P0);
% Plot
plot3(P(:,1),P(:,2),P(:,3),'LineWidth',5);
hold on
end
view(45,5);
[bN, bD] = rat(b);
axis equal off
colororder(flipud(parula))
Remix Tree
Load full remix tree