dx1/dt = -delta,NE*x1*x4 - mu,N*x1; dx2/dt = -mu,LE*x2*x5 + delta,EL*x2 - x2^2*mu,LL - x2*delta,LM; dx3/dt = x2*delta,LM - x3*mu,M; dx4/dt = rho,P*x4^2 - x2*x4*mu,PL - x5*mu,PE*x4 - x4*mu,P; dx5/dt = -mu,EE*x5^2 + delta,NE*x1*x4 - delta,EL*x5 + x5*x4*rho,E; y1 = x1; y2 = x3; y3 = x2 + x5