- /
-
Roller Coaster
on 23 Oct 2021
- 7
- 75
- 0
- 0
- 277
%% Attractor twist roller coaster
t = 0:.05:80;
%% solving dynamics
[~,o]=ode45(@f,t,[.92 -1 0]); %Initical conditon ? .1 0 0
%set(gca,'Color','k');
%Gray scale, [black, white, black, white, black, white]
g=gray;h=flip(g);
C=[h;g;h;g;h;g];
C=C(20:1475,:);
colormap(C);
shading interp
patch(o(:,1),o(:,2),o(:,3),t,EdgeC='b',LineW=4);
view(-40,20)
axis off
%% odinary differential
function S=f(t,X)
x=X(1);
y=X(2);
z=X(3);
a=.2;b=4;%
E=8; d=1;%
S=[-y^2-z^2-a*x+a*E;x*y-b*x*z-x+d;b*x*y+x*z-z];
end