import sympy as sp
from IPython.display import display
x = sp.Symbol('x')
f = sp.Function('f')(x)
eq1 = sp.diff(f,x,2)-2*sp.diff(f,x)+2*f
display(eq1)
display(sp.dsolve(eq1))
display(sp.dsolve(eq1,ics={f.subs(x,0):0,
        sp.diff(f,x,1).subs(x,0):1,sp.diff(f,x,2).subs(x,0):2}))