% Brownian motion %
X=randn(2000,3); S=zeros(2000,3);
for i=2:2000; S(i,:)=sum(X(1:i,:)); end;
comet(S(:,1),S(:,2)); pause
comet3(S(:,1),S(:,2),S(:,3));