In this paper, we present a dynamic model of syphilis. The goal of this work is to provide a mathematical model for syphilis disease that contains information about the suspect, infected, recovered, and exposed. Boundedness, positivity, and unique solutions at equilibrium points are used in the qualitative and quantitative study of a proposed model using the power law kernel. The Routh-Hurwitz stability criteria are used to address the model's local and global stability. Chaos control will use the regulate for linear responses approach to bring the system to stabilize according to its points of equilibrium. Using Lyapunov function tests, the global stability of free and endemic equilibrium points is confirmed. The Lipschitz condition and fixed point theory are utilized to satisfy the requirements for the existence and uniqueness of the exact solution. The generalized form of the power law kernel is solved using a two-step Lagrange polynomial method to investigate the impact of the fractional operator using numerical simulations that illustrate the effects of various parameters on the illness.