#declare max_particles = 25;
        #declare space = 20;
        #declare setup = 1;
        #declare xstart = 0-(sqrt(max_particles)/2)*space;
        #declare ystart = xstart;
        #declare num = 0;
        #declare step_size = 0.3;
        #declare grav_const = 1;
        #declare perimeter = 10000;
        #declare frames = 50;
        #declare alive = array[max_particles+1]
        #declare xvec = array[max_particles+1]
        #declare yvec = array[max_particles+1]
        #declare xpos = array[max_particles+1]
        #declare ypos = array[max_particles+1]
        #declare mass = array[max_particles+1]
        #declare rad = array[max_particles+1]

        #declare rdm = seed(12345);

        #while (num < max_particles)


                        #declare alive[num] = 1;
                        #declare xpos[num] = xstart;
                        #declare ypos[num] = ystart;
                        #declare xvec[num] = (rand(rdm)*100)-50;
                        #declare yvec[num] = (rand(rdm)*100)-50;
                        #declare mass[num] = 10;
                        #declare rad[num] = 1;


                        #declare setup = setup + 1;
                        #declare ystart = ystart + space;

                        #if (setup > sqrt(max_particles))
                                #declare ystart =
0-(sqrt(max_particles)/2)*space;
                                #declare xstart = xstart + space;
                                #declare setup = 1;
                        #end



                        #declare num = num + 1;




        #end

        #declare write_part = 1;

        #fopen particle_text "particle.txt" write
        #write (particle_text ,max_particles,"\n",)
        #write (particle_text ,step_size,"\n",)
        #write (particle_text ,grav_const,"\n",)
        #write (particle_text ,perimeter,"\n",)
        #write (particle_text ,frames,"\n",)
        #while (write_part < max_particles)

                #write (particle_text ,alive[write_part],"\n",)
                #write (particle_text ,xvec[write_part],"\n",)
                #write (particle_text ,yvec[write_part],"\n",)
                #write (particle_text ,xpos[write_part],"\n",)
                #write (particle_text ,ypos[write_part],"\n",)
                #write (particle_text ,mass[write_part],"\n",)
                #write (particle_text ,rad[write_part],"\n",)

                #declare write_part = write_part + 1;

        #end

        #fclose particle_text

Post a reply to this message