From 81ac7ce307b0a931e17b1a7b40cce3f4bc714e54 Mon Sep 17 00:00:00 2001
From: Thibault Barnouin <thibault.barnouin2@etu.unistra.fr>
Date: Tue, 18 Jan 2022 22:34:32 +0100
Subject: [PATCH] final main

---
 main.py | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/main.py b/main.py
index d2c02d1..db9ef05 100755
--- a/main.py
+++ b/main.py
@@ -18,7 +18,7 @@ def main():
             dtype=np.longdouble)*au # Semi-major axis in astronomical units
     e = np.array([0., 0., 0.10],
             dtype=np.longdouble)    # Eccentricity
-    psi = np.array([0., 0., 35.],
+    psi = np.array([0., 0., 80.],
             dtype=np.longdouble)*np.pi/180. # Inclination of the orbital plane in degrees
 
     x1 = a[0]*(1.+e[0])*np.array([0., -1., 0.],
@@ -38,14 +38,14 @@ def main():
     v = np.array([v1, v2, v3],dtype=np.longdouble)
 
     #integration parameters
-    duration = 2000*yr    #integration time in seconds
+    duration = 5000*yr    #integration time in seconds
     step = np.longdouble(1.0/1.*86400.) #integration step in seconds
     integrator = "leapfrog"
     n_bodies = 3
     display = False
     gif = False
     blackstyle = True
-    savename = "{0:d}bodies_test_{1:s}".format(n_bodies, integrator)
+    savename = "{0:d}bodies_{1:s}".format(n_bodies, integrator)
     display_param = True
 
     #simulation start