dsolve({3*diff(x(t),t)+7*diff(y(t),t)+8*x(t)+24*y(t) = exp(2*t), 4*diff(x(t),t)+9*diff(y(t),t)+11*x(t)+31*y(t) = exp(t)})