The perturbed restricted four-body problem has been investigated numerically when the primary is oblate/ prolate spheroid and source of radiation as well. The movement of the positions of the libration points and their stability are unveiled numerically. The regions of possible motion and the basins of convergence associated with the libration points are illustrated by applying the well known multivariate version of the Newton-Raphson iterative scheme.