function Convergentfield
	title('Convergent field');

	a=10; r=28; b=8/3; dt=0.01;
	X=0.01; Y=0.01; Z=0.01;
	n=600;
	XYZ=zeros(n, 3);
	cols=1:n;
	t=0;

	for i=1:n
		X1=X;
		Y1=Y;
		Z1=Z;

		X=X1+(-a*X1+a*Y1)*dt;
		Y=Y1+(-X1*Z1+r*X1-Y1)*dt;
		Z=Z1+(X1*Y1-b*Z1)*dt;
		#t=t+dt
		XYZ(i,:)=[X,Y,Z];
		cols(i)=i/n;
	endfor

	x=XYZ(:,1);
	y=XYZ(:,2);
	z=XYZ(:,3);

	#disp (size(x));
	#disp (cols(1));

	scatter3(x, y, z, [], cols);
	#colormap(jet);
endfunction