In this paper, an analytical method based on the Bernoulli differential equation for extracting new complex soliton solutions to the Gilson–Pickering model is applied. A set of new complex soliton solutions to the Gilson–Pickering model are successfully constructed. In addition, 2D and 3D graphs and contour simulations to the complex soliton solutions are plotted with the help of computational programs. Finally, at the end of the manuscript a conclusion about new complex soliton solutions is given.