Modele: cinetique du premier ordre
step=0.0001;
vk1 = k1*MEK*I_MEK;
vk2 = k2*MEK_i;
vk3 = k3*ERK*MEK;
vk4 = k4*ERK_p*DUSP_a;
vk5 = k5*DUSP_i*ERK_p;
vk6 = k6*DUSP_a;
vk7 = k7*ERK*I_ERK;
vk8 = k8*ERK_i;
vk9 = k9*ERK_p*I_ERK_p;
vk10 = k10*ERK_p_i;
MEK = MEK+step*(vk2-vk1);
I_MEK = I_MEK+step*(vk2-vk1);
MEK_i = MEK_i+step*(vk1-vk2);
ERK = ERK+step*(vk4+vk8-vk3-vk7);
I_ERK = I_ERK+step*(vk8-vk7);
ERK_i = ERK_i+step*(vk7-vk8);
ERK_p = ERK_p+step*(vk3+vk10-vk4-vk9);
I_ERK_p = I_ERK_p+step*(vk10-vk9);
ERK_p_i = ERK_p_i+step*(vk9-vk10);
DUSP_i = DUSP_i+step*(vk6-vk5);
DUSP_a = DUSP_a+step*(vk5-vk6);
MEK
I_MEK
MEK_i
ERK
I_ERK
ERK_i
ERK_p
I_ERK_p
ERK_p_i
DUSP_i
DUSP_a