dx1/dt = -x1.x3 + k12.x2 + EGP0(1 - x5) - F - Fr + u1
dx2/dt = x1.x3 - k12.x2 - x2.x4
dx3/dt = -ka1.x3 + kb1.x6
dx4/dt = -ka1.x4 + kb1.x6
dx5/dt = -ka1.x5 + kb1.x6
dx6/dt = -ke.x6 + ki.x7
dx7/dt = -ki.x7 + ki.x8
dx8/dt = -ki.x8 + u2
states: [x1 x2 x3 x4 x5 x6 x7 x8]
parameters: [k12 EGP0 ka1 ka2 ka3 kb1 kb2 kb3 ki ke]
dis-continuous functions: [1] Fr = 0.003(x1-9), when x1>=9
Fr = 0, [otherwise]
[2] F = 0.97 , when x1>=4.5
F = 0.97.x1/4.5 [otherwise]