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