Animation of the Lorenz Attractor

Hi everybody,
I'm supposed to modify this code so there's a single dot tracing the plot. I think I'm pretty close I just can't completely get it.
function loren3
clear;clf
global A B R
A = 10;
B = 8/3;
R = 28;
u0 = 100*(rand(3,1) - 0.5);
[t,u] = ode45(@lor2,[0,100],u0);
N = find(t>10); v = u(N,:);
x = v(:,1);
y = v(:,2);
plot(x,y);
view(90,90)
M(1)=getframe;
dt=1/100;
for j = 2:1000
t(j)=t(j-1)+dt;
plot(x*t(j),y*t(j),'o','MarkerFaceColor','b','MarkerSize',8)
M(j)=getframe;
if y<=0, break, end
end
pause
movie(M,1)
function uprime = lor2(t,u)
global A B R
uprime = zeros(3,1);
uprime(1) = -A*u(1) + A*u(2);
uprime(2) = R*u(1) - u(2) - u(1)*u(3);
uprime(3) = -B*u(3) + u(1)*u(2);

1 Comment

Jan
Jan on 9 Dec 2016
Edited: Jan on 9 Dec 2016
The youtube link is not working for me, so I cannot guess,what you want to change. It would be efficient, if you explain this directly instead of letting the readers get this most important detail of your question by using an external web service.

Sign in to comment.

 Accepted Answer

KSSV
KSSV on 9 Dec 2016
doc comet

2 Comments

How do I make the comet head go slower?
Read the documentation...you can mention time factor as input.

Sign in to comment.

More Answers (0)

Tags

No tags entered yet.

Asked:

on 9 Dec 2016

Commented:

on 9 Dec 2016

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!