
NFkB

seed_ := 65$

NUMBER OF EQUATIONS$

n_ := 22$

VARIABLES VECTOR$

b_ := {x1,
x2,
x3,
x4,
x5,
x6,
x7,
x8,
x9,
x10,
x11,
x12,
x13,
x14,
x15,
u1,
y1,
y2,
y3,
y4,
y5,
y6}$

UNKNOWN PARAMETER(S) VECTOR$

b1_ := {t1,
a3,
t2,
c1a,
c2a,
c3a,
c4a,
c5a,
c6a,
c1,
c2,
c3,
c4,
c5,
k1,
k3,
kprod,
kdeg,
kv,
i1,
e2a,
i1a,
e1a,
c1c,
c2c,
c3c}$

RANKING AMONG THE VARIABLES$

bb_ := {x1,
x2,
x3,
x4,
x5,
x6,
x7,
df(x1,t),
df(x2,t),
df(x3,t),
df(x4,t),
df(x5,t),
df(x6,t),
df(x7,t),
df(x1,t,2),
df(x2,t,2),
df(x3,t,2),
df(x4,t,2),
df(x5,t,2),
df(x6,t,2),
df(x7,t,2),
df(x1,t,3),
df(x2,t,3),
df(x3,t,3),
df(x4,t,3),
df(x5,t,3),
df(x6,t,3),
df(x7,t,3),
df(x1,t,4),
df(x2,t,4),
df(x3,t,4),
df(x4,t,4),
df(x5,t,4),
df(x6,t,4),
df(x7,t,4),
df(x1,t,5),
df(x2,t,5),
df(x3,t,5),
df(x4,t,5),
df(x5,t,5),
df(x6,t,5),
df(x7,t,5),
df(x1,t,6),
df(x2,t,6),
df(x3,t,6),
df(x4,t,6),
df(x5,t,6),
df(x6,t,6),
df(x7,t,6),
df(x1,t,7),
df(x2,t,7),
df(x3,t,7),
df(x4,t,7),
df(x5,t,7),
df(x6,t,7),
df(x7,t,7),
df(x1,t,8),
df(x2,t,8),
df(x3,t,8),
df(x4,t,8),
df(x5,t,8),
df(x6,t,8),
df(x7,t,8),
df(x1,t,9),
df(x2,t,9),
df(x3,t,9),
df(x4,t,9),
df(x5,t,9),
df(x6,t,9),
df(x7,t,9),
df(x1,t,10),
df(x2,t,10),
df(x3,t,10),
df(x4,t,10),
df(x5,t,10),
df(x6,t,10),
df(x7,t,10),
df(x1,t,11),
df(x2,t,11),
df(x3,t,11),
df(x4,t,11),
df(x5,t,11),
df(x6,t,11),
df(x7,t,11),
df(x1,t,12),
df(x2,t,12),
df(x3,t,12),
df(x4,t,12),
df(x5,t,12),
df(x6,t,12),
df(x7,t,12),
df(x1,t,13),
df(x2,t,13),
df(x3,t,13),
df(x4,t,13),
df(x5,t,13),
df(x6,t,13),
df(x7,t,13),
df(x1,t,14),
df(x2,t,14),
df(x3,t,14),
df(x4,t,14),
df(x5,t,14),
df(x6,t,14),
df(x7,t,14),
df(x1,t,15),
df(x2,t,15),
df(x3,t,15),
df(x4,t,15),
df(x5,t,15),
df(x6,t,15),
df(x7,t,15),
x8,
x9,
x10,
x11,
x12,
x13,
x14,
x15,
u1,
y1,
y2,
y3,
y4,
y5,
y6,
df(x8,t),
df(x9,t),
df(x10,t),
df(x11,t),
df(x12,t),
df(x13,t),
df(x14,t),
df(x15,t),
df(u1,t),
df(y1,t),
df(y2,t),
df(y3,t),
df(y4,t),
df(y5,t),
df(y6,t)}$

NUMBER OF INPUT(S)$

nu_ := 1$

NUMBER OF OUTPUT(S)$

ny_ := 6$

NUMBER OF STATE(S) $

nx_ := 15$

MODEL EQUATION(S)$

c_ := {df(x1,t)=kprod + u1 - kdeg*x1,
df(x2,t)=t1*x4 + t2*x5 - kdeg*x2 - k3*x2 - a3*x13*x2 - a2*x10*x2,
df(x3,t)=k3*x2 - kdeg*x3,
df(x4,t)=a2*x10*x2 - t1*x4,
df(x5,t)=a3*x13*x2 - t2*x5,
df(x6,t)= - (i1*x6 - t2*x5 - c6a*x13) - a1*x10*x6,
df(x7,t)= - a1*x11*x7 + i1*kv*x6,
df(x8,t)=c4*x9 - c5*x8,
df(x9,t)=c2 - c3*x9 + c1*x7,
df(x10,t)=e1a*x11 - i1a*x10 - c5a*x10 + c4a*x12 - a2*x10*x2 - a1*x10*x6,
df(x11,t)= - (e1a*x11 - i1a*x10)*kv - a1*x11*x7,
df(x12,t)=c2a - c3a*x12 + c1a*x7,
df(x13,t)= - (c6a*x13 - e2a*x14 + a3*x13*x2) + a1*x10*x6,
df(x14,t)=a1*x11*x7 - e2a*kv*x14,
df(x15,t)=c2c - c3c*x15 + c1c*x7,
 - k1*x1=u1,
x7=y1,
x10 + x13=y2,
x9=y3,
x2 + x3 + x1=y4,
x2=y5,
x12=y6}$

