diff --git a/lib/objects.py b/lib/objects.py index 0fcbc90..818df52 100755 --- a/lib/objects.py +++ b/lib/objects.py @@ -129,7 +129,7 @@ class System(Body): return mu @property - def ecc(self): #exentricity of sub system of + def ecc(self): #exentricity of two body sub system if len(self.bodylist) == 2 : ecc = (2.*self.E*(np.linalg.norm(self.L)**2))/((G**2)*(self.M**2)*(self.mu**3)) + 1. else : @@ -137,7 +137,7 @@ class System(Body): return ecc @property - def sma(self): #semi major axis of system (if composed of 2 bodies) + def sma(self): #semi major axis of two body sub system if len(self.bodylist) == 2 : sma = -G*self.M*self.mu/(2.*self.E) else : diff --git a/plots/2bodies_leapfrog_a_e.png b/plots/2bodies_leapfrog_a_e.png new file mode 100644 index 0000000..d477a79 Binary files /dev/null and b/plots/2bodies_leapfrog_a_e.png differ diff --git a/plots/2bodies_leapfrog_dEm.png b/plots/2bodies_leapfrog_dEm.png index 61428be..309586d 100644 Binary files a/plots/2bodies_leapfrog_dEm.png and b/plots/2bodies_leapfrog_dEm.png differ diff --git a/plots/2bodies_leapfrog_dL2.png b/plots/2bodies_leapfrog_dL2.png index d7f1abc..e9e06b7 100644 Binary files a/plots/2bodies_leapfrog_dL2.png and b/plots/2bodies_leapfrog_dL2.png differ