dx1/dt = a*x1 + b*x1 - x2*c*x1; dx2/dt = -a*b*x2 + d*x2*x1; y1 = x1