%Solution of the Lorenz ODE.
clear all, close all, clc

[t,y]=ode45('lorenzflow',[0 100], [1 1 1]);

slowdownfactor=30;

y1l=repmat(y(:,1)',slowdownfactor,1);
y2l=repmat(y(:,2)',slowdownfactor,1);
y3l=repmat(y(:,3)',slowdownfactor,1);

yslow1=reshape(y1l,length(y(:,1))*slowdownfactor,1);
yslow2=reshape(y2l,length(y(:,1))*slowdownfactor,1);
yslow3=reshape(y3l,length(y(:,1))*slowdownfactor,1);

comet3(yslow1,yslow3,yslow2,0.2)

