Abstract
This paper displays numerical simulation for bright and dark optical solitons that emerge from Fokas-Lenells equation which is studied in the context of dispersive solitons in polarization-preserving fibers. The Laplace-Adomian decomposition scheme is the numerical tool adopted in the paper. The numerical results, for bright and dark solitons, are expository and therefore supplement the analytical developments, thus far.