% Initialisation des tableaux :
% V = vitesses et P = positions. 
% NOMS contient le nom des étoiles. 
% C donne leurs couleurs.
% M = masses des étoiles.

% 18 étoiles
% les plus proches du système solaire)

% jbr. 9-11-97

V=[-5.69 -8.75 -2.78 7.32 2.38 -.87 4.6 -.82 1.89 2.08 2.51 0 .52 -1.95 -6.75 4.7 1.79 -3.54
   4.76 1.13 22.03 -.47 .75 24.2 .69 9.86 -2.21 10.8 -2.32 0 -6.62 4.68 10.81 6.16 1.36 17.71
   3.35 -15.45 .02 -20.11 -3.65 16.78 -.5 -5.94 -2.59 -.41 -4.09 0 3.92 4.51 10.56 .51 -.11 2.28];
P=[68 464 394 -404 -295 -7 408 -462 -98 487 -683 0 646 -106 608 718 111 334
   -365 -42 -377 107 658 -371 534 136 514 219 44 0 307 -86 -235 227 -536 -194
   631 450 433 307 68 30 -114 62 -157 -175 13 0 -208 -243 -182 -233 -241 -594];
noms=str2mat('Struve 2398','Ross 248','61 Cygni','Lalande 21185',...
	'Procyon 5', 'Etoile de Barnard','Epsilon eridani','Wolf 259',...
	'Sirius','Luyten 726-8','Ross 128','Soleil','Tau ceti','Alpha du centaure',...
	'Luyten 789-6','Luyten 725-32','Ross 154','Epsilon indi');
% couleurs artificelles
C=[1 0 0
   1 0 0
   .5 0 0
   1 0 0 
   0 0 1
   1 0 0
   .5 0 0
   1 0 0
   0 0 1
   1 0 0
   1 0 0
   1 1 0	
   1 0 0    
   0 1 1	
   0 0 1	
   1 0 1	
   0 1 0	
   .5 0 0];
M=[.26 .17 .69 .39 1.29 .21 .74 .1 2.96 .19 .21 1 .85 1.03 .13 .21 .24 .69];
ind=12:17;
amas(P(:,ind),V(:,ind),M(ind),C(ind,:),2000,[-1000 2000 -1000 2000],2000,1)

   
