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