Abstract
In this paper, we investigate the (2+1)-dimensional Kadomtsev-Petviashvili-Benjamin-Bona Mahony equation using two effective methods: the unified scheme and the advanced auxiliary equation scheme, aiming to derive precise wave solutions. These solutions are expressed as combinations of trigonometric, rational, hyperbolic, and exponential functions. Visual representations, including three-dimensional (3D) and two-dimensional (2D) combined charts, are provided for some of these solutions. The influence of the nonlinear parameter p on the wave type is thoroughly examined through diverse figures, illustrating the profound impact of nonlinearity. Additionally, we briefly investigate the Hamiltonian function and the stability of the model using a planar dynamical system approach. This involves examining trajectories, isoclines, and nullclines to illustrate stable solution paths for the wave variables. Numerical results demonstrate that these methods are reliable, straightforward, and potent tools for analyzing various nonlinear evolution equations found in physics, applied mathematics, and engineering.