dsolve({diff(x(t),t) = (a*y(t)+b)*x(t), diff(y(t),t) = (c*x(t)+d)*y(t)})