Abstract
This research paper investigates the numerical solutions of the predator-prey model through five recent numerical schemes (Adomian decomposition, El Kalla, cubic B-spline, extended cubic B-spline, exponential cubic B-spline). We investigate the obtained computational solutions via the modified Khater methods. This model is considered as a well-known bimathematical model to describe the prototype of an excitable system. The obtained solitary solutions emerge the localized wave packet as a persistent and dominant feature. The accuracy of the obtained numerical solutions is investigated by calculating the absolute error between the exact and numerical solutions. Many sketches are given to illustrate the matching between the exact and numerical solutions.