dsolve({a1*diff(diff(x(t),t),t)+b1*diff(x(t),t)+c1*x(t)-A*diff(y(t),t) = B*exp(I*omega*t), a2*diff(diff(y(t),t),t)+b2*diff(y(t),t)+c2*y(t)+A*diff(x(t),t) = 0})