syms y(t) Dy = diff(y); y(t) = dsolve(diff(y, t, 2) == 2*exp(-4*t) - 3*Dy + 4*y, y(0) == 2, Dy(0) == -3.4)