In this research we model the distribution of mass in simulated galaxies by solving the Vlasov-Poisson system of equations. We've recently expanded our simulations to include multiple species of matter. This allows us to visualize the individual spatial density distributions of, for example, stars and dark matter as well as the joint gravitational potential. We have developed a library of
computational tools to allow us to investigate a number of the physical properties of these galaxies. In future work we will use these tools to compare the characteristics of our model galaxies to those of observed galaxies.