The present paper deals with the periodic orbits generated by Lagrangian solutions of the restricted three-body problem when both the primaries are oblate bodies. We have illustrated the periodic orbits for different values of $\mu, h,\sigma_1$ and $\sigma_2$ ($h$ is energy constant, $\mu$ mass ratio of the two primaries, $\sigma_1$ and $\sigma_2$ are oblateness factors). These orbits have been determined by giving displacements along the tangent and normal to the mobile coordinates as defined by Karimov and Sokolsky \cite{Kari}. We have applied the predictor-corrector algorithm to construct the periodic orbits in an attempt to unveil the effect of oblateness of the primaries by taking the fixed values of parameters $\mu, h, \sigma_1$ and $\sigma_2$.

In the framework of restricted five bodies problem, the existence and stability of the libration points are explored and analysed numerically, under the effect of non--isotropic mass variation of the fifth body (test particle or infinitesimal body). The evolution of the positions of these points and the possible regions of motion are illustrated, as a function of the perturbation parameter. We perform a systematic investigation in an attempt to understand how the perturbation parameter due to variable mass of the fifth body, affects the positions, movement and stability of the libration points. In addition, we have revealed how the domain of the basins of convergence associated with the libration points are substantially influenced by the perturbation parameter.

The analysis of the affect of angular velocity on the geometry of the basins of convergence (BoC) linked to the equilibrium points in the restricted three-body problem is illustrated when the primaries are source of radiation. The bivariate scheme of the Newton-Raphson (N-R) iterative method has been used to discuss the topology of the basins of convergence. The parametric evolution of the fractality of the convergence plane is also presented where the degree of fractality is illustrated by evaluating the basin entropy of the convergence plane.

In the present work a systematic study has been presented in the context of the existence of libration points, their linear stability, the regions of motion where the third particle can orbit and the domain of basins of convergence linked to libration points in the spatial configuration of the collinear restricted four-body problem with non-spherical primaries (i.e., the primaries are oblate or prolate spheroid). The parametric evolution of the positions of the libration points as function of the oblateness and prolateness parameters of the primaries and the stability of these points in linear sense are illustrated numerically. Moreover, the numerical investigation shows that the only libration points which lie on either of the axes are linearly stable for several combinations of the oblateness parameter and mass parameter whereas the non-collinear libration points are found linearly unstable, consequently unstable in nonlinear sense also, for studied value of mass parameter and oblateness parameter. Moreover, the regions of possible motion are also depicted, where the infinitesimal mass is free to orbit, as function of Jacobian constant. In addition, the basins of convergence (BoC) linked to the libration points are illustrated by using the multivariate version of the Newton-Raphson (NR) iterative scheme.

In the framework of photogravitational version of the restricted five-body problem, the existence and stability of the in-plane equilibrium points, the possible regions for motion are explored and analysed numerically, under the combined effect of small perturbations in the Coriolis and centrifugal forces. Moreover, the multivariate version of the Newton-Raphson iterative scheme is applied in an attempt to unveil the topology of the basins of convergence linked with the libration points as function of radiation parameters, and the parameters corresponding to Coriolis and centrifugal forces.