This code is to repeat the Figure 2 in the paper "Unlimited multistability in multisite phosphorylation systems (https://www.nature.com/articles/nature08102). The simulation was performed using Spyder for Tellurium 3.7. Please note that certain intitial conditions are necessary in order to repeat the exact results in the Figure 2b.
If you use any of the codes, please cite the article (https://doi.org/10.1371/journal.pone.0314875). This work has been submitted to BioModels (https://www.ebi.ac.uk/biomodels/MODEL2002110001#Overview).