This manuscript presents a numerical model to evaluate the stability of slopes located along the bank (Go Cong) with four different plant species, namely Eucalyptus, Stipa, Artemisia, and Rosmarinus. The analyses were performed using shear strength reduction method combined with finite element method implemented in the commercial Plaxis 2D software. In those analyses, plant-induced additional cohesion, which was a function of root depth and root tensile strength, was initially calculated and subsequently added into the original cohesion of bare soil. The Factor of Safety (FoS) was then calculated to examine the effectiveness of the plant species in improving the slope stability, in various scenarios of slope height and angle. Besides, parametric studies were conducted to investigate the increment of FoS due to vegetation reinforcement for several cases of internal friction angle and cohesion of bare soil.