x = var('x')   
y = function('y',x)   
sol=desolve(diff(y,x) - y - 4, [y,x])
sol