The nonlinear Schrödinger equation (NLSE) and its various forms have significant applications in the field of soliton theory. The Fokas-Lenells (FL) equation stands as a cornerstone in deepening our understanding of nonlinear wave dynamics within optical systems, particularly concerning the behavior of ultrashort pulses across different media. Its significance lies in providing a comprehensive framework to study and analyze complex phenomena, ultimately contributing to advancements in optical technology and applications. The FL equation is an integrable extension of the NLSE that provides a description of the nonlinear propagation of pulses in optical fiber. This paper seeks to discover optical soliton solutions for the FL equation by employing a modified sub-equation method. Additionally, the sensitivity analysis is described by using the various initial conditions. The main novelty of this paper lies in conducting a sensitivity analysis of the FL equation by examining the effects of various initial conditions, providing deeper insights into how these conditions influence the behavior of soliton solutions. For the physical behavior of the models, some solutions are graphically shown in 2D, 3D, and contour graphs by assigning specific values to the parameters under the provided situation at each solution. As a result, we discovered several new families of exact traveling wave solutions, such as bright solitons, dark solitons, and combined bright and dark solitons. This research opens numerous avenues for further exploration in the field of nonlinear wave dynamics and optical soliton theory. The discovery of exact soliton solutions for the FL equation through a modified sub-equation method paves the way for deeper investigations for newcomer researchers. The results of this study will contribute further to the field of mathematical physics, particularly in enhancing the understanding of nonlinear wave propagation and soliton theory in optical and other physical systems.