
pharmacokinetics

seed_ := 25$

NUMBER OF EQUATIONS$

n_ := 6$

VARIABLES VECTOR$

b_ := {x1,
x2,
x3,
x4,
y1,
y2}$

UNKNOWN PARAMETER(S) VECTOR$

b1_ := {a1,
a2,
b1,
b2,
ka,
kc,
vm}$

RANKING AMONG THE VARIABLES$

bb_ := {x1,
x2,
df(x1,t),
df(x2,t),
df(x1,t,2),
df(x2,t,2),
df(x1,t,3),
df(x2,t,3),
df(x1,t,4),
df(x2,t,4),
x3,
x4,
y1,
y2,
df(x3,t),
df(x4,t),
df(y1,t),
df(y2,t)}$

NUMBER OF INPUT(S)$

nu_ := 0$

NUMBER OF OUTPUT(S)$

ny_ := 2$

NUMBER OF STATE(S) $

nx_ := 4$

MODEL EQUATION(S)$

c_ := {df(x1,t)=( - ((kc + x1)*ka + kc*x3)*(x1 - x2)*a1 - ka*vm*x1)/((kc + x1)*ka + kc*x3),
df(x2,t)=(x1 - x2)*a2,
df(x3,t)=( - ((kc + x1)*ka + kc*x3)*(x3 - x4)*b1 - kc*vm*x3)/((kc + x1)*ka + kc*x3),
df(x4,t)=(x3 - x4)*b2,
y1=x1,
y2=x4}$

CHARACTERISTIC SET$

aa_(1) := df(x2,t) - x1*a2 + x2*a2$

