Many changes. Replaced incorrect spherical bessel function with regular bessel function, using the gsl_sf_bessel implementation. Nondimensionalized simulation parameters. Added "particle" struct.