n=200; 	% nombre de points
s=.02;	% température ou vitesse
x=rand(n,1)-.5;	% n points aléatoires entre -1/2 et 1/2
y=rand(n,1)-.5;
h=plot(x,y,'.');
shg
axis([-1 1 -1 1])
axis square
%%%% axis off
set(h,'EraseMode','xor','Markersize',18)
i=1;
nbi=10;
option='nesauvepas';
while 1
	x=x+s*randn(n,1);
	y=y+s*randn(n,1);
	set(h,'XData',x,'YData',y)
	drawnow
	i=i+1;
	if strcmp(option,'sauve')
	 if ~rem(i,nbi)
		[X,MAP]=capture;
		no=['000',num2str(fix(i/nbi))];
		ln=length(no);
		gifwrite(X,MAP,['brown',no(ln-3:ln)]);
	end
	end
end