aa_(2) := df(x1,t,3)*df(x1,t)**3*x1**3*ka*vm + 3*df(x1,t,3)*df(x1,t)**2*x1**4*a1*ka*vm - 3*df(x1,t,3)*df(x1,t)**2*x1**3*
x2*a1*ka*vm + 3*df(x1,t,3)*df(x1,t)*x1**5*a1**2*ka*vm - 6*df(x1,t,3)*df(x1,t)*x1**4*x2*a1**2*ka*vm + 3*df(x1,t,3)*df(x1,
t)*x1**3*x2**2*a1**2*ka*vm + df(x1,t,3)*x1**6*a1**3*ka*vm - 3*df(x1,t,3)*x1**5*x2*a1**3*ka*vm + 3*df(x1,t,3)*x1**4*x2**2
*a1**3*ka*vm - df(x1,t,3)*x1**3*x2**3*a1**3*ka*vm - 2*df(x1,t,2)**2*df(x1,t)**2*x1**3*ka*vm - 4*df(x1,t,2)**2*df(x1,t)*
x1**4*a1*ka*vm + 4*df(x1,t,2)**2*df(x1,t)*x1**3*x2*a1*ka*vm - 2*df(x1,t,2)**2*x1**5*a1**2*ka*vm + 4*df(x1,t,2)**2*x1**4*
x2*a1**2*ka*vm - 2*df(x1,t,2)**2*x1**3*x2**2*a1**2*ka*vm + df(x1,t,2)*df(x1,t)**5*x1**2*( - ka + kc) + df(x1,t,2)*df(x1,
t)**5*x1*kc**2 + 5*df(x1,t,2)*df(x1,t)**4*x1**3*a1*( - ka + kc) + 5*df(x1,t,2)*df(x1,t)**4*x1**2*x2*a1*(ka - kc) + df(x1
,t,2)*df(x1,t)**4*x1**2*(5*a1*kc**2 + ka*vm) - 5*df(x1,t,2)*df(x1,t)**4*x1*x2*a1*kc**2 + 10*df(x1,t,2)*df(x1,t)**3*x1**4
*a1**2*( - ka + kc) + 20*df(x1,t,2)*df(x1,t)**3*x1**3*x2*a1**2*(ka - kc) + df(x1,t,2)*df(x1,t)**3*x1**3*(10*a1**2*kc**2 
- a1*ka*vm + b1*ka*vm + b2*ka*vm) + 10*df(x1,t,2)*df(x1,t)**3*x1**2*x2**2*a1**2*( - ka + kc) + 2*df(x1,t,2)*df(x1,t)**3*
x1**2*x2*a1*( - 10*a1*kc**2 - ka*vm) + 10*df(x1,t,2)*df(x1,t)**3*x1*x2**2*a1**2*kc**2 + 10*df(x1,t,2)*df(x1,t)**2*x1**5*
a1**3*( - ka + kc) + 30*df(x1,t,2)*df(x1,t)**2*x1**4*x2*a1**3*(ka - kc) + df(x1,t,2)*df(x1,t)**2*x1**4*a1*(10*a1**2*kc**
2 - 5*a1*ka*vm + 4*a2*ka*vm + 3*b1*ka*vm + 3*b2*ka*vm) + 30*df(x1,t,2)*df(x1,t)**2*x1**3*x2**2*a1**3*( - ka + kc) + df(
x1,t,2)*df(x1,t)**2*x1**3*x2*a1*( - 30*a1**2*kc**2 + 5*a1*ka*vm - 4*a2*ka*vm - 3*b1*ka*vm - 3*b2*ka*vm) + 10*df(x1,t,2)*
df(x1,t)**2*x1**2*x2**3*a1**3*(ka - kc) + 30*df(x1,t,2)*df(x1,t)**2*x1**2*x2**2*a1**3*kc**2 - 10*df(x1,t,2)*df(x1,t)**2*
x1*x2**3*a1**3*kc**2 + 5*df(x1,t,2)*df(x1,t)*x1**6*a1**4*( - ka + kc) + 20*df(x1,t,2)*df(x1,t)*x1**5*x2*a1**4*(ka - kc) 
+ df(x1,t,2)*df(x1,t)*x1**5*a1**2*(5*a1**2*kc**2 - 3*a1*ka*vm + 8*a2*ka*vm + 3*b1*ka*vm + 3*b2*ka*vm) + 30*df(x1,t,2)*df
(x1,t)*x1**4*x2**2*a1**4*( - ka + kc) + 2*df(x1,t,2)*df(x1,t)*x1**4*x2*a1**2*( - 10*a1**2*kc**2 + 4*a1*ka*vm - 8*a2*ka*
vm - 3*b1*ka*vm - 3*b2*ka*vm) + 20*df(x1,t,2)*df(x1,t)*x1**3*x2**3*a1**4*(ka - kc) + df(x1,t,2)*df(x1,t)*x1**3*x2**2*a1
**2*(30*a1**2*kc**2 - 7*a1*ka*vm + 8*a2*ka*vm + 3*b1*ka*vm + 3*b2*ka*vm) + 5*df(x1,t,2)*df(x1,t)*x1**2*x2**4*a1**4*( - 
ka + kc) + 2*df(x1,t,2)*df(x1,t)*x1**2*x2**3*a1**3*( - 10*a1*kc**2 + ka*vm) + 5*df(x1,t,2)*df(x1,t)*x1*x2**4*a1**4*kc**2
 + df(x1,t,2)*x1**7*a1**5*( - ka + kc) + 5*df(x1,t,2)*x1**6*x2*a1**5*(ka - kc) + df(x1,t,2)*x1**6*a1**3*(a1**2*kc**2 + 4
*a2*ka*vm + b1*ka*vm + b2*ka*vm) + 10*df(x1,t,2)*x1**5*x2**2*a1**5*( - ka + kc) + df(x1,t,2)*x1**5*x2*a1**3*( - 5*a1**2*
kc**2 + a1*ka*vm - 12*a2*ka*vm - 3*b1*ka*vm - 3*b2*ka*vm) + 10*df(x1,t,2)*x1**4*x2**3*a1**5*(ka - kc) + df(x1,t,2)*x1**4
*x2**2*a1**3*(10*a1**2*kc**2 - 3*a1*ka*vm + 12*a2*ka*vm + 3*b1*ka*vm + 3*b2*ka*vm) + 5*df(x1,t,2)*x1**3*x2**4*a1**5*( - 
ka + kc) + df(x1,t,2)*x1**3*x2**3*a1**3*( - 10*a1**2*kc**2 + 3*a1*ka*vm - 4*a2*ka*vm - b1*ka*vm - b2*ka*vm) + df(x1,t,2)
*x1**2*x2**5*a1**5*(ka - kc) + df(x1,t,2)*x1**2*x2**4*a1**4*(5*a1*kc**2 - ka*vm) - df(x1,t,2)*x1*x2**5*a1**5*kc**2 - df(
x1,t)**7*kc**2 + df(x1,t)**6*x1**2*(a1*kc - b1*ka - b2*ka + b2*kc) + df(x1,t)**6*x1*kc**2*( - 5*a1 + b2) + 6*df(x1,t)**6
*x2*a1*kc**2 + df(x1,t)**5*x1**3*a1*(5*a1*kc - a2*kc - 5*b1*ka - 5*b2*ka + 6*b2*kc) + df(x1,t)**5*x1**2*x2*a1*( - 5*a1*
kc + a2*kc + 5*b1*ka + 5*b2*ka - 6*b2*kc) + df(x1,t)**5*x1**2*( - 10*a1**2*kc**2 - a1*a2*kc**2 + 6*a1*b2*kc**2 + 2*a1*ka
*vm - b1*ka*vm - b2*ka*vm + b2*kc*vm) + df(x1,t)**5*x1*x2*a1*kc**2*(25*a1 + a2 - 6*b2) - 15*df(x1,t)**5*x2**2*a1**2*kc**
2 + 5*df(x1,t)**4*x1**4*a1**2*(2*a1*kc - a2*kc - 2*b1*ka - 2*b2*ka + 3*b2*kc) + 10*df(x1,t)**4*x1**3*x2*a1**2*( - 2*a1*
kc + a2*kc + 2*b1*ka + 2*b2*ka - 3*b2*kc) + df(x1,t)**4*x1**3*a1*( - 10*a1**2*kc**2 - 5*a1*a2*kc**2 + 15*a1*b2*kc**2 + 4
*a1*ka*vm - 3*a2*ka*vm - 3*b1*ka*vm - 3*b2*ka*vm + 5*b2*kc*vm) + 5*df(x1,t)**4*x1**2*x2**2*a1**2*(2*a1*kc - a2*kc - 2*b1
*ka - 2*b2*ka + 3*b2*kc) + df(x1,t)**4*x1**2*x2*a1*(40*a1**2*kc**2 + 10*a1*a2*kc**2 - 30*a1*b2*kc**2 - 6*a1*ka*vm + 2*a2
*ka*vm + 4*b1*ka*vm + 4*b2*ka*vm - 5*b2*kc*vm) + 5*df(x1,t)**4*x1*x2**2*a1**2*kc**2*( - 10*a1 - a2 + 3*b2) + 20*df(x1,t)
**4*x2**3*a1**3*kc**2 + 10*df(x1,t)**3*x1**5*a1**3*(a1*kc - a2*kc - b1*ka - b2*ka + 2*b2*kc) + 30*df(x1,t)**3*x1**4*x2*
a1**3*( - a1*kc + a2*kc + b1*ka + b2*ka - 2*b2*kc) + df(x1,t)**3*x1**4*a1*( - 5*a1**3*kc**2 - 10*a1**2*a2*kc**2 + 20*a1
**2*b2*kc**2 + 2*a1**2*ka*vm - 5*a1*a2*ka*vm - 3*a1*b1*ka*vm - 3*a1*b2*ka*vm + 10*a1*b2*kc*vm + a2**2*ka*vm - a2*b1*ka*
vm - a2*b2*ka*vm) + 30*df(x1,t)**3*x1**3*x2**2*a1**3*(a1*kc - a2*kc - b1*ka - b2*ka + 2*b2*kc) + df(x1,t)**3*x1**3*x2*a1
*(30*a1**3*kc**2 + 30*a1**2*a2*kc**2 - 60*a1**2*b2*kc**2 - 8*a1**2*ka*vm + 11*a1*a2*ka*vm + 9*a1*b1*ka*vm + 9*a1*b2*ka*
vm - 20*a1*b2*kc*vm - a2**2*ka*vm + a2*b1*ka*vm + a2*b2*ka*vm) + 10*df(x1,t)**3*x1**2*x2**3*a1**3*( - a1*kc + a2*kc + b1
*ka + b2*ka - 2*b2*kc) + 2*df(x1,t)**3*x1**2*x2**2*a1**2*( - 30*a1**2*kc**2 - 15*a1*a2*kc**2 + 30*a1*b2*kc**2 + 3*a1*ka*
vm - 3*a2*ka*vm - 3*b1*ka*vm - 3*b2*ka*vm + 5*b2*kc*vm) + 10*df(x1,t)**3*x1*x2**3*a1**3*kc**2*(5*a1 + a2 - 2*b2) - 15*df
(x1,t)**3*x2**4*a1**4*kc**2 + 5*df(x1,t)**2*x1**6*a1**4*(a1*kc - 2*a2*kc - b1*ka - b2*ka + 3*b2*kc) + 20*df(x1,t)**2*x1
**5*x2*a1**4*( - a1*kc + 2*a2*kc + b1*ka + b2*ka - 3*b2*kc) + df(x1,t)**2*x1**5*a1**2*( - a1**3*kc**2 - 10*a1**2*a2*kc**
2 + 15*a1**2*b2*kc**2 - a1*a2*ka*vm - a1*b1*ka*vm - a1*b2*ka*vm + 10*a1*b2*kc*vm + a2**2*ka*vm - 3*a2*b1*ka*vm - 3*a2*b2
*ka*vm) + 30*df(x1,t)**2*x1**4*x2**2*a1**4*(a1*kc - 2*a2*kc - b1*ka - b2*ka + 3*b2*kc) + 2*df(x1,t)**2*x1**4*x2*a1**2*(5
*a1**3*kc**2 + 20*a1**2*a2*kc**2 - 30*a1**2*b2*kc**2 - a1**2*ka*vm + 4*a1*a2*ka*vm + 3*a1*b1*ka*vm + 3*a1*b2*ka*vm - 15*
a1*b2*kc*vm - a2**2*ka*vm + 3*a2*b1*ka*vm + 3*a2*b2*ka*vm) + 20*df(x1,t)**2*x1**3*x2**3*a1**4*( - a1*kc + 2*a2*kc + b1*
ka + b2*ka - 3*b2*kc) + df(x1,t)**2*x1**3*x2**2*a1**2*( - 30*a1**3*kc**2 - 60*a1**2*a2*kc**2 + 90*a1**2*b2*kc**2 + 4*a1
**2*ka*vm - 13*a1*a2*ka*vm - 9*a1*b1*ka*vm - 9*a1*b2*ka*vm + 30*a1*b2*kc*vm + a2**2*ka*vm - 3*a2*b1*ka*vm - 3*a2*b2*ka*
vm) + 5*df(x1,t)**2*x1**2*x2**4*a1**4*(a1*kc - 2*a2*kc - b1*ka - b2*ka + 3*b2*kc) + 2*df(x1,t)**2*x1**2*x2**3*a1**3*(20*
a1**2*kc**2 + 20*a1*a2*kc**2 - 30*a1*b2*kc**2 - a1*ka*vm + 3*a2*ka*vm + 2*b1*ka*vm + 2*b2*ka*vm - 5*b2*kc*vm) + 5*df(x1,
t)**2*x1*x2**4*a1**4*kc**2*( - 5*a1 - 2*a2 + 3*b2) + 6*df(x1,t)**2*x2**5*a1**5*kc**2 + df(x1,t)*x1**7*a1**5*(a1*kc - 5*
a2*kc - b1*ka - b2*ka + 6*b2*kc) + 5*df(x1,t)*x1**6*x2*a1**5*( - a1*kc + 5*a2*kc + b1*ka + b2*ka - 6*b2*kc) + df(x1,t)*
x1**6*a1**3*( - 5*a1**2*a2*kc**2 + 6*a1**2*b2*kc**2 + a1*a2*ka*vm + 5*a1*b2*kc*vm - a2**2*ka*vm - 3*a2*b1*ka*vm - 3*a2*
b2*ka*vm) + 10*df(x1,t)*x1**5*x2**2*a1**5*(a1*kc - 5*a2*kc - b1*ka - b2*ka + 6*b2*kc) + df(x1,t)*x1**5*x2*a1**3*(a1**3*
kc**2 + 25*a1**2*a2*kc**2 - 30*a1**2*b2*kc**2 - a1*a2*ka*vm + a1*b1*ka*vm + a1*b2*ka*vm - 20*a1*b2*kc*vm + 3*a2**2*ka*vm
 + 9*a2*b1*ka*vm + 9*a2*b2*ka*vm) + 10*df(x1,t)*x1**4*x2**3*a1**5*( - a1*kc + 5*a2*kc + b1*ka + b2*ka - 6*b2*kc) + df(x1
,t)*x1**4*x2**2*a1**3*( - 5*a1**3*kc**2 - 50*a1**2*a2*kc**2 + 60*a1**2*b2*kc**2 - 3*a1*a2*ka*vm - 3*a1*b1*ka*vm - 3*a1*
b2*ka*vm + 30*a1*b2*kc*vm - 3*a2**2*ka*vm - 9*a2*b1*ka*vm - 9*a2*b2*ka*vm) + 5*df(x1,t)*x1**3*x2**4*a1**5*(a1*kc - 5*a2*
kc - b1*ka - b2*ka + 6*b2*kc) + df(x1,t)*x1**3*x2**3*a1**3*(10*a1**3*kc**2 + 50*a1**2*a2*kc**2 - 60*a1**2*b2*kc**2 + 5*
a1*a2*ka*vm + 3*a1*b1*ka*vm + 3*a1*b2*ka*vm - 20*a1*b2*kc*vm + a2**2*ka*vm + 3*a2*b1*ka*vm + 3*a2*b2*ka*vm) + df(x1,t)*
x1**2*x2**5*a1**5*( - a1*kc + 5*a2*kc + b1*ka + b2*ka - 6*b2*kc) + df(x1,t)*x1**2*x2**4*a1**4*( - 10*a1**2*kc**2 - 25*a1
*a2*kc**2 + 30*a1*b2*kc**2 - 2*a2*ka*vm - b1*ka*vm - b2*ka*vm + 5*b2*kc*vm) + df(x1,t)*x1*x2**5*a1**5*kc**2*(5*a1 + 5*a2
 - 6*b2) - df(x1,t)*x2**6*a1**6*kc**2 + x1**8*a1**6*kc*( - a2 + b2) + 6*x1**7*x2*a1**6*kc*(a2 - b2) + x1**7*a1**4*( - a1
**2*a2*kc**2 + a1**2*b2*kc**2 + a1*b2*kc*vm - a2**2*ka*vm - a2*b1*ka*vm - a2*b2*ka*vm) + 15*x1**6*x2**2*a1**6*kc*( - a2 
+ b2) + x1**6*x2*a1**4*(6*a1**2*a2*kc**2 - 6*a1**2*b2*kc**2 - 5*a1*b2*kc*vm + 4*a2**2*ka*vm + 4*a2*b1*ka*vm + 4*a2*b2*ka
*vm) + 20*x1**5*x2**3*a1**6*kc*(a2 - b2) + x1**5*x2**2*a1**4*( - 15*a1**2*a2*kc**2 + 15*a1**2*b2*kc**2 + 10*a1*b2*kc*vm 
- 6*a2**2*ka*vm - 6*a2*b1*ka*vm - 6*a2*b2*ka*vm) + 15*x1**4*x2**4*a1**6*kc*( - a2 + b2) + 2*x1**4*x2**3*a1**4*(10*a1**2*
a2*kc**2 - 10*a1**2*b2*kc**2 - 5*a1*b2*kc*vm + 2*a2**2*ka*vm + 2*a2*b1*ka*vm + 2*a2*b2*ka*vm) + 6*x1**3*x2**5*a1**6*kc*(
a2 - b2) + x1**3*x2**4*a1**4*( - 15*a1**2*a2*kc**2 + 15*a1**2*b2*kc**2 + 5*a1*b2*kc*vm - a2**2*ka*vm - a2*b1*ka*vm - a2*
b2*ka*vm) + x1**2*x2**6*a1**6*kc*( - a2 + b2) + x1**2*x2**5*a1**5*kc*(6*a1*a2*kc - 6*a1*b2*kc - b2*vm) + x1*x2**6*a1**6*
kc**2*( - a2 + b2)$

