PS Numerical Estimate of the Precession in Quasar OJ287

With such a large precession of 39 degrees per orbit, many orders of magnitude bigger than in the solar system, it should be possible to adjust conditions to use the experimentally observed mass M to find the precession numerically, and to measure it graphically. The new equation of motion is:

r bold double dot = (gamma MG / r cubed) ( v bold (v bold dot r bold) / c squared – r bold)

so only M is used. Here gamma is the Lorentz factor. It is not clear how the astronomers estimated M. The above equation can be used to find M by matching the experimentally observed precession of 39 degrees per orbit to the theoretical result. The equation is not soluble analytically but this is no problem because the numerical solution is very precise.


%d bloggers like this: