apoptosis_1=function(t,y,parms) { # # Transfer dependent variable vector to problem variables y_hif =y[1]; y_o2 =y[2]; y_p300=y[3]; y_p53 =y[4]; y_casp=y[5]; y_kp =y[6]; # # Apoptosis model ODEs yt_hif =a_hif-a_3*y_o2*y_hif-a_4*y_hif*y_p300-a_7*y_p53*y_hif; yt_o2 =a_o2 -a_3*y_o2*y_hif+a_4*y_hif*y_p300-a_11*y_o2; yt_p300=a_8 -a_4*y_hif*y_p300-a_5*y_p300*y_p53; yt_p53 =a_p53-a_9*y_p53 -a_5*y_p300*y_p53; yt_casp=a_12 +a_9*y_p53-a_13*y_casp; yt_kp =-a_10*y_casp*y_kp-a_14*y_kp +a_11*y_o2; # # Increment calls to apoptosis_1 ncall <<- ncall+1 # # Return derivative vector return(list(c(yt_hif,yt_o2,yt_p300,yt_p53,yt_casp,yt_kp))) }