function amas(P,V,M,c,nbiter,ax,fact,dt)
%Simulation d'un amas d'étoiles
% La fonction AMAS(T) représente la simulation des amas d'étoiles
% décrites par leurs positions P, leurs vitesses V et leurs masses M.
% C est une table de couleurs.
% NBITER indique la durée de la simulation. 
% AX donnent les axes
% FACT donne le facteur d'échelle pour l'accélération.
% DT est le pas pour le temps.

% jbr. 9-11-97.

[li,nbe]=size(P);
A=zeros(3,nbe);
g=9.81;				% gravitation

cla
for i=1:nbe
	h(i)=plot(P(1,i),P(2,i),'.','markersize',4,'erasemode','none','color',c(i,:));
	hold on
end
axis(ax)
% axis on
axis off
figure(gcf)

for t=1:nbiter
	for i=1:nbe
		A(:,i)=0;
		j=[1:i-1,i+1:nbe];
		d=[P(1,i)-P(1,j);P(2,i)-P(2,j);P(3,i)-P(3,j)];
		d=sqrt(sum(d.^2));
		f=g*M(i)*M(j)./(d.^2);
		A(1,i)=sum(f.*(P(1,j)-P(1,i))./d);
		A(2,i)=sum(f.*(P(2,j)-P(2,i))./d);
		A(3,i)=sum(f.*(P(3,j)-P(3,i))./d);
	end
	V=V+dt*A*fact;
	P=P+dt*V;
	for k=1:nbe
		set(h(k),'XData',P(1,k),'YData',P(2,k))
	end
	drawnow
end

hold off
			
		
		
		
			