aa_(3) := df(x1,t)*x1*ka + df(x1,t)*x3*kc + df(x1,t)*ka*kc + x1**2*a1*ka - x1*x2*a1*ka + x1*x3*a1*kc + x1*ka*(a1*kc + vm
) - x2*x3*a1*kc - x2*a1*ka*kc$

aa_(4) :=  - df(x1,t,2)*x1**2*ka*vm + df(x1,t)**3*x1*(ka - kc) - df(x1,t)**3*kc**2 + df(x1,t)**2*x1**2*(2*a1*ka - 3*a1*
kc + b1*ka) + df(x1,t)**2*x1*x2*a1*( - 2*ka + 3*kc) + df(x1,t)**2*x1*x4*b1*kc + df(x1,t)**2*x1*( - 3*a1*kc**2 + b1*ka*kc
 + ka*vm - kc*vm) + 3*df(x1,t)**2*x2*a1*kc**2 + df(x1,t)*x1**3*a1*(a1*ka - 3*a1*kc + 2*b1*ka) + 2*df(x1,t)*x1**2*x2*a1*(
 - a1*ka + 3*a1*kc - b1*ka) + 2*df(x1,t)*x1**2*x4*a1*b1*kc + df(x1,t)*x1**2*( - 3*a1**2*kc**2 + 2*a1*b1*ka*kc - 2*a1*kc*
vm + b1*ka*vm) + df(x1,t)*x1*x2**2*a1**2*(ka - 3*kc) - 2*df(x1,t)*x1*x2*x4*a1*b1*kc + df(x1,t)*x1*x2*a1*(6*a1*kc**2 - 2*
b1*ka*kc - ka*vm + 2*kc*vm) - 3*df(x1,t)*x2**2*a1**2*kc**2 + x1**4*a1**2*( - a1*kc + b1*ka) + x1**3*x2*a1**2*(3*a1*kc - 
2*b1*ka) + x1**3*x4*a1**2*b1*kc + x1**3*a1*( - a1**2*kc**2 + a1*b1*ka*kc - a1*kc*vm + a2*ka*vm + b1*ka*vm) + x1**2*x2**2
*a1**2*( - 3*a1*kc + b1*ka) - 2*x1**2*x2*x4*a1**2*b1*kc + x1**2*x2*a1*(3*a1**2*kc**2 - 2*a1*b1*ka*kc + 2*a1*kc*vm - a2*
ka*vm - b1*ka*vm) + x1*x2**3*a1**3*kc + x1*x2**2*x4*a1**2*b1*kc + x1*x2**2*a1**2*kc*( - 3*a1*kc + b1*ka - vm) + x2**3*a1
**3*kc**2$

aa_(5) :=  - x1 + y1$

aa_(6) :=  - df(x1,t,2)*x1**2*ka*vm + df(x1,t)**3*x1*(ka - kc) - df(x1,t)**3*kc**2 + df(x1,t)**2*x1**2*(2*a1*ka - 3*a1*
kc + b1*ka) + df(x1,t)**2*x1*x2*a1*( - 2*ka + 3*kc) + df(x1,t)**2*x1*y2*b1*kc + df(x1,t)**2*x1*( - 3*a1*kc**2 + b1*ka*kc
 + ka*vm - kc*vm) + 3*df(x1,t)**2*x2*a1*kc**2 + df(x1,t)*x1**3*a1*(a1*ka - 3*a1*kc + 2*b1*ka) + 2*df(x1,t)*x1**2*x2*a1*(
 - a1*ka + 3*a1*kc - b1*ka) + 2*df(x1,t)*x1**2*y2*a1*b1*kc + df(x1,t)*x1**2*( - 3*a1**2*kc**2 + 2*a1*b1*ka*kc - 2*a1*kc*
vm + b1*ka*vm) + df(x1,t)*x1*x2**2*a1**2*(ka - 3*kc) - 2*df(x1,t)*x1*x2*y2*a1*b1*kc + df(x1,t)*x1*x2*a1*(6*a1*kc**2 - 2*
b1*ka*kc - ka*vm + 2*kc*vm) - 3*df(x1,t)*x2**2*a1**2*kc**2 + x1**4*a1**2*( - a1*kc + b1*ka) + x1**3*x2*a1**2*(3*a1*kc - 
2*b1*ka) + x1**3*y2*a1**2*b1*kc + x1**3*a1*( - a1**2*kc**2 + a1*b1*ka*kc - a1*kc*vm + a2*ka*vm + b1*ka*vm) + x1**2*x2**2
*a1**2*( - 3*a1*kc + b1*ka) - 2*x1**2*x2*y2*a1**2*b1*kc + x1**2*x2*a1*(3*a1**2*kc**2 - 2*a1*b1*ka*kc + 2*a1*kc*vm - a2*
ka*vm - b1*ka*vm) + x1*x2**3*a1**3*kc + x1*x2**2*y2*a1**2*b1*kc + x1*x2**2*a1**2*kc*( - 3*a1*kc + b1*ka - vm) + x2**3*a1
**3*kc**2$

