Corn Lorenz Attractor

on 30 Oct 2021

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))

Load full remix tree