dsolve({a*diff(x(t),t)+b*diff(y(t),t) = alpha*x(t)+beta*y(t), b*diff(x(t),t)-a*diff(y(t),t) = beta*x(t)-alpha*y(t)})