import sympy as sp from IPython.display import display v,a,t = sp.symbols('v a t',real=True) x = sp.Function('x')(t) eq1 = sp.Eq(sp.diff(x,t,2),a) display(eq1) display(sp.dsolve(eq1)) display(sp.dsolve(eq1,ics={x.subs(t,0):0, sp.diff(x,t,1).subs(t,0):v}))