This manuscript adds a divorce term to a divorce-marriage model that has a fear effect. With MATLAB simulations (ode45), the entire time-dependent system of ordinary differential equations is solved numerically. Line and surface plots with varying sets of parameter values illustrate population variation over time. Equilibrium analysis is done for different cases (such as populations with no divorce, female-only, and male-only populations), and the stability analysis is made using Jacobian eigenvalues. The innovation of this study lies in the explicit incorporation of a divorce compartment into the developed marriage–divorce–fear model, providing an improved representation of social dynamics. Our expansion offers a deeper comprehension of the interplay between fear, divorce, and remarriage in shap- ing long-term population dynamics, unlike previous models that examined fear effects without directly incorporating divorce. This contribution enhances the mathematical modeling of social systems by elucidating novel bifurcation thresholds and stability requirements that encapsulate pivotal transitions in marital dynamics.