Here is a link to the program:

I am trying to model a galaxy so that I have multiple stars certain distances away from each other and with different masses.

The accretion and than cylinder expansion happens with just about every galaxy.

I have a few questions.

First My computer only runs this if the number of particles is in the 500s. If that's the case How can I increase the probability that I will when I get to it have a red giant?

Second how can I modify this so that I have a timescale during which white dwarfs form from some stars and you see in the single stars either no supernova or nova just degeneration or a type 2 supernova and in binaries a type 1a supernova which is like a type 2 supernova of 1 star when it is a red giant, that giant becoming a white dwarf, the other star becoming a red giant, and then the white dwarf sucking gas from the giant until it explodes and ejects its companion?

Third how can I make this compacted and taking up less space in case I want to expand it to a few galaxies or even start with a big bang?