WRITE "NGF_Erk"$

% B_ IS THE VARIABLE VECTOR 
B_:={k_3__TrkA_NGF,k_6__RasGTP,k_8__pRaf,k_10__pMek,s__pErk,y1}$

FOR EACH EL_ IN B_ DO DEPEND EL_,T$

%B1_ IS THE UNKNOWN PARAMETER VECTOR
B1_:={k_1,k_2,k_4,k_5,k_7,k_9,k_11,k_3__TrkA_0,k_6__Ras_0,k_8__Raf_0,k_10__Mek_0,s__Erk_0,s_K,n,NGF_0}$

%NUMBER OF STATES 
NX_:=5$
%NUMBER OF INPUTS 
NU_:=0$
%NUMBER OF OUTPUTS 
NY_:=1$

%MODEL EQUATIONS
C_:={df(k_3__TrkA_NGF,t)=-k_3__TrkA_NGF*k_2-NGF_0*k_1*(k_3__TrkA_NGF-k_3__TrkA_0),
     df(k_6__RasGTP,t)=-k_6__RasGTP*k_5-(k_6__RasGTP-k_6__Ras_0)*(k_4+(k_3__TrkA_NGF*s_K^n)/(s__pErk^n+s_K^n)),
     df(k_8__pRaf,t)=k_6__RasGTP*(k_8__Raf_0-k_8__pRaf)-k_7*k_8__pRaf,
     df(k_10__pMek,t)=k_8__pRaf*(k_10__Mek_0-k_10__pMek)-k_9*k_10__pMek,
     df(s__pErk,t)=-k_11*s__pErk-k_10__pMek*(s__pErk-s__Erk_0),
     y1=s__pErk}$

SEED_:=25$
DAISY()$

% INITIAL CONDITIONS 
IC_:={0,
     (k_4*k_6__Ras_0)/(k_4 + k_5),
     (k_4*k_6__Ras_0*k_8__Raf_0)/((k_4 + k_5)*(k_7 + (k_4*k_6__Ras_0)/(k_4 + k_5))),
     (k_4*k_6__Ras_0*k_8__Raf_0*k_10__Mek_0)/((k_4 + k_5)*(k_7 + (k_4*k_6__Ras_0)/(k_4 + k_5))*(k_9 + (k_4*k_6__Ras_0*k_8__Raf_0)/((k_4 + k_5)*(k_7 + (k_4*k_6__Ras_0)/(k_4 + k_5))))),
     (k_4*k_6__Ras_0*k_8__Raf_0*k_10__Mek_0*s__Erk_0)/((k_4 + k_5)*(k_7 + (k_4*k_6__Ras_0)/(k_4 + k_5))*(k_11 + (k_4*k_6__Ras_0*k_8__Raf_0*k_10__Mek_0)/((k_4 + k_5)*(k_7 + (k_4*k_6__Ras_0)/(k_4 + k_5))*(k_9 + (k_4*k_6__Ras_0*k_8__Raf_0)/((k_4 + k_5)*(k_7 + (k_4*k_6__Ras_0)/(k_4 + k_5))))))*(k_9 + (k_4*k_6__Ras_0*k_8__Raf_0)/((k_4 + k_5)*(k_7 + (k_4*k_6__Ras_0)/(k_4 + k_5)))))}$
CONDINIZ()$
END$



 

