#declare N_Star=200; // Number of stars
#declare A_Star=array[N_Star][2] // Array storing all stars x,z coordonate
#declare Radius=.1;
#declare R1=seed (5789);
#declare R2=seed (624);
#declare N1=0; 
#declare NS=0;
#while (N1<N_Star)
        #declare X1=34*rand(R1)-8;
        #declare Z1=20*rand(R2)-18;
        #if (NS=0) // The first time
                #declare NS=1;
                #declare A_Star[NS-1][0]=X1;
                #declare A_Star[NS-1][1]=Z1;
                #declare N1=N1+1;
        #else
                #declare N2=0;
                #declare Bool=0;
                #while (N2< NS)
                        #declare Dist=vlength(<X1,80,Z1>-<A_Star[N2][0],80,A_Star[N2][1]>);
                        #if (Dist>2*Radius) // the min distance beetween 2 stars must be 2*Radius
                                #declare Bool=1;
                       #end
                       #declare N2=N2+1;
                #end
                #if (Bool=1)
                        #declare NS=NS+1;
                        #declare A_Star[NS-1][0]=X1;
                        #declare A_Star[NS-1][1]=Z1;
                        #declare N1=N1+1;
                #end
        #end
        #render concat("NS=",str(NS,0,0),"\n")

#end
#declare N1=0;
#while (N1<N_Star)
        sphere {
                <A_Star[N1][0],80,A_Star[N1][1]>,Radius
                pigment {color rgb 1}
        }
        #declare N1=N1+1;
#end
camera {
        location <9,40,-8>
        look_at <9,100,-8>
}
light_source { <10,200,10> color rgb 1}
light_source { <-10,-200,-10> color rgb 1}