Work by Roberto Molina Ballesteros and Sofia Llàcer Caro for the couse Deterministic Modelling in the MSc Modelling for Science and Engineering at the Universitat Autònoma de Barcelona.
Code under Scripts can be run by foowing these insctructions as incorporated in the files:
# States: uncomment the desired one
states=[0,1]
#states = [1,2,3]
#states = list(range(20))