ClearAll["Global`*"]
G=6.6743*
-11
10
;​​c=2.99792458*
8
10
;​​Period=5.599829;​​P=Period*86400;​​ecc=0.018;​​a=
3
2
P
*G*(m1+m2)
4*
2
π
​​μ=
m1m2
m1+m2
;​​aph=a(1+ecc);​​x0=
m2*aph
m1+m2
​​y0=0;​​X0=-
m1*aph
m1+m2
​​Y0=0;​​Vel1=
1/2
Gμ
a
m2
m1
1-ecc
1+ecc
(*initialvelocitystar1*)​​Vel2=-
1/2
Gμ
a
m2
m1
1-ecc
1+ecc
m1
m2
(*initialvelocitystar2*)​​​​​​system1=​​x1''[t]-
Gm1m2(x1[t]-x2[t])
m1
3/2
(
2
(x1[t]-x2[t])
+
2
(y1[t]-y2[t])
)
,​​y1''[t]-
Gm1m2(y1[t]-y2[t])
m1
3/2
(
2
(x1[t]-x2[t])
+
2
(y1[t]-y2[t])
)
,​​x2''[t]
Gm1m2(x1[t]-x2[t])
m2
3/2
(
2
(x1[t]-x2[t])
+
2
(y1[t]-y2[t])
)
,​​y2''[t]
Gm1m2(y1[t]-y2[t])
m2
3/2
(
2
(x1[t]-x2[t])
+
2
(y1[t]-y2[t])
)
;​​​​initials={​​x1[0]x0,​​x1'[0]0,​​y1[0]0,​​y1'[0]Vel1,​​x2[0]X0,​​x2'[0]0,​​y2[0]0,​​y2'[0]Vel2};​​​​{xbh1,ybh1,xbh2,ybh2}=DSolveValue[{system1,initials},{x1,y1,x2,y2},t];​​bh1={(xbh1[t*86400]/(1.5*10^11)),(ybh1[t*86400]/(1.5*10^11))};(*largerstarpositionvectors*)​​bh2={(xbh2[t*86400]/(1.5*10^11)),(ybh2[t*86400]/(1.5*10^11))};(*smallerstarpositionvectors*)
0.734189
1/3
(m1+m2)
0.747404m2
2/3
(m1+m2)
-
0.747404m1
2/3
(m1+m2)
9.36442×
-6
10
2
m2
4/3
(m1+m2)
-
9.36442×
-6
10
m1
2
m2
4/3
(m1+m2)
m2
Set
:Lists {xbh1,ybh1,xbh2,ybh2} and DSolveValue
′′
x1
[t]-
6.6743×
-11
10
m2(x1[t]-x2[1])
3/2
(Power[2]+Power[2])
,
′′
y1
[t]-
6.6743×
-11
10
m2(y1[t]-y2[1])
3/2
(Power[2]+Power[2])
,
′′
x2
[t]
6.6743×
-11
10
m1(x1[t]-x2[1])
3/2
(Power[2]+Power[2])
,
′′
y2
[t]
6.6743×
-11
10
m1(y1[t]-y2[1])
3/2
(Power[2]+Power[2])
,x1[0]
0.747404m2
2/3
(m1+m2)
,
′
x1
[0]0,y1[0]0,
′
y1
[0]9.36442×
-6
10
Power[2]Power[2]
,x2[0]-
0.747404m1
2/3
(m1+m2)
,
′
x2
[0]0,y2[0]0,
′
y2
[0]-
9.36442×
-6
10
m1
Power[2]Power[2]
m2
,{x1,y1,x2,y2},t are not the same shape.