dx1/dt = reaction,1,k1*x1 - reaction,9,k1*x1 - reaction,1,k2*x1; dx2/dt = -reaction,1,k1*x1 + reaction,1,k2*x1 - x2*EGFR,turnover + EGFR,turnover*u1; dx3/dt = -reaction,4,k1*x3 + reaction,9,k1*x1 - x3*x9*reaction,2,k1 + reaction,3,k1*x5 + x5*reaction,2,k2; dx4/dt = x6*reaction,5,k1*x7 - reaction,6,k1*x4 - x4*reaction,5,k2; dx5/dt = x3*x9*reaction,2,k1 - reaction,3,k1*x5 - x5*reaction,2,k2; dx6/dt = -x6*reaction,7,k1 - x6*reaction,5,k1*x7 + reaction,6,k1*x4 + reaction,3,k1*x5 + x4*reaction,5,k2; dx7/dt = x8*reaction,8,k1 - x6*reaction,5,k1*x7 + x4*reaction,5,k2; dx8/dt = -x8*reaction,8,k1 + reaction,6,k1*x4; dx9/dt = x6*reaction,7,k1 - x3*x9*reaction,2,k1 + x5*reaction,2,k2; y1 = x3*a1 + a1*x5; y2 = a2*x6 + a2*x4; y3 = x8*a3
