%
% LORENZ ATTRACTOR
%
% simulation data
  clear all
  xold=-5.0;yold=+5.0;zold=+5.0;
  a=10.0;b=27.0;c=2.50;delt=0.01;
% simulation loop
  for i=1:5000
  xnew=xold+delt*(a*yold-a*xold);
  ynew=yold+delt*(b*xold-xold*zold-yold);
  znew=zold+delt*(xold*yold-c*zold);
  xold=xnew;yold=ynew;zold=znew;
  x(i)=xnew;y(i)=ynew;z(i)=znew;  
  plot3(x,z,y)
  sketch(i)=getframe;
  end
  movie(sketch)