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