x = var('x')

#y = function('y',x)

#diff(v,x) = diff(y,x) + 4

v = function('v', x)

DE = diff(v,x) - 4 - v^2

desolve(DE, [v,x])