CHARACTERISTIC SET$

aa_(1) :=  - df(x1,t) - x1*(k1 + kdeg) + kprod$

aa_(2) := df(x3,t) - x2*k3 + x3*kdeg$

aa_(3) := df(x2,t) + df(x4,t) + df(x5,t) + x2*(k3 + kdeg)$

aa_(4) := df(x2,t)*a2*c6a + df(x4,t)*x6*a1*a3 + df(x4,t)*a2*c6a + df(x6,t)*x2*a2*a3 - x2*x5*a2*a3*t2 + x2*x6*a2*a3*i1 + 
x2*a2*c6a*(k3 + kdeg) + x4*x6*a1*a3*t1 - x5*a2*c6a*t2$

aa_(5) := df(x2,t)*x7*a2*c6a*i1*kv + df(x4,t)*x6*x7*a1*a3*i1*kv + df(x4,t)*x7**2*a1*a3*i1a*kv + df(x4,t)*x7*a2*c6a*i1*kv
 + df(x7,t,2)*x2*x7*a2*a3 - df(x7,t)**2*x2*a2*a3 + df(x7,t)*x2*x6*a2*a3*i1*kv + df(x7,t)*x2*x7**2*a1*a2*a3 + df(x7,t)*x2
*x7*a2*a3*e1a*kv - x2*x5*x7*a2*a3*i1*kv*t2 - x2*x6*x7**2*a1*a2*a3*i1*kv + x2*x6*x7*a2*a3*i1*kv*( - e1a*kv + i1) + x2*x7*
a2*c6a*i1*kv*(k3 + kdeg) + x4*x6*x7*a1*a3*i1*kv*t1 + x4*x7**2*a1*a3*i1a*kv*t1 - x5*x7*a2*c6a*i1*kv*t2$

aa_(6) :=  - df(x2,t,3)*x2**2*x7*a1*a2**2*a3 + 3*df(x2,t,2)*df(x2,t)*x2*x7*a1*a2**2*a3 - df(x2,t,2)*x2**3*x7*a1*a2**2*a3
**2 - df(x2,t,2)*x2**2*x7*a1*a2**2*a3*(c6a + e2a*kv + t2) - df(x2,t,2)*x2*x4*x7*a1*a2**2*a3*t1 - df(x2,t,2)*x2*x5*x7*a1*
a2**2*a3*t2 - 2*df(x2,t)**3*x7*a1*a2**2*a3 + df(x2,t)**2*x2*x7*a1*a2**2*a3*(c6a + e2a*kv + 2*t2) + 2*df(x2,t)**2*x4*x7*
a1*a2**2*a3*t1 + 2*df(x2,t)**2*x5*x7*a1*a2**2*a3*t2 + df(x2,t)*df(x4,t)*x2*x6*x7*a1**2*a2*a3*( - a2 + a3) + df(x2,t)*df(
x4,t)*x2*x7*a1*a2*( - a1*a2*c6a + a1*a3*c6a - a2*a3*c3a - a2*a3*c5a + a2*a3*c6a + a2*a3*e2a*kv - a2*a3*i1a - 2*a2*a3*t1 
+ 2*a2*a3*t2) - df(x2,t)*x2**3*x7*a1*a2**2*a3**2*(e2a*kv + k3 + kdeg + t2) + df(x2,t)*x2**2*x7*a1*a2**2*a3*( - c6a*e2a*
kv - c6a*t2 - e2a*kv*t2 + k3*t2 + kdeg*t2) + df(x2,t)*x2*x4*x6*x7*a1**2*a2*a3*t1*( - a2 + a3) + df(x2,t)*x2*x4*x7*a1*a2*
t1*( - a1*a2*c6a + a1*a3*c6a - a2*a3*c3a - a2*a3*c5a - a2*a3*i1a) - df(x2,t)*x2*x5*x7*a1*a2**2*a3*t2*(c6a + e2a*kv) + df
(x4,t,2)*x2**3*x7*a1*a2**2*a3*(a2 - a3) + df(x4,t,2)*x2**2*x6*x7*a1**2*a2*a3*(a2 - a3) + df(x4,t,2)*x2**2*x7*a1*a2**2*a3
*(c3a + c5a - c6a - e2a*kv + i1a + t1 - t2) + df(x4,t)**2*x2*x6*x7*a1**3*a3*( - a2 + a3) + df(x4,t)**2*x2*x7*a1**2*a2*
c6a*( - a2 + a3) + df(x4,t)*x2**3*x7*a1*a2**2*a3*(a2*c3a + a2*t1 - a3*e2a*kv - a3*t2) + df(x4,t)*x2**2*x5*x7*a1**2*a2*a3
*t2*(a2 - a3) + df(x4,t)*x2**2*x6*x7*a1**2*a2*a3*(a2*c3a - a2*i1 + a2*t1 - a3*e2a*kv + a3*i1 - a3*t1) + df(x4,t)*x2**2*
x7*a1*a2*( - a1*a2*c6a*k3 - a1*a2*c6a*kdeg + a1*a3*c6a*k3 + a1*a3*c6a*kdeg + a2*a3*c3a*c5a + a2*a3*c3a*i1a + a2*a3*c3a*
t1 + a2*a3*c5a*t1 - a2*a3*c6a*e2a*kv - a2*a3*c6a*t2 - a2*a3*e1a*i1a*kv - a2*a3*e2a*kv*t2 + a2*a3*i1a*t1) + 2*df(x4,t)*x2
*x4*x6*x7*a1**3*a3*t1*( - a2 + a3) + df(x4,t)*x2*x4*x7*a1**2*a2*c6a*t1*( - a2 + a3) + df(x4,t)*x2*x5*x7*a1**2*a2*c6a*t2*
(a2 - a3) + df(x7,t)*x2**3*x7*a1*a2**2*a3*( - a2*e1a + a3*e2a) + df(x7,t)*x2**3*a2**3*a3*e1a*(c3a - e1a*kv) - x2**4*x7*
a1*a2**2*a3**2*(e2a*k3*kv + e2a*kdeg*kv + k3*t2 + kdeg*t2) + x2**3*x4*x7*a1*a2**3*a3*c3a*t1 + x2**3*x5*x7*a1*a2**2*a3**2
*e2a*kv*t2 + x2**3*x6*x7*a1*a2**2*a3*i1*kv*(a2*e1a - a3*e2a) + x2**3*x6*a2**3*a3*e1a*i1*kv*( - c3a + e1a*kv) - x2**3*x7
**2*a1*a2**3*a3*c1a*c4a - x2**3*x7*a1*a2**2*a3*(a2*c2a*c4a + c6a*e2a*k3*kv + c6a*e2a*kdeg*kv + c6a*k3*t2 + c6a*kdeg*t2 +
 e2a*k3*kv*t2 + e2a*kdeg*kv*t2) + x2**2*x4*x5*x7*a1**2*a2*a3*t1*t2*(a2 - a3) + x2**2*x4*x6*x7*a1**2*a2*a3*t1*(a2*c3a - 
a2*i1 - a3*e2a*kv + a3*i1) + x2**2*x4*x7*a1*a2*t1*( - a1*a2*c6a*k3 - a1*a2*c6a*kdeg + a1*a3*c6a*k3 + a1*a3*c6a*kdeg + a2
*a3*c3a*c5a + a2*a3*c3a*i1a - a2*a3*e1a*i1a*kv) + x2**2*x5*x7*a1*a2**2*a3*c6a*e2a*kv*t2 + x2*x4**2*x6*x7*a1**3*a3*t1**2*
( - a2 + a3) + x2*x4*x5*x7*a1**2*a2*c6a*t1*t2*(a2 - a3)$

