Loading [MathJax]/jax/output/HTML-CSS/jax.js

Sunday, July 11, 2021

I love it when the planets come together

The Xiddlerian astronomers have identified three non-Keplerian planets that circularly orbit a neighboring star. Planet A is three astronomical units away from its star and completes its orbit in three years. Planet B is four astronomical units away from the star and completes its orbit in four years. Finally, Planet C is five astronomical units away from the star and completes its orbit in five years.

They report their findings to Xiddler’s Grand Minister, along with the auspicious news that all three planets are currently lined up with their star. However, the Grand Minister is far more interested in the three planets than the star and wants to know how long it will be until the planets are next aligned. How many years will it be until the three planets are again collinear (not necessarily including the star)?

Let's denote the positions of planets A, B and C at time t (in years) as A(t)=(3cos2πt3,3sin2πt3),B(t)=(4cosπt2,4sinπt2),C(t)=(5cos2πt5,5sin2πt5), respectively. In order to test for collinearity of A(t), B(t) and C(t) we will want to test whether the area of the triangle with ABC has zero area. Given two vectors u,vR2 the area of the parallelogram spanned by u and v is given by the norm of their cross product u×v therefore we can calculate the area of the triangle spanned by A(t), B(t) and C(t) is given by S(t)=12(B(t)A(t))×(C(t)A(t))=12abs(|4cosπt23cos2πt34sinπt23sin2πt35cos2πt53cos2πt35sin2πt53sin2πt3|)=12|(4cosπt23cos2πt3)(5sin2πt53sin2πt3)(5cos2πt53cos2πt3)(4sinπt23sin2πt3)|=12|20cosπt2sin2πt512cosπt2sin2πt315cos2πt3sin2πt5+9cos2πt3sin2πt320cos2πt5sinπt2+12cos2πt3sinπt2+15cos2πt5sin2πt39cos2πt3sin2πt3|=12|20sinπt1012sinπt6+15sin4πt15|

Obviously in the context of area, we would want to keep the absolute value, though since we are trying to find zeros of this trigonometric function, it might be easier to just deal with the signed version ˜S(t)=12(15sin4πt1520sinπt1012sinπt6). The signed area is shown in the figure below for the first full period (until t=60 years).

We see from eyeballing it, as all good Xiddlerian astronomers are wont to do, that the first zero of ˜S is between t=5 and t=10 years, with ˜S(5)=19.4952 and ˜S(10)=11.6913. From here we can resort to whatever friendly neighborhood root solving technique floats our Xiddlerian boats to deduce that the next alignments of Planets A, B and C will occur in ˆt=7.76678 years. For instance, the secant method gets with initial a1=5,a0=10 and an=an2˜S(an1)an1˜S(an2)˜S(an1)˜S(an2) with stopping criterion |˜S(an)|<106 solution quoted above within 6 steps. The figure below shows the alignment of the planets with respect to their shared sun at the center of their concentric orbits.

No comments:

Post a Comment