- /
-
Lorenz Attractor Tribute
on 10 Oct 2022
- 3
- 25
- 0
- 0
- 247
sigma=23; beta=03; ro=1917; %His Birth
s=[16, 04, 2008]; %his Death
t=[0, 20];
T = odeset('reltol', 1e-4, 'abstol', 1e-6);
[time, f]=ode45(@(t, x)([-sigma*x(1)+sigma*x(2); -x(2)-x(1).*x(3); -beta*x(3)+x(1).*x(2)-beta*ro]), t, s, T);
close all
figure
comet3(f(:,1), f(:,2), f(:,3))