dx1/dt = -b1*b0*x1*x5*x3 - b0*x1*x3 - x1*mu + mu + g*x4; dx2/dt = M*x5; dx3/dt = -nu*x3 + b1*b0*x1*x5*x3 + b0*x1*x3 - mu*x3; dx4/dt = nu*x3 - mu*x4 - g*x4; dx5/dt = -M*x2; y1 = x3; y2 = x4