aa_(7) :=  - df(x2,t,2)*df(x4,t)*x2*x7*a1*a2*a3 - df(x2,t,2)*x2*x4*x7*a1*a2*a3*t1 + 2*df(x2,t)**2*df(x4,t)*x7*a1*a2*a3 +
 2*df(x2,t)**2*x4*x7*a1*a2*a3*t1 - 2*df(x2,t)*df(x4,t,2)*x2*x7*a1*a2*a3 - df(x2,t)*df(x4,t)*x2*x6*x7*a1**2*a2*a3 + df(x2
,t)*df(x4,t)*x2*x7*a1*a2*( - a1*c6a - a3*c3a - a3*c5a - a3*i1a - 2*a3*t1) - df(x2,t)*x2*x4*x6*x7*a1**2*a2*a3*t1 - df(x2,
t)*x2*x4*x7*a1*a2*t1*(a1*c6a + a3*c3a + a3*c5a + a3*i1a) + df(x4,t,3)*x2**2*x7*a1*a2*a3 + df(x4,t,2)*x2**3*x7*a1*a2**2*
a3 + df(x4,t,2)*x2**2*x6*x7*a1**2*a2*a3 + df(x4,t,2)*x2**2*x7*a1*a2*a3*(c3a + c5a + i1a + t1) - df(x4,t)**2*x2*x6*x7*a1
**3*a3 - df(x4,t)**2*x2*x7*a1**2*a2*c6a + df(x4,t)*x2**3*x7*a1*a2**2*a3*(c3a + t1) + df(x4,t)*x2**2*x5*x7*a1**2*a2*a3*t2
 + df(x4,t)*x2**2*x6*x7*a1**2*a2*a3*(c3a - i1 + t1) + df(x4,t)*x2**2*x7*a1*a2*( - a1*c6a*k3 - a1*c6a*kdeg + a3*c3a*c5a +
 a3*c3a*i1a + a3*c3a*t1 + a3*c5a*t1 - a3*e1a*i1a*kv + a3*i1a*t1) - 2*df(x4,t)*x2*x4*x6*x7*a1**3*a3*t1 - df(x4,t)*x2*x4*
