Basins of attraction of equilibrium points in the planar circular restricted five-body problem

We numerically explore the Newton-Raphson basins of convergence, related to the libration points (which act as attractors), in the planar circular restricted five-body problem (CR5BP). The evolution of the position and the linear stability of the equilibrium points is determined, as a function of the value of the mass parameter. The attracting regions, on several types of two dimensional planes, are revealed by using the multivariate version of the classical Newton-Raphson iterative method. We perform a systematic investigation in an attempt to understand how the mass parameter affects the geometry as well as the degree of fractality of the basins of attraction. The regions of convergence are also related with the required number of iterations and also with the corresponding probability distributions.