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