dx1/dt = k*x5 - x1*u; dx2/dt = lm - x2*d - x2*x1*beta; dx3/dt = c*x4*q*x5 - h*x3; dx4/dt = -b*x4 + c*x4*x2*x5 - c*x4*q*x5; dx5/dt = x2*x1*beta - a*x5; y1 = x3; y2 = x4