syms x(t) Dx = diff(x); dsolve(diff(Dx) == exp(-x^2), Dx(0) == 1)