diff --git a/main.py b/main.py
index f63e80b..5858c22 100755
--- a/main.py
+++ b/main.py
@@ -9,7 +9,7 @@ from lib.units import *
 
 def main():
     #initialisation
-    m = np.array([1., 1., 0])*Ms  # Masses in Solar mass
+    m = np.array([1., 1., 1e-5])*Ms  # Masses in Solar mass
     a = np.array([1., 1., 5.])*au   # Semi-major axis in astronomical units
     e = np.array([0., 0., 1./4.])   # Eccentricity
     psi = np.array([0., 0., 0.])*np.pi/180.    # Inclination of the orbital plane in degrees
@@ -31,12 +31,12 @@ def main():
     dyn_syst.COMShift()
 
     duration, step1, step2 = 100*yr, 1e4, 1e5
-    E1, L1 = dyn_syst.leapfrog(duration, step1, recover_param=True)#, display=True)
+    E1, L1 = dyn_syst.leapfrog(duration, step1, recover_param=True, display=True)
     E2, L2 = dyn_syst.leapfrog(duration, step2, recover_param=True)#, display=True)
     #E1, L1 = dyn_syst.hermite(duration, step1, recover_param=True)#, display=True)
     #E2, L2 = dyn_syst.hermite(duration, step2, recover_param=True)#, display=True)
     parameters = [duration, [step1, step2], dyn_syst, "leapfrog"]
-    display_parameters([E1, E2], [L1, L2], parameters=parameters, savename="3bodies_leapfrog")
+    display_parameters([E1, E2], [L1, L2], parameters=parameters, savename="3bodies_mass_leapfrog")
 
     return 0
 
diff --git a/plots/3bodies_mass_leapfrog_dEm.png b/plots/3bodies_mass_leapfrog_dEm.png
new file mode 100644
index 0000000..9c59db8
Binary files /dev/null and b/plots/3bodies_mass_leapfrog_dEm.png differ
diff --git a/plots/3bodies_mass_leapfrog_dL2.png b/plots/3bodies_mass_leapfrog_dL2.png
new file mode 100644
index 0000000..12e140f
Binary files /dev/null and b/plots/3bodies_mass_leapfrog_dL2.png differ
diff --git a/plots/3bodies_mass_leapfrog_dynsyst.gif b/plots/3bodies_mass_leapfrog_dynsyst.gif
new file mode 100644
index 0000000..2e39d2c
Binary files /dev/null and b/plots/3bodies_mass_leapfrog_dynsyst.gif differ