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