diff --git a/main.py b/main.py index 738a70e..8683319 100755 --- a/main.py +++ b/main.py @@ -12,9 +12,9 @@ from lib.units import * def main(): #initialisation - m = np.array([1., 1., 1e-1],dtype=np.longdouble)*Ms#/Ms # Masses in Solar mass - a = np.array([1.0, 1.0, 7.0],dtype=np.longdouble)*au#/au # Semi-major axis in astronomical units - e = np.array([0., 0., 0.10],dtype=np.longdouble) # Eccentricity + m = np.array([1., 1., 0.],dtype=np.longdouble)*Ms#/Ms # Masses in Solar mass + a = np.array([1.00, 1.00, 10.0],dtype=np.longdouble)*au#/au # Semi-major axis in astronomical units + e = np.array([0., 0., 0.25],dtype=np.longdouble) # Eccentricity psi = np.array([0., 0., 80.],dtype=np.longdouble)*np.pi/180. # Inclination of the orbital plane in degrees x1 = np.array([0., -1., 0.],dtype=np.longdouble)*a[0]*(1.+e[0]) @@ -28,12 +28,12 @@ def main(): v = np.array([v1, v2, v3],dtype=np.longdouble) #integration parameters - duration, step = 5000*yr, np.longdouble(1./1.*86400.) #integration time and step in seconds + duration, step = 5000*yr, np.longdouble(1./2.*86400.) #integration time and step in seconds integrator = "leapfrog" n_bodies = 3 display = False gif = False - savename = "{0:d}bodies_{1:s}".format(n_bodies, integrator) + savename = "{0:d}bodies_massless_{1:s}".format(n_bodies, integrator) display_param = True #simulation start diff --git a/plots/3bodies_1_leapfrog_E.png b/plots/3bodies_1_leapfrog_E.png new file mode 100644 index 0000000..092afdc Binary files /dev/null and b/plots/3bodies_1_leapfrog_E.png differ diff --git a/plots/3bodies_1_leapfrog_L.png b/plots/3bodies_1_leapfrog_L.png new file mode 100644 index 0000000..cae1a6b Binary files /dev/null and b/plots/3bodies_1_leapfrog_L.png differ diff --git a/plots/3bodies_1_leapfrog_a_e.png b/plots/3bodies_1_leapfrog_a_e.png new file mode 100644 index 0000000..793b580 Binary files /dev/null and b/plots/3bodies_1_leapfrog_a_e.png differ diff --git a/plots/3bodies_1_leapfrog_dEm.png b/plots/3bodies_1_leapfrog_dEm.png new file mode 100644 index 0000000..1f2a156 Binary files /dev/null and b/plots/3bodies_1_leapfrog_dEm.png differ diff --git a/plots/3bodies_1_leapfrog_dL2.png b/plots/3bodies_1_leapfrog_dL2.png new file mode 100644 index 0000000..13ea865 Binary files /dev/null and b/plots/3bodies_1_leapfrog_dL2.png differ diff --git a/plots/3bodies_1_leapfrog_phi.png b/plots/3bodies_1_leapfrog_phi.png new file mode 100644 index 0000000..003b266 Binary files /dev/null and b/plots/3bodies_1_leapfrog_phi.png differ diff --git a/plots/3bodies_massless_leapfrog_E.png b/plots/3bodies_massless_leapfrog_E.png new file mode 100644 index 0000000..00e602e Binary files /dev/null and b/plots/3bodies_massless_leapfrog_E.png differ diff --git a/plots/3bodies_massless_leapfrog_L.png b/plots/3bodies_massless_leapfrog_L.png new file mode 100644 index 0000000..93192fa Binary files /dev/null and b/plots/3bodies_massless_leapfrog_L.png differ diff --git a/plots/3bodies_massless_leapfrog_a_e.png b/plots/3bodies_massless_leapfrog_a_e.png new file mode 100644 index 0000000..6ba65aa Binary files /dev/null and b/plots/3bodies_massless_leapfrog_a_e.png differ diff --git a/plots/3bodies_massless_leapfrog_dEm.png b/plots/3bodies_massless_leapfrog_dEm.png new file mode 100644 index 0000000..2e4c0a4 Binary files /dev/null and b/plots/3bodies_massless_leapfrog_dEm.png differ diff --git a/plots/3bodies_massless_leapfrog_dL2.png b/plots/3bodies_massless_leapfrog_dL2.png new file mode 100644 index 0000000..956b8ba Binary files /dev/null and b/plots/3bodies_massless_leapfrog_dL2.png differ diff --git a/plots/3bodies_massless_leapfrog_phi.png b/plots/3bodies_massless_leapfrog_phi.png new file mode 100644 index 0000000..51d3fe7 Binary files /dev/null and b/plots/3bodies_massless_leapfrog_phi.png differ