A system of three non-linear differential equations with exponential feedback terms is proposed to model the self-regulating cortisol secretion system and explain the fluctuation patterns observed in clinical data. It is shown that the model exhibits bifurcation and chaos patterns for a certain range of parametric values. This helps us to explain clinical observations and characterize different dynamic behaviors of the self-regulative system.