Reproduced plots of Figure 2A in the original publication. Variable values used are as follows:
r = 0.03
a = 0.02 
mu = 47
All other values are as indicated in the paper.

Model simulated and plot produced using COPASI 4.24 (Build 197).