This final project discusses the mathematical model of love and happiness with affair. In this model, system of differential equations is used to describe quality of each individual. There are four conditions observed, in which each case is analyzed using Routh-Hurwitz criterion. Furthermore, simulation is given to show the behaviour and the stability of each case. Keywords: Differential e…