Basins of convergence in the circular Sitnikov four-body problem with non-spherical primaries

The Newton-Raphson basins of convergence, related to the equilibrium points, in the Sitnikov four-body problem with non-spherical primaries are numerically investigated. We monitor the parametric evolution of the positions of the roots, as a function of the oblateness coefficient. The classical Newton-Raphson optimal method is used for revealing the basins of convergence, by classifying dense grids of initial conditions in several types of two-dimensional planes. We perform a systematic and thorough analysis in an attempt to understand how the oblateness coefficient affects the geometry as well as the basin entropy of the convergence regions. The convergence areas are related with the required number of iterations and also with the corresponding probability distributions.