dsolve({diff(x(t),t)-diff(y(t),t)+x(t) = 2*t, diff(diff(x(t),t),t)+diff(y(t),t)-9*x(t)+3*y(t) = sin(2*t)})