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