Back to Gallery
  • Remix
  • Share
  • New Entry

  • Elkmay

  • /

  • Lorenz attractor - Paints

on 7 Oct 2021
  • 11
  • 37
  • 3
  • 0
  • 224
% 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
Remix Tree
Load full remix tree