diff --git a/jobs/src/Galaxy_Collision_Simulation/3D/simulate.py b/jobs/src/Galaxy_Collision_Simulation/3D/simulate.py
new file mode 100644
index 0000000000000000000000000000000000000000..0d2dd5cafc41e67eadf2f255c3eb165a67fbbc7a
--- /dev/null
+++ b/jobs/src/Galaxy_Collision_Simulation/3D/simulate.py
@@ -0,0 +1,12 @@
+"""
+This code is used to run the simulation of the galaxy collision.
+Adjust the number of Particles (num_particles), the number of time steps (num_steps),
+the time step (time_step), the number of snapshots to display and if the axes of 
+the plots shoulb be fixed or not (fix_axes), and then run the simulation
+"""
+
+from galaxysimulation import GalaxyCollisionSimulation
+
+sim = GalaxyCollisionSimulation(num_particles=2)
+sim.simulate(num_steps=10000, time_step=0.001)
+sim.display_snapshots(200, fix_axes=True)
\ No newline at end of file