x = var('x')
y = function('y', x)

DE = diff(y,x,x) + 3*diff(y,x) + 2*y

ans = desolve(DE, [y,x])
print ans