- /
-
Cold Lorenz Attractor
on 30 Oct 2021
- 5
- 28
- 1
- 0
- 256
figure(Color='k')
% 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(jet)
