mass 1.989e30[kg] "Star mass"
gravity G_const "Gravitation constant"
au 149597870700[m] "Astronomical unit"
distance 30*au "Initial distance between stars"
tilt -14.06884013142[deg] "Tilt angle for figure-8 configuration"
x1 distance*cos(tilt) "Star 1 initial x-coordinate"
y1 distance*sin(tilt) "Star 1 initial y-coordinate"
x2 -x1 "Star 2 initial x-coordinate"
y2 -y1 "Star 2 initial y-coordinate"
x3 0 "Star 3 initial x-coordinate"
y3 0 "Star 3 initial y-coordinate"
u1 -u3/2 "Star 1 initial velocity, x-component"
v1 -v3/2 "Star 1 initial velocity, y-component"
u2 -u3/2 "Star 2 initial velocity, x-component"
v2 -v3/2 "Star 2 initial velocity, y-component"
u3 -0.93240737*sqrt(mass*gravity/distance) "Star 3 initial velocity, x-component"
v3 -0.86473146*sqrt(mass*gravity/distance) "Star 3 initial velocity, y-component"
T1 6.32591398*sqrt(distance^3/(mass*gravity)) "Orbital period, figure 8 configuration"
x4 distance*cos(0[deg]) "Star 4 initial x-coordinate"
y4 distance*sin(0[deg]) "Star 4 initial y-coordinate"
x5 distance*cos(120[deg]) "Star 5 initial x-coordinate"
y5 distance*sin(120[deg]) "Star 5 initial y-coordinate"
x6 distance*cos(240[deg]) "Star 6 initial x-coordinate"
y6 distance*sin(240[deg]) "Star 6 initial y-coordinate"
V_Lagrange sqrt(gravity*mass/distance)/(3^0.25) "Initial speed, Lagrange configuration"
u4 -V_Lagrange*sin(0[deg]) "Star 4 initial velocity, x-component"
v4 V_Lagrange*cos(0[deg]) "Star 4 initial velocity, y-component"
u5 -V_Lagrange*sin(120[deg]) "Star 5 initial velocity, x-component"
v5 V_Lagrange*cos(120[deg]) "Star 5 initial velocity, y-component"
u6 -V_Lagrange*sin(240[deg]) "Star 6 initial velocity, x-component"
v6 V_Lagrange*cos(240[deg]) "Star 6 initial velocity, y-component"
T2 2*pi*distance/V_Lagrange "Orbital period, Lagrange configuration"