This paper is a numerical simulation of highly dispersive optical solitons in birefringent fibers with polynomial nonlinear form, which achieved for the first time. The algorithmic approach applied usage Laplace–Adomian decomposition scheme. Dark and bright soliton simulations are presented. error measure has very low count, thus, almost an exact replica such that analytically arise from govern...