Simulation of 1D self gravitating object

Phase space evolution of self gravitation system of particles in one dimensions (Total number of particles: N=128)

Initial Conditions :

Positions: a= [0,1,2,3,4......., 127] (in x-axis).

Momentum: P_a = sin(2πa/N); N=128

Phase space evolution of self gravitation system of particles in one dimensions (Total number of particles: N=128)

Initial Conditions :

Positions: a= [0,1,2,3,4......., 127] (in x-axis).

Momentum: sin(θ), where θ = [−π/2, π/2] over x=[0,127]

Evolution of Momentum of system of particles

Evolution of positions of system of particles