syms y(t) Dy = diff(y); y(t) = dsolve(diff(y, t, 2) == 27 - 6*Dy - 9*y, y(0) == 5, Dy(0) == -5)