Abstract
This paper is a numerical simulation of highly dispersive optical solitons in birefringent fibers with polynomial nonlinear form, which is achieved for the first time. The algorithmic approach is applied with the usage of the Laplace–Adomian decomposition scheme. Dark and bright soliton simulations are presented. The error measure has a very low count, and thus, the simulations are almost an exact replica of such solitons that analytically arise from the governing system. The suggested iterative scheme finds the solution without any discretization, linearization, or restrictive assumptions.