%%
beta=8/3;
sig=10;
rho=28;
y = [-1 0 0]';
s = [-1.1 0 0]';

h = 0.01;

hold off
plot3(y(1),y(2),y(3),'.', 'MarkerSize',6,'MarkerEdgeColor','red');axis([-10 10 -10 10 -10 50])
hold on
plot3(s(1),s(2),s(3),'.','MarkerSize',6,'MarkerEdgeColor','green'); 

for i=1:4000
    A = [ -sig sig 0; rho -1 -y(1); y(2) 0 -beta ];
    S = [ -sig sig 0; rho -1 -s(1); s(2) 0 -beta ];
    ydot = A*y;
    sdot = S*s;
    y = y + h*ydot;
    s = s + h*sdot;

plot3(y(1),y(2),y(3),'.', 'MarkerSize',6,'MarkerEdgeColor','red'); 
plot3(s(1),s(2),s(3),'.','MarkerSize',6,'MarkerEdgeColor','green'); 
F(i) = getframe;    

end

