reset()
x = var('x')
y = function('y', x)
DE = desolve(diff(y,x)-1-0.001*y^2 , y)
DE.simplify()