MODEL ALGEBRAICALLY OBSERVABLE$

NORMALIZED  INPUT /OUTPUT RELATION(S) $

aan_(1) := df(x2,t) - x1*a2 + x2*a2$

aan_(2) := (df(x1,t,3)*df(x1,t)**3*x1**3*ka*vm + 3*df(x1,t,3)*df(x1,t)**2*x1**4*a1*ka*vm - 3*df(x1,t,3)*df(x1,t)**2*x1**
3*x2*a1*ka*vm + 3*df(x1,t,3)*df(x1,t)*x1**5*a1**2*ka*vm - 6*df(x1,t,3)*df(x1,t)*x1**4*x2*a1**2*ka*vm + 3*df(x1,t,3)*df(
x1,t)*x1**3*x2**2*a1**2*ka*vm + df(x1,t,3)*x1**6*a1**3*ka*vm - 3*df(x1,t,3)*x1**5*x2*a1**3*ka*vm + 3*df(x1,t,3)*x1**4*x2
**2*a1**3*ka*vm - df(x1,t,3)*x1**3*x2**3*a1**3*ka*vm - 2*df(x1,t,2)**2*df(x1,t)**2*x1**3*ka*vm - 4*df(x1,t,2)**2*df(x1,t
)*x1**4*a1*ka*vm + 4*df(x1,t,2)**2*df(x1,t)*x1**3*x2*a1*ka*vm - 2*df(x1,t,2)**2*x1**5*a1**2*ka*vm + 4*df(x1,t,2)**2*x1**
4*x2*a1**2*ka*vm - 2*df(x1,t,2)**2*x1**3*x2**2*a1**2*ka*vm + df(x1,t,2)*df(x1,t)**5*x1**2*( - ka + kc) + df(x1,t,2)*df(
x1,t)**5*x1*kc**2 + 5*df(x1,t,2)*df(x1,t)**4*x1**3*a1*( - ka + kc) + 5*df(x1,t,2)*df(x1,t)**4*x1**2*x2*a1*(ka - kc) + df
(x1,t,2)*df(x1,t)**4*x1**2*(5*a1*kc**2 + ka*vm) - 5*df(x1,t,2)*df(x1,t)**4*x1*x2*a1*kc**2 + 10*df(x1,t,2)*df(x1,t)**3*x1
**4*a1**2*( - ka + kc) + 20*df(x1,t,2)*df(x1,t)**3*x1**3*x2*a1**2*(ka - kc) + df(x1,t,2)*df(x1,t)**3*x1**3*(10*a1**2*kc
**2 - a1*ka*vm + b1*ka*vm + b2*ka*vm) + 10*df(x1,t,2)*df(x1,t)**3*x1**2*x2**2*a1**2*( - ka + kc) + 2*df(x1,t,2)*df(x1,t)
**3*x1**2*x2*a1*( - 10*a1*kc**2 - ka*vm) + 10*df(x1,t,2)*df(x1,t)**3*x1*x2**2*a1**2*kc**2 + 10*df(x1,t,2)*df(x1,t)**2*x1
**5*a1**3*( - ka + kc) + 30*df(x1,t,2)*df(x1,t)**2*x1**4*x2*a1**3*(ka - kc) + df(x1,t,2)*df(x1,t)**2*x1**4*a1*(10*a1**2*
kc**2 - 5*a1*ka*vm + 4*a2*ka*vm + 3*b1*ka*vm + 3*b2*ka*vm) + 30*df(x1,t,2)*df(x1,t)**2*x1**3*x2**2*a1**3*( - ka + kc) + 
df(x1,t,2)*df(x1,t)**2*x1**3*x2*a1*( - 30*a1**2*kc**2 + 5*a1*ka*vm - 4*a2*ka*vm - 3*b1*ka*vm - 3*b2*ka*vm) + 10*df(x1,t,
2)*df(x1,t)**2*x1**2*x2**3*a1**3*(ka - kc) + 30*df(x1,t,2)*df(x1,t)**2*x1**2*x2**2*a1**3*kc**2 - 10*df(x1,t,2)*df(x1,t)
**2*x1*x2**3*a1**3*kc**2 + 5*df(x1,t,2)*df(x1,t)*x1**6*a1**4*( - ka + kc) + 20*df(x1,t,2)*df(x1,t)*x1**5*x2*a1**4*(ka - 
kc) + df(x1,t,2)*df(x1,t)*x1**5*a1**2*(5*a1**2*kc**2 - 3*a1*ka*vm + 8*a2*ka*vm + 3*b1*ka*vm + 3*b2*ka*vm) + 30*df(x1,t,2
)*df(x1,t)*x1**4*x2**2*a1**4*( - ka + kc) + 2*df(x1,t,2)*df(x1,t)*x1**4*x2*a1**2*( - 10*a1**2*kc**2 + 4*a1*ka*vm - 8*a2*
ka*vm - 3*b1*ka*vm - 3*b2*ka*vm) + 20*df(x1,t,2)*df(x1,t)*x1**3*x2**3*a1**4*(ka - kc) + df(x1,t,2)*df(x1,t)*x1**3*x2**2*
a1**2*(30*a1**2*kc**2 - 7*a1*ka*vm + 8*a2*ka*vm + 3*b1*ka*vm + 3*b2*ka*vm) + 5*df(x1,t,2)*df(x1,t)*x1**2*x2**4*a1**4*( -
 ka + kc) + 2*df(x1,t,2)*df(x1,t)*x1**2*x2**3*a1**3*( - 10*a1*kc**2 + ka*vm) + 5*df(x1,t,2)*df(x1,t)*x1*x2**4*a1**4*kc**
2 + df(x1,t,2)*x1**7*a1**5*( - ka + kc) + 5*df(x1,t,2)*x1**6*x2*a1**5*(ka - kc) + df(x1,t,2)*x1**6*a1**3*(a1**2*kc**2 + 
4*a2*ka*vm + b1*ka*vm + b2*ka*vm) + 10*df(x1,t,2)*x1**5*x2**2*a1**5*( - ka + kc) + df(x1,t,2)*x1**5*x2*a1**3*( - 5*a1**2
*kc**2 + a1*ka*vm - 12*a2*ka*vm - 3*b1*ka*vm - 3*b2*ka*vm) + 10*df(x1,t,2)*x1**4*x2**3*a1**5*(ka - kc) + df(x1,t,2)*x1**
4*x2**2*a1**3*(10*a1**2*kc**2 - 3*a1*ka*vm + 12*a2*ka*vm + 3*b1*ka*vm + 3*b2*ka*vm) + 5*df(x1,t,2)*x1**3*x2**4*a1**5*( -
 ka + kc) + df(x1,t,2)*x1**3*x2**3*a1**3*( - 10*a1**2*kc**2 + 3*a1*ka*vm - 4*a2*ka*vm - b1*ka*vm - b2*ka*vm) + df(x1,t,2
)*x1**2*x2**5*a1**5*(ka - kc) + df(x1,t,2)*x1**2*x2**4*a1**4*(5*a1*kc**2 - ka*vm) - df(x1,t,2)*x1*x2**5*a1**5*kc**2 - df
(x1,t)**7*kc**2 + df(x1,t)**6*x1**2*(a1*kc - b1*ka - b2*ka + b2*kc) + df(x1,t)**6*x1*kc**2*( - 5*a1 + b2) + 6*df(x1,t)**
6*x2*a1*kc**2 + df(x1,t)**5*x1**3*a1*(5*a1*kc - a2*kc - 5*b1*ka - 5*b2*ka + 6*b2*kc) + df(x1,t)**5*x1**2*x2*a1*( - 5*a1*
kc + a2*kc + 5*b1*ka + 5*b2*ka - 6*b2*kc) + df(x1,t)**5*x1**2*( - 10*a1**2*kc**2 - a1*a2*kc**2 + 6*a1*b2*kc**2 + 2*a1*ka
*vm - b1*ka*vm - b2*ka*vm + b2*kc*vm) + df(x1,t)**5*x1*x2*a1*kc**2*(25*a1 + a2 - 6*b2) - 15*df(x1,t)**5*x2**2*a1**2*kc**
2 + 5*df(x1,t)**4*x1**4*a1**2*(2*a1*kc - a2*kc - 2*b1*ka - 2*b2*ka + 3*b2*kc) + 10*df(x1,t)**4*x1**3*x2*a1**2*( - 2*a1*
kc + a2*kc + 2*b1*ka + 2*b2*ka - 3*b2*kc) + df(x1,t)**4*x1**3*a1*( - 10*a1**2*kc**2 - 5*a1*a2*kc**2 + 15*a1*b2*kc**2 + 4
*a1*ka*vm - 3*a2*ka*vm - 3*b1*ka*vm - 3*b2*ka*vm + 5*b2*kc*vm) + 5*df(x1,t)**4*x1**2*x2**2*a1**2*(2*a1*kc - a2*kc - 2*b1
*ka - 2*b2*ka + 3*b2*kc) + df(x1,t)**4*x1**2*x2*a1*(40*a1**2*kc**2 + 10*a1*a2*kc**2 - 30*a1*b2*kc**2 - 6*a1*ka*vm + 2*a2
*ka*vm + 4*b1*ka*vm + 4*b2*ka*vm - 5*b2*kc*vm) + 5*df(x1,t)**4*x1*x2**2*a1**2*kc**2*( - 10*a1 - a2 + 3*b2) + 20*df(x1,t)
**4*x2**3*a1**3*kc**2 + 10*df(x1,t)**3*x1**5*a1**3*(a1*kc - a2*kc - b1*ka - b2*ka + 2*b2*kc) + 30*df(x1,t)**3*x1**4*x2*
a1**3*( - a1*kc + a2*kc + b1*ka + b2*ka - 2*b2*kc) + df(x1,t)**3*x1**4*a1*( - 5*a1**3*kc**2 - 10*a1**2*a2*kc**2 + 20*a1
**2*b2*kc**2 + 2*a1**2*ka*vm - 5*a1*a2*ka*vm - 3*a1*b1*ka*vm - 3*a1*b2*ka*vm + 10*a1*b2*kc*vm + a2**2*ka*vm - a2*b1*ka*
vm - a2*b2*ka*vm) + 30*df(x1,t)**3*x1**3*x2**2*a1**3*(a1*kc - a2*kc - b1*ka - b2*ka + 2*b2*kc) + df(x1,t)**3*x1**3*x2*a1
*(30*a1**3*kc**2 + 30*a1**2*a2*kc**2 - 60*a1**2*b2*kc**2 - 8*a1**2*ka*vm + 11*a1*a2*ka*vm + 9*a1*b1*ka*vm + 9*a1*b2*ka*
vm - 20*a1*b2*kc*vm - a2**2*ka*vm + a2*b1*ka*vm + a2*b2*ka*vm) + 10*df(x1,t)**3*x1**2*x2**3*a1**3*( - a1*kc + a2*kc + b1
*ka + b2*ka - 2*b2*kc) + 2*df(x1,t)**3*x1**2*x2**2*a1**2*( - 30*a1**2*kc**2 - 15*a1*a2*kc**2 + 30*a1*b2*kc**2 + 3*a1*ka*
vm - 3*a2*ka*vm - 3*b1*ka*vm - 3*b2*ka*vm + 5*b2*kc*vm) + 10*df(x1,t)**3*x1*x2**3*a1**3*kc**2*(5*a1 + a2 - 2*b2) - 15*df
(x1,t)**3*x2**4*a1**4*kc**2 + 5*df(x1,t)**2*x1**6*a1**4*(a1*kc - 2*a2*kc - b1*ka - b2*ka + 3*b2*kc) + 20*df(x1,t)**2*x1
**5*x2*a1**4*( - a1*kc + 2*a2*kc + b1*ka + b2*ka - 3*b2*kc) + df(x1,t)**2*x1**5*a1**2*( - a1**3*kc**2 - 10*a1**2*a2*kc**
2 + 15*a1**2*b2*kc**2 - a1*a2*ka*vm - a1*b1*ka*vm - a1*b2*ka*vm + 10*a1*b2*kc*vm + a2**2*ka*vm - 3*a2*b1*ka*vm - 3*a2*b2
*ka*vm) + 30*df(x1,t)**2*x1**4*x2**2*a1**4*(a1*kc - 2*a2*kc - b1*ka - b2*ka + 3*b2*kc) + 2*df(x1,t)**2*x1**4*x2*a1**2*(5
*a1**3*kc**2 + 20*a1**2*a2*kc**2 - 30*a1**2*b2*kc**2 - a1**2*ka*vm + 4*a1*a2*ka*vm + 3*a1*b1*ka*vm + 3*a1*b2*ka*vm - 15*
a1*b2*kc*vm - a2**2*ka*vm + 3*a2*b1*ka*vm + 3*a2*b2*ka*vm) + 20*df(x1,t)**2*x1**3*x2**3*a1**4*( - a1*kc + 2*a2*kc + b1*
ka + b2*ka - 3*b2*kc) + df(x1,t)**2*x1**3*x2**2*a1**2*( - 30*a1**3*kc**2 - 60*a1**2*a2*kc**2 + 90*a1**2*b2*kc**2 + 4*a1
**2*ka*vm - 13*a1*a2*ka*vm - 9*a1*b1*ka*vm - 9*a1*b2*ka*vm + 30*a1*b2*kc*vm + a2**2*ka*vm - 3*a2*b1*ka*vm - 3*a2*b2*ka*
vm) + 5*df(x1,t)**2*x1**2*x2**4*a1**4*(a1*kc - 2*a2*kc - b1*ka - b2*ka + 3*b2*kc) + 2*df(x1,t)**2*x1**2*x2**3*a1**3*(20*
a1**2*kc**2 + 20*a1*a2*kc**2 - 30*a1*b2*kc**2 - a1*ka*vm + 3*a2*ka*vm + 2*b1*ka*vm + 2*b2*ka*vm - 5*b2*kc*vm) + 5*df(x1,
t)**2*x1*x2**4*a1**4*kc**2*( - 5*a1 - 2*a2 + 3*b2) + 6*df(x1,t)**2*x2**5*a1**5*kc**2 + df(x1,t)*x1**7*a1**5*(a1*kc - 5*
a2*kc - b1*ka - b2*ka + 6*b2*kc) + 5*df(x1,t)*x1**6*x2*a1**5*( - a1*kc + 5*a2*kc + b1*ka + b2*ka - 6*b2*kc) + df(x1,t)*
x1**6*a1**3*( - 5*a1**2*a2*kc**2 + 6*a1**2*b2*kc**2 + a1*a2*ka*vm + 5*a1*b2*kc*vm - a2**2*ka*vm - 3*a2*b1*ka*vm - 3*a2*
b2*ka*vm) + 10*df(x1,t)*x1**5*x2**2*a1**5*(a1*kc - 5*a2*kc - b1*ka - b2*ka + 6*b2*kc) + df(x1,t)*x1**5*x2*a1**3*(a1**3*
kc**2 + 25*a1**2*a2*kc**2 - 30*a1**2*b2*kc**2 - a1*a2*ka*vm + a1*b1*ka*vm + a1*b2*ka*vm - 20*a1*b2*kc*vm + 3*a2**2*ka*vm
 + 9*a2*b1*ka*vm + 9*a2*b2*ka*vm) + 10*df(x1,t)*x1**4*x2**3*a1**5*( - a1*kc + 5*a2*kc + b1*ka + b2*ka - 6*b2*kc) + df(x1
,t)*x1**4*x2**2*a1**3*( - 5*a1**3*kc**2 - 50*a1**2*a2*kc**2 + 60*a1**2*b2*kc**2 - 3*a1*a2*ka*vm - 3*a1*b1*ka*vm - 3*a1*
b2*ka*vm + 30*a1*b2*kc*vm - 3*a2**2*ka*vm - 9*a2*b1*ka*vm - 9*a2*b2*ka*vm) + 5*df(x1,t)*x1**3*x2**4*a1**5*(a1*kc - 5*a2*
kc - b1*ka - b2*ka + 6*b2*kc) + df(x1,t)*x1**3*x2**3*a1**3*(10*a1**3*kc**2 + 50*a1**2*a2*kc**2 - 60*a1**2*b2*kc**2 + 5*
a1*a2*ka*vm + 3*a1*b1*ka*vm + 3*a1*b2*ka*vm - 20*a1*b2*kc*vm + a2**2*ka*vm + 3*a2*b1*ka*vm + 3*a2*b2*ka*vm) + df(x1,t)*
x1**2*x2**5*a1**5*( - a1*kc + 5*a2*kc + b1*ka + b2*ka - 6*b2*kc) + df(x1,t)*x1**2*x2**4*a1**4*( - 10*a1**2*kc**2 - 25*a1
*a2*kc**2 + 30*a1*b2*kc**2 - 2*a2*ka*vm - b1*ka*vm - b2*ka*vm + 5*b2*kc*vm) + df(x1,t)*x1*x2**5*a1**5*kc**2*(5*a1 + 5*a2
 - 6*b2) - df(x1,t)*x2**6*a1**6*kc**2 + x1**8*a1**6*kc*( - a2 + b2) + 6*x1**7*x2*a1**6*kc*(a2 - b2) + x1**7*a1**4*( - a1
**2*a2*kc**2 + a1**2*b2*kc**2 + a1*b2*kc*vm - a2**2*ka*vm - a2*b1*ka*vm - a2*b2*ka*vm) + 15*x1**6*x2**2*a1**6*kc*( - a2 
+ b2) + x1**6*x2*a1**4*(6*a1**2*a2*kc**2 - 6*a1**2*b2*kc**2 - 5*a1*b2*kc*vm + 4*a2**2*ka*vm + 4*a2*b1*ka*vm + 4*a2*b2*ka
*vm) + 20*x1**5*x2**3*a1**6*kc*(a2 - b2) + x1**5*x2**2*a1**4*( - 15*a1**2*a2*kc**2 + 15*a1**2*b2*kc**2 + 10*a1*b2*kc*vm 
- 6*a2**2*ka*vm - 6*a2*b1*ka*vm - 6*a2*b2*ka*vm) + 15*x1**4*x2**4*a1**6*kc*( - a2 + b2) + 2*x1**4*x2**3*a1**4*(10*a1**2*
a2*kc**2 - 10*a1**2*b2*kc**2 - 5*a1*b2*kc*vm + 2*a2**2*ka*vm + 2*a2*b1*ka*vm + 2*a2*b2*ka*vm) + 6*x1**3*x2**5*a1**6*kc*(
a2 - b2) + x1**3*x2**4*a1**4*( - 15*a1**2*a2*kc**2 + 15*a1**2*b2*kc**2 + 5*a1*b2*kc*vm - a2**2*ka*vm - a2*b1*ka*vm - a2*
b2*ka*vm) + x1**2*x2**6*a1**6*kc*( - a2 + b2) + x1**2*x2**5*a1**5*kc*(6*a1*a2*kc - 6*a1*b2*kc - b2*vm) + x1*x2**6*a1**6*
kc**2*( - a2 + b2))/(ka*vm)$

