Reproduction of fig 3b from the article Integration performed using SBML ODESolver (20080507) Plot created using Xmgrace