x7*a1**2*a2*c6a*t1 + df(x4,t)*x2*x5*x7*a1**2*a2*c6a*t2 - df(x7,t)*x2**3*x7*a1*a2**2*a3*e1a + df(x7,t)*x2**3*a2**2*a3*e1a
*(c3a - e1a*kv) + x2**3*x4*x7*a1*a2**2*a3*c3a*t1 + x2**3*x6*x7*a1*a2**2*a3*e1a*i1*kv + x2**3*x6*a2**2*a3*e1a*i1*kv*( - 
c3a + e1a*kv) - x2**3*x7**2*a1*a2**2*a3*c1a*c4a - x2**3*x7*a1*a2**2*a3*c2a*c4a + x2**2*x4*x5*x7*a1**2*a2*a3*t1*t2 + x2**
2*x4*x6*x7*a1**2*a2*a3*t1*(c3a - i1) + x2**2*x4*x7*a1*a2*t1*( - a1*c6a*k3 - a1*c6a*kdeg + a3*c3a*c5a + a3*c3a*i1a - a3*
e1a*i1a*kv) - x2*x4**2*x6*x7*a1**3*a3*t1**2 + x2*x4*x5*x7*a1**2*a2*c6a*t1*t2$

aa_(8) := df(x4,t) - x10*x2*a2 + x4*t1$

aa_(9) := df(x7,t) + x11*x7*a1 - x6*i1*kv$

aa_(10) :=  - df(x2,t)*df(x4,t)*x7*a1 - df(x2,t)*x4*x7*a1*t1 + df(x4,t,2)*x2*x7*a1 + df(x4,t)*x2**2*x7*a1*a2 + df(x4,t)*
x2*x6*x7*a1**2 + df(x4,t)*x2*x7*a1*(c5a + i1a + t1) + df(x7,t)*x2**2*a2*e1a - x12*x2**2*x7*a1*a2*c4a + x2**2*x4*x7*a1*a2
*t1 - x2**2*x6*a2*e1a*i1*kv + x2*x4*x6*x7*a1**2*t1 + x2*x4*x7*a1*t1*(c5a + i1a)$

aa_(11) :=  - df(x2,t) - df(x4,t) - x13*x2*a3 - x2*(k3 + kdeg) + x5*t2$

aa_(12) := df(x2,t,2)*x2*a2 - df(x2,t)**2*a2 - df(x2,t)*df(x4,t)*a2 + df(x2,t)*x2**2*a2*a3 + df(x2,t)*x2*a2*(c6a + t2) +
 df(x2,t)*x5*a2*t2 + df(x4,t,2)*x2*a2 + df(x4,t)*x2**2*a2*a3 + df(x4,t)*x2*x6*a1*a3 + df(x4,t)*x2*a2*(c6a + t2) + x14*x2
**2*a2*a3*e2a + x2**3*a2*a3*(k3 + kdeg) - x2**2*x5*a2*a3*t2 + x2**2*a2*(c6a*k3 + c6a*kdeg + k3*t2 + kdeg*t2) + x2*x4*x6*
a1*a3*t1 - x2*x5*a2*c6a*t2$

aa_(13) :=  - u1 - x1*k1$

aa_(14) := x7 - y1$

aa_(15) :=  - df(x2,t)*a2 + df(x4,t)*( - a2 + a3) - x2*y2*a2*a3 - x2*a2*(k3 + kdeg) + x4*a3*t1 + x5*a2*t2$

aa_(16) := x9 - y3$

aa_(17) := x1 + x2 + x3 - y4$

aa_(18) := x2 - y5$

aa_(19) := df(x2,t)*df(x4,t)*x7*a1 + df(x2,t)*x4*x7*a1*t1 - df(x4,t,2)*x2*x7*a1 - df(x4,t)*x2**2*x7*a1*a2 - df(x4,t)*x2*
x6*x7*a1**2 - df(x4,t)*x2*x7*a1*(c5a + i1a + t1) - df(x7,t)*x2**2*a2*e1a - x2**2*x4*x7*a1*a2*t1 + x2**2*x6*a2*e1a*i1*kv 
+ x2**2*x7*y6*a1*a2*c4a - x2*x4*x6*x7*a1**2*t1 - x2*x4*x7*a1*t1*(c5a + i1a)$

aa_(20) := df(x8,t) + x8*c5 - x9*c4$

aa_(21) := df(x9,t) - x7*c1 + x9*c3 - c2$

aa_(22) := df(x15,t) + x15*c3c - x7*c1c - c2c$

MODEL NOT ALGEBRAICALLY OBSERVABLE$

NORMALIZED  INPUT /OUTPUT RELATION(S) $

aan_(1) :=  - df(x1,t) - x1*(k1 + kdeg) + kprod$

aan_(2) := df(x3,t) - x2*k3 + x3*kdeg$

aan_(3) := df(x2,t) + df(x4,t) + df(x5,t) + x2*(k3 + kdeg)$

aan_(4) := (df(x2,t)*a2*c6a + df(x4,t)*x6*a1*a3 + df(x4,t)*a2*c6a + df(x6,t)*x2*a2*a3 - x2*x5*a2*a3*t2 + x2*x6*a2*a3*i1 
+ x2*a2*c6a*(k3 + kdeg) + x4*x6*a1*a3*t1 - x5*a2*c6a*t2)/(a2*c6a)$

