y'=-y^2/x + xy/(x+1) + x/(x+1)