Estos códigos se dessrrollaron usando el lenguaje fortran, en junio del 2019, para reproducir los resultados del articulo de investigación Sprott, J. C. (2005). Dynamical models of happiness. Nonlinear Dynamics, Psychology, and Life Sciences, 9(1), 23-36.
Se seleccionan ecuaciones diferenciales de interes, para resolverlas se usa el método de runge-Kutta de orden 2, ya que la precisión cuantitativa no era necesaria. Las ecuaciones diferenciales se transforman en un sistema acoplado de ecuaciones lineales y es este el que se resuelve.