diff --git a/examples/symba_chambers_2013/init_cond.py b/examples/symba_chambers_2013/init_cond.py
index 6a6dd8526..d02830a11 100755
--- a/examples/symba_chambers_2013/init_cond.py
+++ b/examples/symba_chambers_2013/init_cond.py
@@ -60,10 +60,10 @@
 # Define the initial orbital elements of the big and small bodies
 avalb = default_rng().uniform(0.3, 2.0, Nb)
 avals = default_rng().uniform(0.3, 2.0, Ns)
-evalb = default_rng().uniform(0.0, 0.01, Nb)
-evals = default_rng().uniform(0.0, 0.01, Ns)
-incvalb = default_rng().uniform(0.0, 0.005 * 180 / np.pi, Nb)
-incvals = default_rng().uniform(0.0, 0.005 * 180 / np.pi, Ns)
+evalb = np.abs(default_rng().normal(loc=0.0, scale=0.13, size=Nb))
+evals = np.abs(default_rng().normal(loc=0.0, scale=0.13, size=Ns))
+incvalb = default_rng().normal(loc=0.0, scale=1.0, size=Nb)
+incvals = default_rng().normal(loc=0.0, scale=1.0, size=Ns)
 capomvalb = default_rng().uniform(0.0, 360.0, Nb)
 capomvals = default_rng().uniform(0.0, 360.0, Ns)
 omegavalb = default_rng().uniform(0.0, 360.0, Nb)