In this study, a novel method known as the Haar wavelet collocation method (HWCM) is used to analyze the mathematical model of Chlamydia transmission in the United States. We use dependent variables with various parameters, such as birth rate, mortality rate, and recovery rate, etc., to explore the nature of the Chlamydia disease in susceptible unvaccinated individuals, susceptible vaccinated individuals, infected, treated, and recovered individuals. The ordinary differential equations (ODEs) in this model are coupled nonlinear. In this method, the Chlamydia model is converted into a system of non-linear algebraic equations using properties of the operational matrix of Haar wavelets. Later, the Newton–Raphson approach is used to extract the unknown coefficients. Mathematica software has been used for all calculations. Tables and graphs contain tabulated results of calculations. As a result, it can be seen that the current approach is more precise than those used in the literature. Also, we discuss the effect of different parameters on susceptible unvaccinated individuals, susceptible vaccinated individuals, infected, treated, and recovered individuals through graphs.