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