dsolve({diff(x(t),t) = a1*x(t)+b1*y(t)+c1, diff(y(t),t) = a2*x(t)+b2*y(t)+c2})