The model has been encoded in COPASI 4.24 (Build 197) and the figure 1 lower panel has been reproduced using COPASI The value for c was used to be 0.02 which produces the exact figure