“…It should be observed that the perturbed coupled nonlinear Schrödinger differential equations systems, which describe wave propagation in real optical media, do not present analytical solution (Korpinar, 2019). In literature there are several numerical approaches, whose objective is to describe the propagation of perturbed solitons in dielectric environments, most of them using the finite difference method (Cirilo, Natti, Romeiro, Natti, & Oliveira, 2010;Mahboub & Bekhti, 2014;Li, Li, Cao, & Li, 2017) or the finite element method (Queiroz, Natti, Romeiro, & Natti, 2006;Saka, 2012;Karczewska, Rozmej, Szczecinski, & Boguniewicz, 2016). On the other hand, to solve numerically the system of generated equations, authors use various methods like Newton's method (Syafwan, Susanto, Coix, & Malomed, 2012), Crank-Nicolson method (Liang, Khaliq, & Sheng, 2014), Runge-Kutta Method (Shirazi, Solaimani, Farnam, Ghalandari, & Aleomraninejad, 2017), among others (Dehghan & Taleei, 2010;Gupta, Dalal, & Mishra, 2014;Gupta, Dalal, & Mishra, 2015;Choudlhary, Jha, Mishra, & Vandana, 2018).…”