aan_(5) := (df(x2,t)*x7*a2*c6a*i1*kv + df(x4,t)*x6*x7*a1*a3*i1*kv + df(x4,t)*x7**2*a1*a3*i1a*kv + df(x4,t)*x7*a2*c6a*i1*
kv + df(x7,t,2)*x2*x7*a2*a3 - df(x7,t)**2*x2*a2*a3 + df(x7,t)*x2*x6*a2*a3*i1*kv + df(x7,t)*x2*x7**2*a1*a2*a3 + df(x7,t)*
x2*x7*a2*a3*e1a*kv - x2*x5*x7*a2*a3*i1*kv*t2 - x2*x6*x7**2*a1*a2*a3*i1*kv + x2*x6*x7*a2*a3*i1*kv*( - e1a*kv + i1) + x2*
x7*a2*c6a*i1*kv*(k3 + kdeg) + x4*x6*x7*a1*a3*i1*kv*t1 + x4*x7**2*a1*a3*i1a*kv*t1 - x5*x7*a2*c6a*i1*kv*t2)/(a2*c6a*i1*kv)
$

aan_(6) := (df(x2,t,3)*x2**2*x7*a1*a2**2*a3 - 3*df(x2,t,2)*df(x2,t)*x2*x7*a1*a2**2*a3 + df(x2,t,2)*x2**3*x7*a1*a2**2*a3
**2 + df(x2,t,2)*x2**2*x7*a1*a2**2*a3*(c6a + e2a*kv + t2) + df(x2,t,2)*x2*x4*x7*a1*a2**2*a3*t1 + df(x2,t,2)*x2*x5*x7*a1*
a2**2*a3*t2 + 2*df(x2,t)**3*x7*a1*a2**2*a3 + df(x2,t)**2*x2*x7*a1*a2**2*a3*( - c6a - e2a*kv - 2*t2) - 2*df(x2,t)**2*x4*
x7*a1*a2**2*a3*t1 - 2*df(x2,t)**2*x5*x7*a1*a2**2*a3*t2 + df(x2,t)*df(x4,t)*x2*x6*x7*a1**2*a2*a3*(a2 - a3) + df(x2,t)*df(
x4,t)*x2*x7*a1*a2*(a1*a2*c6a - a1*a3*c6a + a2*a3*c3a + a2*a3*c5a - a2*a3*c6a - a2*a3*e2a*kv + a2*a3*i1a + 2*a2*a3*t1 - 2
*a2*a3*t2) + df(x2,t)*x2**3*x7*a1*a2**2*a3**2*(e2a*kv + k3 + kdeg + t2) + df(x2,t)*x2**2*x7*a1*a2**2*a3*(c6a*e2a*kv + 
c6a*t2 + e2a*kv*t2 - k3*t2 - kdeg*t2) + df(x2,t)*x2*x4*x6*x7*a1**2*a2*a3*t1*(a2 - a3) + df(x2,t)*x2*x4*x7*a1*a2*t1*(a1*
a2*c6a - a1*a3*c6a + a2*a3*c3a + a2*a3*c5a + a2*a3*i1a) + df(x2,t)*x2*x5*x7*a1*a2**2*a3*t2*(c6a + e2a*kv) + df(x4,t,2)*
x2**3*x7*a1*a2**2*a3*( - a2 + a3) + df(x4,t,2)*x2**2*x6*x7*a1**2*a2*a3*( - a2 + a3) + df(x4,t,2)*x2**2*x7*a1*a2**2*a3*( 
- c3a - c5a + c6a + e2a*kv - i1a - t1 + t2) + df(x4,t)**2*x2*x6*x7*a1**3*a3*(a2 - a3) + df(x4,t)**2*x2*x7*a1**2*a2*c6a*(
a2 - a3) + df(x4,t)*x2**3*x7*a1*a2**2*a3*( - a2*c3a - a2*t1 + a3*e2a*kv + a3*t2) + df(x4,t)*x2**2*x5*x7*a1**2*a2*a3*t2*(
 - a2 + a3) + df(x4,t)*x2**2*x6*x7*a1**2*a2*a3*( - a2*c3a + a2*i1 - a2*t1 + a3*e2a*kv - a3*i1 + a3*t1) + df(x4,t)*x2**2*
x7*a1*a2*(a1*a2*c6a*k3 + a1*a2*c6a*kdeg - a1*a3*c6a*k3 - a1*a3*c6a*kdeg - a2*a3*c3a*c5a - a2*a3*c3a*i1a - a2*a3*c3a*t1 -
 a2*a3*c5a*t1 + a2*a3*c6a*e2a*kv + a2*a3*c6a*t2 + a2*a3*e1a*i1a*kv + a2*a3*e2a*kv*t2 - a2*a3*i1a*t1) + 2*df(x4,t)*x2*x4*
x6*x7*a1**3*a3*t1*(a2 - a3) + df(x4,t)*x2*x4*x7*a1**2*a2*c6a*t1*(a2 - a3) + df(x4,t)*x2*x5*x7*a1**2*a2*c6a*t2*( - a2 + 
a3) + df(x7,t)*x2**3*x7*a1*a2**2*a3*(a2*e1a - a3*e2a) + df(x7,t)*x2**3*a2**3*a3*e1a*( - c3a + e1a*kv) + x2**4*x7*a1*a2**
2*a3**2*(e2a*k3*kv + e2a*kdeg*kv + k3*t2 + kdeg*t2) - x2**3*x4*x7*a1*a2**3*a3*c3a*t1 - x2**3*x5*x7*a1*a2**2*a3**2*e2a*kv
*t2 + x2**3*x6*x7*a1*a2**2*a3*i1*kv*( - a2*e1a + a3*e2a) + x2**3*x6*a2**3*a3*e1a*i1*kv*(c3a - e1a*kv) + x2**3*x7**2*a1*
a2**3*a3*c1a*c4a + x2**3*x7*a1*a2**2*a3*(a2*c2a*c4a + c6a*e2a*k3*kv + c6a*e2a*kdeg*kv + c6a*k3*t2 + c6a*kdeg*t2 + e2a*k3
*kv*t2 + e2a*kdeg*kv*t2) + x2**2*x4*x5*x7*a1**2*a2*a3*t1*t2*( - a2 + a3) + x2**2*x4*x6*x7*a1**2*a2*a3*t1*( - a2*c3a + a2
*i1 + a3*e2a*kv - a3*i1) + x2**2*x4*x7*a1*a2*t1*(a1*a2*c6a*k3 + a1*a2*c6a*kdeg - a1*a3*c6a*k3 - a1*a3*c6a*kdeg - a2*a3*
c3a*c5a - a2*a3*c3a*i1a + a2*a3*e1a*i1a*kv) - x2**2*x5*x7*a1*a2**2*a3*c6a*e2a*kv*t2 + x2*x4**2*x6*x7*a1**3*a3*t1**2*(a2 
- a3) + x2*x4*x5*x7*a1**2*a2*c6a*t1*t2*( - a2 + a3))/(a1*a2**2*a3)$

