syms y(t) Dy = diff(y); y(t) = dsolve(diff(y, t, 2) == -10 - Dy + 2*y, y(0) == 12, Dy(0) == -2)