RANDOMLY CHOSEN NUMERICAL PARAMETER(S) VECTOR$

b2_ := {a1=8,a2=16,b1=3,b2=20,ka=17,kc=13,vm=5}$

EXHAUSTIVE SUMMARY $

flist_ := { - a2 + 16,
a2 - 16,
a1 - 8,
 - a1 + 8,
a1**2 - 64,
 - a1**2 + 64,
a1**2 - 64,
a1**3 - 512,
 - a1**3 + 512,
a1**3 - 512,
 - a1**3 + 512,
 - a1 + 8,
a1 - 8,
 - a1**2 + 64,
a1**2 - 64,
 - a1**2 + 64,
( - 169*ka*vm + 85*kc**2)/(85*ka*vm),
(85*a1*kc**2 - 1352*ka*vm)/(17*ka*vm),
( - 85*a1*kc**2 + 1352*ka*vm)/(17*ka*vm),
(2*(85*a1**2*kc**2 - 10816*ka*vm))/(17*ka*vm),
(6*(85*a1**3*kc**2 - 86528*ka*vm))/(17*ka*vm),
(2*( - 85*a1**3*kc**2 + 86528*ka*vm))/(17*ka*vm),
(85*a1**4*kc**2 - 692224*ka*vm)/(17*ka*vm),
( - 85*a1**5*kc**2 + 5537792*ka*vm)/(85*ka*vm),
(169*ka*vm - 85*kc**2)/(85*ka*vm),
(6*(85*a1*kc**2 - 1352*ka*vm))/(85*ka*vm),
(3*( - 85*a1**2*kc**2 + 10816*ka*vm))/(17*ka*vm),
(4*(85*a1**3*kc**2 - 86528*ka*vm))/(17*ka*vm),
(3*( - 85*a1**4*kc**2 + 692224*ka*vm))/(17*ka*vm),
(6*(85*a1**5*kc**2 - 5537792*ka*vm))/(85*ka*vm),
( - 85*a1**6*kc**2 + 44302336*ka*vm)/(85*ka*vm),
(6*( - 85*a1**3*ka + 85*a1**3*kc + 2048*ka*vm))/(17*ka*vm),
(2*(85*a1**3*ka - 85*a1**3*kc - 2048*ka*vm))/(17*ka*vm),
(4*ka*vm - 85*ka + 85*kc)/(85*ka*vm),
( - 85*a1*ka + 85*a1*kc + 32*ka*vm)/(17*ka*vm),
( - 85*a1**4*ka + 85*a1**4*kc + 16384*ka*vm)/(17*ka*vm),
(4*(85*a1**4*ka - 85*a1**4*kc - 16384*ka*vm))/(17*ka*vm),
(6*( - 85*a1**4*ka + 85*a1**4*kc + 16384*ka*vm))/(17*ka*vm),
(4*(85*a1**4*ka - 85*a1**4*kc - 16384*ka*vm))/(17*ka*vm),
( - 85*a1**4*ka + 85*a1**4*kc + 16384*ka*vm)/(17*ka*vm),
(2*( - 170*a1**4*kc**2 + 17*a1**3*ka*vm + 1375744*ka*vm))/(17*ka*vm),
(2*( - 85*a1**2*ka + 85*a1**2*kc + 256*ka*vm))/(17*ka*vm),
( - 85*a1**5*ka + 85*a1**5*kc + 131072*ka*vm)/(85*ka*vm),
(85*a1**5*ka - 85*a1**5*kc - 131072*ka*vm)/(17*ka*vm),
(2*( - 85*a1**5*ka + 85*a1**5*kc + 131072*ka*vm))/(17*ka*vm),
(2*(85*a1**5*ka - 85*a1**5*kc - 131072*ka*vm))/(17*ka*vm),
( - 85*a1**5*ka + 85*a1**5*kc + 131072*ka*vm)/(17*ka*vm),
(85*a1**5*ka - 85*a1**5*kc - 131072*ka*vm)/(85*ka*vm),
(85*a1**5*kc**2 - 17*a1**4*ka*vm - 5468160*ka*vm)/(17*ka*vm),
(4*(85*a1**2*ka - 85*a1**2*kc - 256*ka*vm))/(17*ka*vm),
( - 85*a1*kc**2 + 17*b2*kc**2 + 676*ka*vm)/(17*ka*vm),
(2*( - 85*a1**2*ka + 85*a1**2*kc + 256*ka*vm))/(17*ka*vm),
(2*( - 170*a1**2*kc**2 - 17*a1*ka*vm + 21768*ka*vm))/(17*ka*vm),
(85*a1*ka - 85*a1*kc - 32*ka*vm)/(17*ka*vm),
(2*( - 85*a1**3*ka + 85*a1**3*kc + 2048*ka*vm))/(17*ka*vm),
(6*(85*a1**3*ka - 85*a1**3*kc - 2048*ka*vm))/(17*ka*vm),
(a1**5*kc**2*(5*a1 + 5*a2 - 6*b2))/(ka*vm),
( - 85*a1**6*a2*kc + 85*a1**6*b2*kc - 13631488*ka*vm)/(85*ka*vm),
(6*(85*a1**6*a2*kc - 85*a1**6*b2*kc + 13631488*ka*vm))/(85*ka*vm),
(3*( - 85*a1**6*a2*kc + 85*a1**6*b2*kc - 13631488*ka*vm))/(17*ka*vm),
(4*(85*a1**6*a2*kc - 85*a1**6*b2*kc + 13631488*ka*vm))/(17*ka*vm),
(3*( - 85*a1**6*a2*kc + 85*a1**6*b2*kc - 13631488*ka*vm))/(17*ka*vm),
(6*(85*a1**6*a2*kc - 85*a1**6*b2*kc + 13631488*ka*vm))/(85*ka*vm),
( - 85*a1**6*a2*kc + 85*a1**6*b2*kc - 13631488*ka*vm)/(85*ka*vm),
( - 85*a1**6*a2*kc**2 + 85*a1**6*b2*kc**2 - 177209344*ka*vm)/(85*ka*vm),
(2125*a1**2*kc**2 + 85*a1*a2*kc**2 - 510*a1*b2*kc**2 - 129792*ka*vm)/(85*ka*vm),
( - 850*a1**3*kc**2 - 85*a1**2*a2*kc**2 + 255*a1**2*b2*kc**2 + 389376*ka*vm)/(17*ka*vm),
(2*(425*a1**4*kc**2 + 85*a1**3*a2*kc**2 - 170*a1**3*b2*kc**2 - 1384448*ka*vm))/(17*ka*vm),
( - 425*a1**5*kc**2 - 170*a1**4*a2*kc**2 + 255*a1**4*b2*kc**2 + 8306688*ka*vm)/(17*ka*vm),
(510*a1**6*a2*kc**2 - 510*a1**6*b2*kc**2 - 85*a1**5*b2*kc*vm + 1105854464*ka*vm)/(85*ka*vm),
(85*a1**5*kc**2 + 340*a1**3*a2*ka*vm + 85*a1**3*b1*ka*vm + 85*a1**3*b2*ka*vm - 9324032*ka*vm)/(85*ka*vm),
(170*a1**2*kc**2 - 17*a1*ka*vm + 17*b1*ka*vm + 17*b2*ka*vm - 21887*ka*vm)/(17*ka*vm),
(85*a1*kc - 85*b1*ka - 85*b2*ka + 85*b2*kc + 27*ka*vm)/(85*ka*vm),
(510*a1**4*kc**2 - 119*a1**3*ka*vm + 136*a1**2*a2*ka*vm + 51*a1**2*b1*ka*vm + 51*a1**2*b2*ka*vm - 4306752*ka*vm)/(17*ka*
vm),
(10*(17*a1**4*kc - 17*a1**3*a2*kc - 17*a1**3*b1*ka - 17*a1**3*b2*ka + 34*a1**3*b2*kc - 2560*ka*vm))/(17*ka*vm),
(30*( - 17*a1**4*kc + 17*a1**3*a2*kc + 17*a1**3*b1*ka + 17*a1**3*b2*ka - 34*a1**3*b2*kc + 2560*ka*vm))/(17*ka*vm),
(30*(17*a1**4*kc - 17*a1**3*a2*kc - 17*a1**3*b1*ka - 17*a1**3*b2*ka + 34*a1**3*b2*kc - 2560*ka*vm))/(17*ka*vm),
(10*( - 17*a1**4*kc + 17*a1**3*a2*kc + 17*a1**3*b1*ka + 17*a1**3*b2*ka - 34*a1**3*b2*kc + 2560*ka*vm))/(17*ka*vm),
(170*a1**3*kc - 85*a1**2*a2*kc - 170*a1**2*b1*ka - 170*a1**2*b2*ka + 255*a1**2*b2*kc + 128*ka*vm)/(17*ka*vm),
( - 170*a1**5*kc**2 + 51*a1**4*ka*vm - 68*a1**3*a2*ka*vm - 17*a1**3*b1*ka*vm - 17*a1**3*b2*ka*vm + 11623936*ka*vm)/(17*
ka*vm),
(85*a1**5*kc - 170*a1**4*a2*kc - 85*a1**4*b1*ka - 85*a1**4*b2*ka + 255*a1**4*b2*kc - 315392*ka*vm)/(17*ka*vm),
(4*( - 85*a1**5*kc + 170*a1**4*a2*kc + 85*a1**4*b1*ka + 85*a1**4*b2*ka - 255*a1**4*b2*kc + 315392*ka*vm))/(17*ka*vm),
(6*(85*a1**5*kc - 170*a1**4*a2*kc - 85*a1**4*b1*ka - 85*a1**4*b2*ka + 255*a1**4*b2*kc - 315392*ka*vm))/(17*ka*vm),
(4*( - 85*a1**5*kc + 170*a1**4*a2*kc + 85*a1**4*b1*ka + 85*a1**4*b2*ka - 255*a1**4*b2*kc + 315392*ka*vm))/(17*ka*vm),
(85*a1**5*kc - 170*a1**4*a2*kc - 85*a1**4*b1*ka - 85*a1**4*b2*ka + 255*a1**4*b2*kc - 315392*ka*vm)/(17*ka*vm),
(2*( - 170*a1**3*kc + 85*a1**2*a2*kc + 170*a1**2*b1*ka + 170*a1**2*b2*ka - 255*a1**2*b2*kc - 128*ka*vm))/(17*ka*vm),
( - 510*a1**3*kc**2 + 85*a1**2*ka*vm - 68*a1*a2*ka*vm - 51*a1*b1*ka*vm - 51*a1*b2*ka*vm + 531816*ka*vm)/(17*ka*vm),
(85*a1**6*kc - 425*a1**5*a2*kc - 85*a1**5*b1*ka - 85*a1**5*b2*ka + 510*a1**5*b2*kc - 7634944*ka*vm)/(85*ka*vm),
( - 85*a1**6*kc + 425*a1**5*a2*kc + 85*a1**5*b1*ka + 85*a1**5*b2*ka - 510*a1**5*b2*kc + 7634944*ka*vm)/(17*ka*vm),
(2*(85*a1**6*kc - 425*a1**5*a2*kc - 85*a1**5*b1*ka - 85*a1**5*b2*ka + 510*a1**5*b2*kc - 7634944*ka*vm))/(17*ka*vm),
(2*( - 85*a1**6*kc + 425*a1**5*a2*kc + 85*a1**5*b1*ka + 85*a1**5*b2*ka - 510*a1**5*b2*kc + 7634944*ka*vm))/(17*ka*vm),
(85*a1**6*kc - 425*a1**5*a2*kc - 85*a1**5*b1*ka - 85*a1**5*b2*ka + 510*a1**5*b2*kc - 7634944*ka*vm)/(17*ka*vm),
( - 85*a1**6*kc + 425*a1**5*a2*kc + 85*a1**5*b1*ka + 85*a1**5*b2*ka - 510*a1**5*b2*kc + 7634944*ka*vm)/(85*ka*vm),
(170*a1**3*kc - 85*a1**2*a2*kc - 170*a1**2*b1*ka - 170*a1**2*b2*ka + 255*a1**2*b2*kc + 128*ka*vm)/(17*ka*vm),
(170*a1**3*kc**2 - 85*a1**2*ka*vm + 68*a1*a2*ka*vm + 51*a1*b1*ka*vm + 51*a1*b2*ka*vm - 185704*ka*vm)/(17*ka*vm),
(2*( - 170*a1**4*kc**2 + 68*a1**3*ka*vm - 136*a1**2*a2*ka*vm - 51*a1**2*b1*ka*vm - 51*a1**2*b2*ka*vm + 1563968*ka*vm))/(
17*ka*vm),
(170*a1**5*kc**2 - 51*a1**4*ka*vm + 204*a1**3*a2*ka*vm + 51*a1**3*b1*ka*vm + 51*a1**3*b2*ka*vm - 13138432*ka*vm)/(17*ka*
vm),
( - 85*a1**5*kc**2 + 17*a1**4*ka*vm - 204*a1**3*a2*ka*vm - 51*a1**3*b1*ka*vm - 51*a1**3*b2*ka*vm + 7739904*ka*vm)/(17*ka
*vm),
(85*a1**4*kc**2 - 51*a1**3*ka*vm + 136*a1**2*a2*ka*vm + 51*a1**2*b1*ka*vm + 51*a1**2*b2*ka*vm - 880448*ka*vm)/(17*ka*vm)
,
(425*a1**2*kc - 85*a1*a2*kc - 425*a1*b1*ka - 425*a1*b2*ka + 510*a1*b2*kc + 664*ka*vm)/(85*ka*vm),
( - 425*a1**2*kc + 85*a1*a2*kc + 425*a1*b1*ka + 425*a1*b2*ka - 510*a1*b2*kc - 664*ka*vm)/(85*ka*vm),
( - 85*a1**6*a2*kc**2 + 85*a1**6*b2*kc**2 + 85*a1**5*b2*kc*vm - 85*a1**4*a2**2*ka*vm - 85*a1**4*a2*b1*ka*vm - 85*a1**4*
a2*b2*ka*vm - 2555904*ka*vm)/(85*ka*vm),
(510*a1**6*a2*kc**2 - 510*a1**6*b2*kc**2 - 425*a1**5*b2*kc*vm + 340*a1**4*a2**2*ka*vm + 340*a1**4*a2*b1*ka*vm + 340*a1**
4*a2*b2*ka*vm + 407240704*ka*vm)/(85*ka*vm),
( - 255*a1**6*a2*kc**2 + 255*a1**6*b2*kc**2 + 170*a1**5*b2*kc*vm - 102*a1**4*a2**2*ka*vm - 102*a1**4*a2*b1*ka*vm - 102*
a1**4*a2*b2*ka*vm - 356122624*ka*vm)/(17*ka*vm),
(2*(170*a1**6*a2*kc**2 - 170*a1**6*b2*kc**2 - 85*a1**5*b2*kc*vm + 34*a1**4*a2**2*ka*vm + 34*a1**4*a2*b1*ka*vm + 34*a1**4
*a2*b2*ka*vm + 310116352*ka*vm))/(17*ka*vm),
( - 255*a1**6*a2*kc**2 + 255*a1**6*b2*kc**2 + 85*a1**5*b2*kc*vm - 17*a1**4*a2**2*ka*vm - 17*a1**4*a2*b1*ka*vm - 17*a1**4
*a2*b2*ka*vm - 530776064*ka*vm)/(17*ka*vm),
( - 85*a1**5*a2*kc**2 + 102*a1**5*b2*kc**2 + 17*a1**4*a2*ka*vm + 85*a1**4*b2*kc*vm - 17*a1**3*a2**2*ka*vm - 51*a1**3*a2*
b1*ka*vm - 51*a1**3*a2*b2*ka*vm - 38903808*ka*vm)/(17*ka*vm),
( - 170*a1**6*kc**2 - 425*a1**5*a2*kc**2 + 510*a1**5*b2*kc**2 - 34*a1**4*a2*ka*vm - 17*a1**4*b1*ka*vm - 17*a1**4*b2*ka*
vm + 85*a1**4*b2*kc*vm - 134402048*ka*vm)/(17*ka*vm),
( - 850*a1**2*kc**2 - 85*a1*a2*kc**2 + 510*a1*b2*kc**2 + 170*a1*ka*vm - 85*b1*ka*vm - 85*b2*ka*vm + 85*b2*kc*vm - 33153*
ka*vm)/(85*ka*vm),
(2*(340*a1**5*kc**2 + 340*a1**4*a2*kc**2 - 510*a1**4*b2*kc**2 - 17*a1**4*ka*vm + 51*a1**3*a2*ka*vm + 34*a1**3*b1*ka*vm +
 34*a1**3*b2*ka*vm - 85*a1**3*b2*kc*vm + 16530432*ka*vm))/(17*ka*vm),
(2*( - 510*a1**4*kc**2 - 255*a1**3*a2*kc**2 + 510*a1**3*b2*kc**2 + 51*a1**3*ka*vm - 51*a1**2*a2*ka*vm - 51*a1**2*b1*ka*
vm - 51*a1**2*b2*ka*vm + 85*a1**2*b2*kc*vm - 2058688*ka*vm))/(17*ka*vm),
( - 170*a1**3*kc**2 - 85*a1**2*a2*kc**2 + 255*a1**2*b2*kc**2 + 68*a1**2*ka*vm - 51*a1*a2*ka*vm - 51*a1*b1*ka*vm - 51*a1*
b2*ka*vm + 85*a1*b2*kc*vm - 301688*ka*vm)/(17*ka*vm),
(680*a1**3*kc**2 + 170*a1**2*a2*kc**2 - 510*a1**2*b2*kc**2 - 102*a1**2*ka*vm + 34*a1*a2*ka*vm + 68*a1*b1*ka*vm + 68*a1*
b2*ka*vm - 85*a1*b2*kc*vm + 259648*ka*vm)/(17*ka*vm),
( - 85*a1**5*kc**2 - 850*a1**4*a2*kc**2 + 1275*a1**4*b2*kc**2 - 85*a1**3*a2*ka*vm - 85*a1**3*b1*ka*vm - 85*a1**3*b2*ka*
vm + 850*a1**3*b2*kc*vm + 85*a1**2*a2**2*ka*vm - 255*a1**2*a2*b1*ka*vm - 255*a1**2*a2*b2*ka*vm - 91719168*ka*vm)/(85*ka*
vm),
(85*a1**6*kc**2 + 2125*a1**5*a2*kc**2 - 2550*a1**5*b2*kc**2 - 85*a1**4*a2*ka*vm + 85*a1**4*b1*ka*vm + 85*a1**4*b2*ka*vm 
- 1700*a1**4*b2*kc*vm + 255*a1**3*a2**2*ka*vm + 765*a1**3*a2*b1*ka*vm + 765*a1**3*a2*b2*ka*vm + 989753344*ka*vm)/(85*ka*
vm),
( - 85*a1**6*kc**2 - 850*a1**5*a2*kc**2 + 1020*a1**5*b2*kc**2 - 51*a1**4*a2*ka*vm - 51*a1**4*b1*ka*vm - 51*a1**4*b2*ka*
vm + 510*a1**4*b2*kc*vm - 51*a1**3*a2**2*ka*vm - 153*a1**3*a2*b1*ka*vm - 153*a1**3*a2*b2*ka*vm - 387010560*ka*vm)/(17*ka
*vm),
(170*a1**6*kc**2 + 850*a1**5*a2*kc**2 - 1020*a1**5*b2*kc**2 + 85*a1**4*a2*ka*vm + 51*a1**4*b1*ka*vm + 51*a1**4*b2*ka*vm 
- 340*a1**4*b2*kc*vm + 17*a1**3*a2**2*ka*vm + 51*a1**3*a2*b1*ka*vm + 51*a1**3*a2*b2*ka*vm + 353505280*ka*vm)/(17*ka*vm),
( - 85*a1**4*kc**2 - 170*a1**3*a2*kc**2 + 340*a1**3*b2*kc**2 + 34*a1**3*ka*vm - 85*a1**2*a2*ka*vm - 51*a1**2*b1*ka*vm - 
51*a1**2*b2*ka*vm + 170*a1**2*b2*kc*vm + 17*a1*a2**2*ka*vm - 17*a1*a2*b1*ka*vm - 17*a1*a2*b2*ka*vm - 3467584*ka*vm)/(17*
ka*vm),
(2*(85*a1**5*kc**2 + 340*a1**4*a2*kc**2 - 510*a1**4*b2*kc**2 - 17*a1**4*ka*vm + 68*a1**3*a2*ka*vm + 51*a1**3*b1*ka*vm + 
51*a1**3*b2*ka*vm - 255*a1**3*b2*kc*vm - 17*a1**2*a2**2*ka*vm + 51*a1**2*a2*b1*ka*vm + 51*a1**2*a2*b2*ka*vm + 33212928*
ka*vm))/(17*ka*vm),
(510*a1**4*kc**2 + 510*a1**3*a2*kc**2 - 1020*a1**3*b2*kc**2 - 136*a1**3*ka*vm + 187*a1**2*a2*ka*vm + 153*a1**2*b1*ka*vm 
+ 153*a1**2*b2*ka*vm - 340*a1**2*b2*kc*vm - 17*a1*a2**2*ka*vm + 17*a1*a2*b1*ka*vm + 17*a1*a2*b2*ka*vm + 8277184*ka*vm)/(
17*ka*vm),
( - 510*a1**5*kc**2 - 1020*a1**4*a2*kc**2 + 1530*a1**4*b2*kc**2 + 68*a1**4*ka*vm - 221*a1**3*a2*ka*vm - 153*a1**3*b1*ka*
vm - 153*a1**3*b2*ka*vm + 510*a1**3*b2*kc*vm + 17*a1**2*a2**2*ka*vm - 51*a1**2*a2*b1*ka*vm - 51*a1**2*a2*b2*ka*vm - 
82804224*ka*vm)/(17*ka*vm)}$

MODEL PARAMETER SOLUTION(S)$

 G_:=GROESOLVE(FLIST_,B1_) $

g_ := {{b1=3,b2=20,vm=5,ka=17,a2=16,kc=13,a1=8}}$

MODEL GLOBALLY IDENTIFIABLE$
Elapsed time for Pharmacokinetics: 3.0434571 seconds
