syms y(t)
Dy = diff(y);
y(t) = dsolve(diff(y, t, 2) == - 5*Dy - 3*y + 6*t^2 - t - 1, y(0) == 10, Dy(0) == -7)