diff --git a/examples/symba_chambers_2013/init_cond.py b/examples/symba_chambers_2013/init_cond.py
index 5f6782491..38b023ed1 100755
--- a/examples/symba_chambers_2013/init_cond.py
+++ b/examples/symba_chambers_2013/init_cond.py
@@ -54,6 +54,8 @@
 dens = 3000.0 / (MU2KG / DU2M**3)
 Rb = (3 * Mb / (4 * np.pi * dens) )**(1.0 / 3.0)
 Rs = (3 * Ms / (4 * np.pi * dens) )**(1.0 / 3.0)
+MTINY = 1e-2 * Ms
+sim.param['GMTINY'] = MTINY * GU
 
 # Define the initial orbital elements of the big and small bodies
 avalb = default_rng().uniform(0.3, 2.0, Nb)