Tibeuleu/KozaiLidov
1
0
Fork 0
Code Activity

modify output name

Browse Source
This commit is contained in:
Thibault Barnouin
2021-11-06 02:50:19 +01:00
parent e5a5601771
commit ca9ca60b92
4 changed files with 10 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
main.py
View File

@@ -16,13 +16,13 @@ def main():
e = np.array([0., 0., 1./4.]) # Eccentricity
psi = np.array([0., 0., 80.])*np.pi/180. # Inclination of the orbital plane in degrees
x1 = np.array([-1., 0., 0.])*a[0]
x2 = np.array([1., 0., 0.])*a[1]
x1 = np.array([0., -1., 0.])*a[0]
x2 = np.array([0., 1., 0.])*a[1]
x3 = np.array([np.cos(psi[2]), 0., np.sin(psi[2])])*a[2]
q = np.array([x1, x2, x3])
v1 = np.array([0., -np.sqrt(G*m[1]**2/((m[0]+m[1])*np.sqrt(np.sum((q[0]-q[1])**2)))), 0.])
v2 = np.array([0., np.sqrt(G*m[0]**2/((m[0]+m[1])*np.sqrt(np.sum((q[0]-q[1])**2)))), 0.])
v1 = np.array([-np.sqrt(G*m[1]**2/((m[0]+m[1])*np.sqrt(np.sum((q[0]-q[1])**2)))), 0., 0.])
v2 = np.array([np.sqrt(G*m[0]**2/((m[0]+m[1])*np.sqrt(np.sum((q[0]-q[1])**2)))), 0., 0.])
v3 = np.array([0., np.sqrt(G*(m[0]+m[1])*(2./np.sqrt(np.sum(q[2]**2))-1./a[2])), 0.])
v = np.array([v1, v2, v3])