WRITE "NFkB"$

% add ICs to B_??!!

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

FOR EACH EL_ IN B_ DO DEPEND EL_,T$

%B1_ IS THE UNKNOWN PARAMETER 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}$

%NUMBER OF STATES 
NX_:=15$
%NUMBER OF INPUTS 
NU_:=1$
%NUMBER OF OUTPUTS 
NY_:=6$

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

% INITIAL CONDITIONS
%IC_:={x1=x01,x2=0,x3=0,x4=0,x5=0,x6=x06,x7=x07,x8=x08,x9=x09,x10=x010,x11=x011,x12=x012,x13=NF,x14=x014,x15=x015}$
%CONDINIZ()$
END$
