The analytical and numerical solutions of the (2+1) dimensional, Fisher-Kolmogorov-Petrovskii-Piskunov ((2+1) D-Fisher-KPP) model are investigated by employing the modified direct algebraic (MDA), modified Kudryashov (MKud.), and trigonometric-quantic B-spline (TQBS) schemes. This model, which arises in population genetics and nematic liquid crystals, describes the reaction–diffusion system by traveling waves in population genetics and the propagation of domain walls, pattern formation in bi-stable systems, and nematic liquid crystals. Many novel analytical solutions are constructed. These solutions are used to evaluate the requested numerical technique’s conditions. The numerical solutions of the considered model are studied, and the absolute value of error between analytical and numerical is calculated to demonstrate the matching between both solutions. Some figures are represented to explain the obtained analytical solutions and the match between analytical and numerical results. The used schemes’ performance shows their effectiveness and power and their ability to handle many nonlinear evolution equations.