aan_(7) := (df(x2,t,2)*df(x4,t)*x2*x7*a1*a2*a3 + df(x2,t,2)*x2*x4*x7*a1*a2*a3*t1 - 2*df(x2,t)**2*df(x4,t)*x7*a1*a2*a3 - 
2*df(x2,t)**2*x4*x7*a1*a2*a3*t1 + 2*df(x2,t)*df(x4,t,2)*x2*x7*a1*a2*a3 + df(x2,t)*df(x4,t)*x2*x6*x7*a1**2*a2*a3 + df(x2,
t)*df(x4,t)*x2*x7*a1*a2*(a1*c6a + a3*c3a + a3*c5a + a3*i1a + 2*a3*t1) + df(x2,t)*x2*x4*x6*x7*a1**2*a2*a3*t1 + df(x2,t)*
x2*x4*x7*a1*a2*t1*(a1*c6a + a3*c3a + a3*c5a + a3*i1a) - df(x4,t,3)*x2**2*x7*a1*a2*a3 - df(x4,t,2)*x2**3*x7*a1*a2**2*a3 -
 df(x4,t,2)*x2**2*x6*x7*a1**2*a2*a3 - df(x4,t,2)*x2**2*x7*a1*a2*a3*(c3a + c5a + i1a + t1) + df(x4,t)**2*x2*x6*x7*a1**3*
a3 + df(x4,t)**2*x2*x7*a1**2*a2*c6a - df(x4,t)*x2**3*x7*a1*a2**2*a3*(c3a + t1) - df(x4,t)*x2**2*x5*x7*a1**2*a2*a3*t2 + 
df(x4,t)*x2**2*x6*x7*a1**2*a2*a3*( - c3a + i1 - t1) + df(x4,t)*x2**2*x7*a1*a2*(a1*c6a*k3 + a1*c6a*kdeg - a3*c3a*c5a - a3
*c3a*i1a - a3*c3a*t1 - a3*c5a*t1 + a3*e1a*i1a*kv - a3*i1a*t1) + 2*df(x4,t)*x2*x4*x6*x7*a1**3*a3*t1 + df(x4,t)*x2*x4*x7*
a1**2*a2*c6a*t1 - df(x4,t)*x2*x5*x7*a1**2*a2*c6a*t2 + df(x7,t)*x2**3*x7*a1*a2**2*a3*e1a + df(x7,t)*x2**3*a2**2*a3*e1a*( 
- c3a + e1a*kv) - x2**3*x4*x7*a1*a2**2*a3*c3a*t1 - x2**3*x6*x7*a1*a2**2*a3*e1a*i1*kv + x2**3*x6*a2**2*a3*e1a*i1*kv*(c3a 
- e1a*kv) + x2**3*x7**2*a1*a2**2*a3*c1a*c4a + x2**3*x7*a1*a2**2*a3*c2a*c4a - x2**2*x4*x5*x7*a1**2*a2*a3*t1*t2 + x2**2*x4
*x6*x7*a1**2*a2*a3*t1*( - c3a + i1) + x2**2*x4*x7*a1*a2*t1*(a1*c6a*k3 + a1*c6a*kdeg - a3*c3a*c5a - a3*c3a*i1a + a3*e1a*
i1a*kv) + x2*x4**2*x6*x7*a1**3*a3*t1**2 - x2*x4*x5*x7*a1**2*a2*c6a*t1*t2)/(a1*a2*a3)$

