CNLN: Phosphorelay

From:
Klipp, E., Herwig, R., Kowald, A., Wierling, C., Lehrach, H.: Systems biology in practice: concepts, implementation and application.
Dreossi T., Dang T., Piazza C.: Reachability Computation for Polynomial Dynamical Systems.

Dynamics:
x1' = -0.4*x1 + 5*x3*x4
x2' = 0.4*x1 - x2
x3' = x2-5*x3*x4
x4' = 5*x5*x6 - 5*x3*x4
x5' = -5*x5*x6 + 5*x3*x4
x6' = 0.5*x7 - 5*x5*x6
x7' = -0.5*x7 + 5*x5*x6

Parameters:
Delta = 0.01 (for Euclid discretization)

Initial set:
x1(0), x2(0), x3(0), x4(0), x5(0), x6(0), x7(0) \in [1.00, 1.01]