Motion of Spinning Particles around Black Holes

The motion of spinning particles around compact objects, for example a rotating stellar object moving around a supermassive black hole, is described by differential equations that are, in general, non-integrable. In this work, we present a computational code that integrates the equations of motion of a spinning particle moving in the equatorial plane of a spherically symmetric spacetime and gives a visual representation of its trajectory. This code is open source, freely available and modular, so that users may extend its application not only to the Schwarzschild metric, but also to other backgrounds.
14 pages, 2 figures