RANDOMLY CHOSEN NUMERICAL PARAMETER(S) VECTOR$

b2_ := {t1=63,
a3=41,
t2=38,
c1a=35,
c2a=32,
c3a=23,
c4a=8,
c5a=10,
c6a=49,
c1=25,
c2=22,
c3=45,
c4=12,
c5=50,
k1=20,
k3=46,
kprod=4,
kdeg=19,
kv=42,
i1=62,
e2a=58,
i1a=60,
e1a=47,
c1c=28,
c2c=11,
c3c=26}$

EXHAUSTIVE SUMMARY $

flist_ := {kprod - 4,
 - k3 + 46,
kdeg - 19,
 - t2 + 38,
 - t2 + 38,
a3 - 41,
t1 - 63,
t2 - 38,
 - t1 + 63,
 - t2 + 38,
a3 - 41,
 - c3a*t1 + 1449,
 - a3*e2a*kv*t2 + 3795288,
c1a*c4a - 280,
 - c6a*e2a*kv*t2 + 4535832,
t1 - 63,
 - t1 + 63,
t1 - 63,
 - t2 + 38,
e1a - 47,
 - c3a*t1 + 1449,
 - e1a*i1*kv + 122388,
c1a*c4a - 280,
c2a*c4a - 256,
 - t1*t2 + 2394,
k3 + kdeg - 65,
(a1*(49*a3 - 41*c6a))/(49*a2*c6a),
(49*a3 - 41*c6a)/(49*c6a),
( - 49*a3*t2 + 1558*c6a)/(49*c6a),
(49*a3*i1 - 2542*c6a)/(49*c6a),
k3 + kdeg - 65,
(a1*( - a3 + 41))/a2,
c6a*t2 + e2a*kv*t2 - 94430,
(a1*(7*a3*t1 - 369*c6a))/(7*a2*c6a),
(a1*(a3 - 41))/a2,
(a1**2*( - a3 + 41))/a2**2,
 - k1 - kdeg + 39,
(a1*(49*a3 - 41*c6a))/(49*a2*c6a),
(a1*(1519*a3*i1a - 1230*c6a*i1))/(1519*a2*c6a*i1),
(127596*a3 - 41*c6a*i1*kv)/(127596*c6a*i1*kv),
( - 127596*a3 + 41*c6a*i1*kv)/(127596*c6a*i1*kv),
(49*a3 - 41*c6a)/(49*c6a),
(a1*(127596*a3 - 41*c6a*i1*kv))/(127596*c6a*i1*kv),
(a1*( - 49*a3 + 41*c6a))/(41*a3),
 - c3a - t1 + 86,
(3038*a3*e1a - 1927*c6a*i1)/(3038*c6a*i1),
(2*a1**2*(t1 - 63))/a2,
(a1*( - 3087*a3 + 41*c6a*t1))/(41*a3),
(a1*(1862*a3 - 41*c6a*t2))/(41*a3),
( - 49*a3*t2 + 1558*c6a)/(49*c6a),
(a1*( - 49*a3 + 41*c6a))/(49*c6a),
k3 + kdeg - 65,
(a1*(7*a3*t1 - 369*c6a))/(7*a2*c6a),
(a1*(217*a3*i1a*t1 - 11070*c6a*i1))/(217*a2*c6a*i1),
 - c3a*t1 + i1*t1 - 2457,
(a1**2*(t1**2 - 3969))/a2,
(a1*(117306*a3 - 41*c6a*t1*t2))/(41*a3),
 - a2*e1a*i1*kv + 122388*a2 + a3*e2a*i1*kv - 6192312,
(a2*(c3a*e1a*i1*kv - e1a**2*i1*kv**2 + 238778988))/a1,
 - c6a - e2a*kv - 2*t2 + 2561,
c6a + e2a*kv + t2 - 2523,
 - c3a + i1 - t1 + 24,
(a2*( - c3a*e1a + e1a**2*kv - 91697))/a1,
( - 49*a3*e1a*kv + 49*a3*i1 + 78392*c6a)/(49*c6a),
(a2*(c3a*e1a*i1*kv - e1a**2*i1*kv**2 + 238778988))/a1,
a2*e1a - 47*a2 - a3*e2a + 2378,
(a2*( - c3a*e1a + e1a**2*kv - 91697))/a1,
 - c3a - c5a - i1a - t1 + 156,
a3*e2a*kv + a3*k3 + a3*kdeg + a3*t2 - 104099,
(a1*( - 49*a2*a3 + 41*a2*c6a - 41*a3*c6a + 2009*a3))/(41*a2*a3),
(a1*( - a2*t1*t2 + 2394*a2 + a3*t1*t2 - 98154))/a2,
(a1*( - a2*t2 + 38*a2 + a3*t2 - 1558))/a2,
(2*a1**2*(a2*t1 - 63*a2 - a3*t1 + 2583))/a2**2,
(a1*( - 3087*a2*a3 + 41*a2*c6a*t1 - 41*a3*c6a*t1 + 126567*a3))/(41*a2*a3),
(a1*(1862*a2*a3 - 41*a2*c6a*t2 + 41*a3*c6a*t2 - 76342*a3))/(41*a2*a3),
(a1**2*(a2*t1**2 - 3969*a2 - a3*t1**2 + 162729))/a2**2,
(a1*(117306*a2*a3 - 41*a2*c6a*t1*t2 + 41*a3*c6a*t1*t2 - 4809546*a3))/(41*a2*a3),
a3*e2a*k3*kv + a3*e2a*kdeg*kv + a3*k3*t2 + a3*kdeg*t2 - 6593210,
(a1*(a2*t1 - 63*a2 - a3*t1 + 2583))/a2,
 - a2*c3a - a2*t1 + 86*a2 + a3*e2a*kv + a3*t2 - 101434,
c6a*e2a*kv + c6a*t2 + e2a*kv*t2 - k3*t2 - kdeg*t2 - 211324,
(a1*( - a2*c3a*t1 + a2*i1*t1 - 2457*a2 + a3*e2a*kv*t1 - a3*i1*t1 - 6132042))/a2,
( - 3087*a1*a3 + 41*a1*c6a*t1 + 41*a3*c3a*t1 + 41*a3*c5a*t1 + 41*a3*i1a*t1 - 240219*a3)/(41*a3),
( - 49*a1*a3 + 41*a1*c6a + 41*a3*c3a + 41*a3*c5a + 41*a3*i1a + 82*a3*t1 - 8979*a3)/(41*a3),
 - c3a - c5a + c6a + e2a*kv - i1a - t1 + t2 - 2367,
( - 200655*a1*a3 + 41*a1*c6a*k3*t1 + 41*a1*c6a*kdeg*t1 - 41*a3*c3a*c5a*t1 - 41*a3*c3a*i1a*t1 + 41*a3*e1a*i1a*kv*t1 - 
301771890*a3)/(41*a3),
(a1*( - a2*c3a + a2*i1 - a2*t1 + 24*a2 + a3*e2a*kv - a3*i1 + a3*t1 - 99917))/a2,
( - 3087*a1*a2*a3 + 41*a1*a2*c6a*t1 - 41*a1*a3*c6a*t1 + 126567*a1*a3 + 41*a2*a3*c3a*t1 + 41*a2*a3*c5a*t1 + 41*a2*a3*i1a*
t1 - 240219*a2*a3)/(41*a2*a3),
a2*c2a*c4a - 256*a2 + c6a*e2a*k3*kv + c6a*e2a*kdeg*kv + c6a*k3*t2 + c6a*kdeg*t2 + e2a*k3*kv*t2 + e2a*kdeg*kv*t2 - 
13896610,
( - 3185*a1*a3 + 41*a1*c6a*k3 + 41*a1*c6a*kdeg - 41*a3*c3a*c5a - 41*a3*c3a*i1a - 41*a3*c3a*t1 - 41*a3*c5a*t1 + 41*a3*e1a
*i1a*kv - 41*a3*i1a*t1 - 4549811*a3)/(41*a3),
( - 200655*a1*a2*a3 + 41*a1*a2*c6a*k3*t1 + 41*a1*a2*c6a*kdeg*t1 - 41*a1*a3*c6a*k3*t1 - 41*a1*a3*c6a*kdeg*t1 + 8226855*a1
*a3 - 41*a2*a3*c3a*c5a*t1 - 41*a2*a3*c3a*i1a*t1 + 41*a2*a3*e1a*i1a*kv*t1 - 301771890*a2*a3)/(41*a2*a3),
( - 49*a1*a2*a3 + 41*a1*a2*c6a - 41*a1*a3*c6a + 2009*a1*a3 + 41*a2*a3*c3a + 41*a2*a3*c5a - 41*a2*a3*c6a - 41*a2*a3*e2a*
kv + 41*a2*a3*i1a + 82*a2*a3*t1 - 82*a2*a3*t2 + 96022*a2*a3)/(41*a2*a3),
( - 3185*a1*a2*a3 + 41*a1*a2*c6a*k3 + 41*a1*a2*c6a*kdeg - 41*a1*a3*c6a*k3 - 41*a1*a3*c6a*kdeg + 130585*a1*a3 - 41*a2*a3*
c3a*c5a - 41*a2*a3*c3a*i1a - 41*a2*a3*c3a*t1 - 41*a2*a3*c5a*t1 + 41*a2*a3*c6a*e2a*kv + 41*a2*a3*c6a*t2 + 41*a2*a3*e1a*
i1a*kv + 41*a2*a3*e2a*kv*t2 - 41*a2*a3*i1a*t1 - 13315365*a2*a3)/(41*a2*a3)}$

MODEL PARAMETER SOLUTION(S)$

 G_:=GROESOLVE(FLIST_,B1_) $

g_ := {{c2a=256/c4a,
k1=20,
kprod=4,
c1a=280/c4a,
c5a=10,
i1a=60,
e1a=47,
c3a=23,
i1=62,
k3=46,
kdeg=19,
e2a=58,
t1=63,
kv=42,
c6a=49,
a3=41,
t2=38}}$

MODEL NON IDENTIFIABLE$
Elapsed time for NFkB: 12.1259494 seconds
