dsolve({diff(x(t),t) = x(t)*(a*(p*x(t)+q*y(t))+alpha), diff(y(t),t) = y(t)*(beta+b*(p*x(t)+q*y(t)))})