dx1/dt = alpha*x3 - beta*x1; dx2/dt = (gama*sigma*x1*x2 - delta*sigma*x4*x2) / (x4); dx3/dt = (-b*c*x3 - b*x3*x2 + 1) / (c + x2); dx4/dt = gama*x1 - delta*x4; y1 = x3